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ABSTRACT 


VENABLE,  CHARLES  J.  Scheduling  of  Multiproducts  with 
Limited  Resources  in  an  UPS  Testing  Facility.  (Under  the 
direction  of  DR.  S.  E.  ELMAGHRABY.) 

Scheduling  procedures  in  an  Uninterruptible  Power 
Supply  (UPS)  testing  facility  were  investigated.  It  was 
determined  that  the  facility  was  experiencing  an  inordinate 
number  of  missed  product  due  dates  and  it  was  the  objective 
of  this  study  to: 

1.  Develop  a  model  to  describe  the  assignment 
functions  with  regards  to  resource  constraints. 

2.  Develop  and  implement  an  algorithm  to  determine  the 
practicability  of  monthly  test  requirements. 

3.  Determine  feasibility  of  a  capacity  loading  model. 

4.  Determine  feasibility  of  an  optimum  scheduling 
model . 

An  analytical  scheduling  model  was  developed.  However, 
due  to  a  unique  shift  organization  the  computational 
requirements  of  the  model  overshadowed  its  usefulness.  In 
order  to  lessen  the  number  of  required  constraints  a  reduced 
formulation  was  developed  which  considered  only  one  shift. 

An  algorithm  was  then  presented  which,  when  used  with  the 
reduced  formulation,  produced  acceptable  feasible  scheduling 
solutions  for  both  shifts. 

It  was  also  determined  that  due  to  the  multiproduct 
nature  of  the  industry,  test  facility  capacity  was  not  con¬ 
sidered  for  specific  product  mixes.  An  analytical  capacity 


model  was  developed  to  be  used  in  conjunction  with  a 
heuristic  to  estimate  capacity  prior  to  schedule 
development.  Again,  the  computational  requirements  of  the 
model  proved  too  large,  and  the  determination  was  made  to 
use  a  simulation  approach.  An  operational  tool  for  the  test 
facility  manager  was  developed  in  the  form  of  a  sequential 
simulation,  to  be  utilized  as  an  aid  in  scheduling  and 
determining  capacity  and  resource  restrictions. 

Data  was  collected  from  the  test  facility  in  order  to 
develop  valid  parameters  to  the  simulation.  Combining  the 
input  data  with  the  sequential  simulation  model  a  valid  ran¬ 
domized  simulator  for  the  test  facility  was  constructed.  A 
simulation  experiment  was  conducted  to  determine  the  pre¬ 
ferred  queue  discipline  and  second  shift  processing  cri¬ 
terion.  It  was  statistically  determined  that  neither  had  an 
effect  on  the  number  of  late  jobs.  However,  it  was  shown 
that  the  dispatch  rule  of  shortest  processing  time  signifi¬ 
cantly  increased  the  mean  lateness,  while  the  second  shift 
processing  criterion  had  no  effect.  The  second  shift  pro¬ 
cessing  criterion  was  also  shown  to  have  no  effect  on  second 
shift  testor  utilization. 

These  results,  thereby,  indicated  that  the  source  of 
the  inordinate  number  of  late  jobs  lies  outside  the  testing 
facility  and  that  the  second  shift  testor  configuration  may 
not  meet  the  future  needs  of  the  system.  This  conclusion 
was  verified  by  the  sparsity  of  the  second  shift 
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"qualification  matrix",  and  a  recommendation  to  management 
was  made  to  this  effect. 


It  was  also  recommended  that  the  reduced  scheduling 
model  and  the  related  algorithm  be  employed  by  the  company 
both  as  a  planning  tool  and  an  operational  tool.  As  a 
planning  tool  it  was  recommended  for  use  as  a  capacity 
model,  with  the  assumption  that  products  are  immediately 
available,  and  as  the  first  step  in  a  backward  planning 
process.  As  an  operational  tool  it  was  recommended  for  use 
by  the  test  facility  manager  on  receipt  of  monthly  test 
requirements.  Additionally,  it  was  recommended  that  the 
test  facility  manager  employ  the  sequential  simulator  as  a 
desk-side  reference  to  provide  immediate  scheduling  update 
as  production  information  is  received. 
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I.  Introduction 


With  the  growing  reliance  in  business  and  government  on 
computers  and  computer-based  equipment  a  new  problem  has 
emerged  -  how  to  protect  this  equipment  against  power  fluc¬ 
tuations  and  outages.  There  are  several  technologies  avail¬ 
able  which  offer  varying  degrees  of  power  protection.  The 
most  complete  solution,  however,  is  an  Uninterruptible  Power 
Supply  (UPS) . 

The  UPS  is  composed  of  a  static  rectif ier/charger  which 
takes  the  incoming  AC  power  and  converts  it  to  DC;  a  static 
inverter  which  converts  DC  back  to  clean,  reliable  AC  power; 
and  a  bank  of  battery  cells  to  supply  full  operating  power, 
through  the  inverter,  for  a  specified  period  of  time  during 
an  extended  power  outage  (Fig.  1.1,  1.2).  Unlike  other 
types  of  standby  power  sources  which  only  take  over  during  a 
utility  failure,  the  UPS  system  actually  becomes  the  princi¬ 
pal  source  of  power  to  the  critical  load,  supplying  smooth, 
regulated  AC  power  regardless  of  the  condition  of  incoming 
utility  power. 

The  demand  for  UPS  protection  has  increased  dramati¬ 
cally  as  computer  technology  has  evolved  and  will  continue 
to  accelerate  in  the  coming  years.  In  order  to  meet  this 
demand  UPS  production  firms  are  continuing  to  introduce  new 
product  lines,  provide  more  custom  engineering  and  increase 
production.  However,  sales  continue  to  outpace  production 
capacity  causing  strained  customer  relations,  possible  market 
loss,  employee  morale  problems  due  to  increased  overtime 
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requirements  and  numerous  other  problems.  In  order  to  mini¬ 
mize  these  problems,  UPS  production  firms  are  increasingly 
applying  operations  research  methods  to  optimize  production 
methods . 

Common  to  all  such  firms  is  the  requirement  to  fully 
test  all  UPS  products,  with  a  battery,  following  final 
assembly.  This  final  testing  provides  a  natural  choke  point 
in  the  production  sequence.  This  area  then  becomes  critical 
in  any  production  optimization  study.  This  study  will, 
therefore,  concentrate  on  the  UPS  system  test  phase. 

In  order  to  properly  diagnose  the  problems  associated 
with  this  phase  of  the  production  sequence  a  generalized 
discussion  of  the  system  test  area  and  the  resource  allo¬ 
cation/planning  systems  in  use  will  follow. 

II .  Problem  Definition 

The  System  Test  Area 

The  system  test  area  must  be  capable  of  fully  testing 
all  UPS  Products,  with  a  battery.  this  requires  a  tremen¬ 
dous  power  versatility.  Products  range  in  size  from  15 
kilowatts  to  600  kilowatts.  They  are  available  in  single  or 
three-phase  with  input  voltages  ranging  from  208  to  480 
volts  and  output  voltages  ranging  from  120  to  480  volts. 
Additionally,  depending  on  product  destination,  it  may  be 
configured  for  50,  60,  or  415  Hz  per  second.  Appendix  A 
lists  72  product  types  and  as  mentioned  previously,  new  UPS 
products  are  continuing  to  be  introduced  to  the  market. 


Systems  can  also  be  obtained  in  both  single  module 
(non-redundant)  configurations. 

System  test  areas  are  divided  into  stations  with  all 
stations  having  four  variables;  maximum  size  unit,  load 
size,  available  voltages,  and  a  battery.  These  variables 
will  dictate  the  product  types  that  may  be  tested  at  that 
station.  Appendix  B  details  these  variables  for  a  given 
test  facility  and  gives  a  matrix  showing  which  of  the  72 
identified  products  can  be  tested  at  which  stations. 

This  test  facility  will  then  provide  a  practical 
example  for  the  discussion  of  the  planning  and  resource 
allocation  that  follows. 

The  System  Test  Forward  Planning  System 

Upon  sale  of  an  UPS  the  firm  gives  their  customer  a 
promise  date  for  delivery  of  the  equipment.  This  promise 
date  is  arrived  at  by  a  mathematical  equation  that  incorpo 
rates  manufacturing  lead  times.  Each  product  line  has  a 
standard  amount  of  time,  lead  time,  for  production.  These 
times  were  developed  by  the  firm's  Industrial  Engineering 
Department  for  each  phase  of  the  production  sequence, 
including  Systems  Test.  These  figures  were  then  totaled 
with  an  additional  "safety"  factor  to  obtain  the  standard. 
Due  dates  are  then  published  in  a  number  of  documents,  chi 
of  which  is  the  Project  Management  System  PM  21  report 
(Appendix  C)  which  would  be  used  by  the  manager  of  the 
Systems  Test  area  to  prepare  his  monthly  plan  for  test 
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(Appendix  H) .  The  PM  21  report  is  by  due  date  and  in  effect 
sets  the  priority  for  the  equipment.  In  order  to  gain  more 
information  on  the  product,  the  manager  would  then 
cross-reference  the  equipment  in  priority  with  the  Project 
Management  System  PM01  report  (Appendix  D) .  Of  primary 
concern  to  the  manager  at  this  time  is  the  date  at  which  the 
item  is  scheduled  to  arrive  at  Systems  Test  (Reference  #1, 
App.  D) ,  and  the  scheduled  completion  date  for  Systems  Test 
(Reference  #2,  App.  D) .  These  dates  give  the  manager  the 
window  within  which  the  equipment  must  be  tested.  The  next 
item  of  information  from  this  report,  of  interest  to  the 
manager,  is  the  type  of  equipment  (Reference  #3,  App.  D) 
which  can  also  be  obtained  from  the  PM21  report,  Appendix  C. 
With  this  information  he  can  schedule  a  test  station 
compatible  with  this  particular  model.  The  manager  must 
then  determine  the  number  of  test  hours  required  for  the 
equipment . 

Each  product  has  a  standard  allotted  test  time.  How¬ 
ever,  additional  hours  may  be  necessary.  For  example,  the 
customer  may  wish  to  observe  or  personally  verify  testing. 
This  is  termed  a  "witness  test."  Or  the  customer  may  have 
ordered  additional  specifications  to  the  product  requiring 
custom  engineering  and,  therefore,  additional  testing. 
Referring  again  to  the  PM01  report  the  manager  determines  if 
a  witness  test  is  required  (Reference  #4,  App.  D)  and  if  its 
order  type  is  custom  or  standard  (Reference  #5,  App.  D) . 
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To  gain  more  technical  information  about  the  piece  of 
equipment,  specifically  if  other  tests  are  required  which 
will  add  test  time,  the  manager  refers  to  the 
Special/Standard  Order  Data  Sheet,  SSPODS  (Appendix  E) . 

Using  the  allotted  test  times  (Appendix  F)  and  the  known 
test  station  constraints  (Appendix  B)  the  manager  schedules 
the  UPS  for  a  particular  date  and  station.  Once  the 
manager's  test  plan  has  been  prepared  the  supervisor  of 
Systems  Test  will  schedule  his  personnel  based  on  their 
productivity  hours  for  the  month  and  their  skill 
qualifications . 

As  the  UPS  moves  through  the  assembly  process,  updates 
will  be  made  to  the  PM01  report  on  its  progress.  Of  primary 
interest  to  the  manager  is  the  Current  Scheduled  Start 
(Reference  #6,  App.  D)  and  the  Current  Scheduled  Completion 
(Reference  #7,  App.  D) .  These  dates  will  advise  the  manager 
as  to  any  changes  in  the  availability  window  for  Systems 
Test  and  allow  him  to  revise  his  plan  of  tests  accordingly. 
The  PM01  report  is  run  each  week  and,  therefore,  provides 
the  manager  with  a  weekly  update.  This  section  is 
summarized  by  Figure  2.1. 


The  manager  of  Systems  Test  currently  uses  the  Monthly 
Sales  Plan  (Appendix  G)  as  his  priority  document.  From  this 
document  he  is  able  to  tell  which  week  of  the  month  the  com¬ 


pany  expects  the  equipment  to  be  delivered.  The  Monthly 
Sales  Plan  also  allows  the  manager  to  identify  sold,  defi¬ 
nite  orders,  versus  pending,  projected  orders.  The  manager 
then  cross-ref erences  the  system  number  for  the  particular 
piece  of  equipment  with  the  system  number  on  the  PM01  report 
(Appendix  D)  and  determines  the  Operations  Original  Promise 
Date  (Reference  #8,  App.  D) .  This  date  will  further  set  his 
priority  for  test  as  this  is  the  last  possible  date  for  com¬ 
pletion  of  production.  Additionally,  he  is  able  to  tell 
from  both  of  these  documents  whether  the  UPS  requires  a  wit¬ 
ness  test.  He  then  refers  to  the  SSPODS  (Appendix  E)  to 
determine  the  date  for  the  witness  test,  if  required.  This 
date  then  becomes  critical  in  the  test  sequence,  since  a 
representative  of  the  customer  will  be  present  on  that  date 
to  verify  testing.  Additionally,  the  planning  office  gives 
a  periodic  priority  update  in  informal  conferences  with  the 
System  Test  Supervisor.  In  this  conference  information  is 
collected  as  to  the  status  and  projected  completion  dates 
for  priority  products.  These  individuals  also  serve  as 
troubleshooters  if  a  problem,  such  as  engineer  support, 
stands  in  the  way  of  test  completion.  A  recent  change  in 
the  working  system  that  has  taken  place  since  the  beginning 
of  this  study  is  the  implementation  of  a  thrice  weekly 
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conference  on  the  current  status  of  production.  In  atten¬ 
dance  at  this  meeting  are  the  major  manufacturing  managers, 
including  the  Systems  Test  manager.  The  result  of  the  con¬ 
ference  is  an  updated  UPS  Production  Plan,  Appendix  I.  This 
plan  tracks  the  planned  and  actual  completion  dates  at  each 
stage  of  the  manufacturing  process  for  each  scheduled 
product.  This  plan  then  serves  to  update  the  current 
planned  arrival  date  of  a  product  to  Systems  Test.  This 
section  is  summarized  by  Figure  2.2. 
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Monthly 
Sales  Plan 

\ _ / 


PMO 1 


Sold  Products  by  Week 


Due  Date  (OPCPR) 
Product  Type 
Witness  Test 
Custom  or  Standard 


SSPODS 


\  UPS  / 
Production 
^  Plan  j 


Tests 

Witness  Test  Date  (If 
Applicable) 


Figure  2.2  The  Working  System 


Once  an  UPS  is  ready  for  Systems  Test  from  the  produc¬ 
tion  line,  the  "lead  hand"  notifies  the  Systems  Test 
Supervisor.  The  supervisor  determines  whether  a  station  is 
available  that  will  handle  the  specific  model.  If  so  it  is 
moved  to  the  station.  System*.  Test  does  not  have  a 
designated  holding  area.  Therefore,  if  space  is  not 
available,  it  may  sit  at  the  end  of  the  production  line,  be 
moved  to  a  floor  space  in  Systems  Test,  or  if  it  is  for 
stock,  moved  to  the  warehouse.  The  supervisor  then 
determines  whether  he  has  a  testor  available  to  begin 
testing  the  equipment.  Testors  are  only  qualified  to  test 
certain  product  lines  and  are  qualified  to  test  all  models 
within  that  product  line.  Additionally,  the  supervisor  must 
determine  whether  to  also  assign  the  UPS  to  second  shift 
based  on  its  priority  and  projected  system  test  time.  This 
decision  may  be  made  anytime  until  the  UPS  is  fully  tested. 
However,  once  the  decision  has  been  made  to  test  the  product 
on  second  shift  as  well  as  first,  it  will  continue  to  be 
tested  by  both  shifts  until  completion.  Similarly,  the 
individual  testors,  once  assigned  to  a  specific  UPS,  will 
stay  with  that  piece  of  equipment  until  completion.  This 
section  is  summarized  by  Figure  2.3. 

When  a  testor  receives  his  assignment  he  reports  to  the 
quality  control  office.  Here  he  is  issued  the  test  packet 
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for  the  particular  UPS  and  the  Test  Data  Sheets.  The  data 
sheets  provide  the  testor  with  the  test  sequence  and  serve 
as  a  final  record  of  testing  when  completed.  As  the  testor 
moves  through  the  test  sequence  he  completes  the  Unit  Status 
Report.  This  report  is  a  record  of  his  progress  and  serves 
as  an  inter-shift  report  if  the  UPS  is  assigned  to  second 
shift  or  if  for  some  unforeseen  reason  the  testor  is  not 
able  to  complete  the  sequence.  Additionally,  the  tester 
records  all  problems/ faults  that  he  encounters  in  the  test 
sequence  on  a  Test  Failure  Report.  This  is  a  quality 
control  form  and  gives  a  description  of  the  failure, 
analysis  of  the  cause,  and  the  action  taken. 

Each  day  the  supervisor  prepares  a  status  report.  This 
report  is  compiled  based  on  the  Unit  Status  Reports  and 
gives  a  percentage  of  completion,  the  current  procedure 
being  tested,  any  any  major  problems.  This  report  is 
forwarded  to  the  Quality  Assurance  Manager.  Once  an  UPS  has 
completed  test  the  Shipping  Department  is  informed  and  the 
UPS  is  moved. 

Existing  System  Problems 

The  manager  of  Systems  Test  does  not  have  confidence  in 
the  outlined  forward  planning  system.  This  is  due  primarily 
to  the  failure  of  the  update  procedure.  In  the  formal  system 
without  the  update  information,  initial  data  becomes  meaning¬ 
less  and  the  manager  is  unable  to  form  a  complete  plan  for 
test.  This  situation  has  been  recognized  and  an  attempt  to 
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rectify  the  problem  has  been  made  by  the  implementation  of 
the  thrice  weekly  conference  on  the  status  of  production. 
Through  this  conference  the  manager  is  kept  abreast  of 
production  dates  and  specif ically ,  scheduled  dates  for 
arrival  at  Systems  Test.  However,  delays  that  occur  in  the 
production  system,  although  now  accounted  for  by  the  UPS 
Production  Plan  (Appendix  I) ,  are  not  reflected  in  the 
original  and  current  dates  of  the  PM01  report.  Since 
Systems  Test  is  the  last  in  the  production  sequence,  the 
manager  finds  the  ideal  lead  times  drastically  cut  forcing 
him  to  make  up  production  lost  time  to  meet  the  promise 
date . 

Monthly  planning  is  based  primarily  on  sales  projec¬ 
tions,  and  even  though  standard  production  lead  times  are 
used  in  determining  due  dates,  consideration  is  not  given  to 
whether  the  specific  monthly  product  mix  will  be  able  to  be 
processed  by  Systems  Test  and  still  fall  within  the  desig¬ 
nated  lead  times,  meeting  the  due  date.  As  a  result, 

Systems  Test  finds  itself  overcommitted.  This  problem  is 
then  amplified  by  the  carryover  from  the  previous  month  of 
the  products  which  were  not  tested.  This  forces  a  severe 
backlog  to  occur  at  the  end  of  the  month  as  the  manager 
attempts  to  meet  the  monthly  projection  figures 
necessitating  numerous  overtime  hours. 
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III .  Project  Objectives 

It  has  become  apparent  during  the  Problem  Definition 
phase  of  this  study  that  two  major  problems  exist.  The  first 
is  the  lack  of  update  in  the  firm's  reporting  procedures. 

As  mentioned  previously,  the  weekly  manager's  conferences  on 
the  status  of  production  is  an  attempt  to  address  the  problem, 
however,  the  results  must  be  integrated  into  the  existing 
formal  planning  framework  to  gain  maximum  benefit.  The  second 
is  the  fact  that  Systems  Test  area  constraints  are  not  con¬ 
sidered  during  the  firm's  planning  phase.  This  is  a  result 
of  a  failure  to  accurately  determine  the  Systems  Test  area 
capacity.  This  capacity,  however,  varies  with  product  mix 
and,  therefore,  may  not  be  able  to  be  explicitly  determined. 

It  is  the  second  problem  on  which  the  remainder  of  this  study 
will  focus. 

It  is  the  objective  of  this  study  to: 

1.  Develop  a  model  to  describe  the  System  Test 
assignment  functions  with  regard  to  resource 
constraints . 

2.  Develop  and  implement  an  algorithm  which  will 
allow  the  manager  of  a  Systems  Test  facility  to 
determine  whether  he  will  be  able  to  meet  monthly  test 
requirements  as  delineated  in  the  firm's  monthly 
plan/ forecast . 

3.  Determine  feasibility  of  a  capacity  loading  model 
for  Systems  Test. 

4.  Determine  feasibility  of  an  optimum  scheduling 


model  for  Systems  Test. 
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IV .  Model  Formulation 

Prior  to  the  statement  of  the  model,  the  following 
aspects  of  the  real  system,  developed  during  the  problem 
definition  phase,  are  recalled: 

Each  product  is  assigned  to  a  test  station  and 
testor(s).  It  remains  at  that  location,  with  that  testor  or 
testors,  until  processing  is  complete.  However,  due  to  the 
limited  number  of  second  shift  testors,  processing  on  second 
shift  is  reserved  for  products  which  would  otherwise  be 
completed  late  or  within  a  certain  factor  of  the  due  date  if 
processed  completely  on  first  shift.  Each  test  station  is 
not  compatible  with  every  product  and  each  operator  is  not 
qualified  to  test  every  unit.  There  are  also  various 
"pairwise"  station  constraints  that  may  exist,  where 
assignments  in  one  station  affect  the  product  compatibility 
in  another  station  due  to  limited  resources  shared  between 
the  stations.  Additionally,  operator  productivity  varies 
with  individual  experience  levels.  Another  aspect  which 
must  be  considered  is  that  parallel  systems  must  have  their 
components  tested  separately  prior  to  system  testing.  The 
final  key  consideration  is  the  fact  that  each  product  has  an 
assigned  due  date. 

Therefore,  for  a  complete  scheduling  model,  equations 
must  be  developed  to  ensure  that  a  schedule  meets  the  fol¬ 
lowing  constraints: 

1.  Limited  resources; 

2.  Precedence  relations  between  jobs; 
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3.  Product  due  dates; 

4.  Limited  use  of  second  shift;  and 

5.  Product  assignment. 

During  the  problem  definition  phase  it  was  noted  that 
there  were  difficulties  in  both  scheduling  within  Systems 
Test  and  capacity  planning.  Models  for  both  scheduling  and 
capacity  will,  therefore,  be  presented. 


Scheduling  Model 
Definitions 


i  = 

j  = 

k  = 


product  number,  i  =  1,  ...,  I;  1=  number  of  sys¬ 
tems  to  be  tested  in  planning  horizon. 


test  station  number,  j  =  1, 
test  stations. 


first  shift  testor  number,  k  =  1, 
number  of  first  shift  testors. 


.  .  ,  J ;  J  =  number  of 

,  . ,  K;  K  = 


m  = 


second  shift  testor  number,  m  =  K+l,  ...,  M;  M  = 
total  number  of  testors. 


t 


time  period,  t  =  1,  ...,  T;  T  =  last  time  period 
of  planning  horizon. 


®i(jk)  " 
6i(jm) 

Xi (jk)t 


a .  = 
l 

FSP .  = 

l 


efficiency  index  for  product  i  tested  on  test  sta 

tion  j  with  first  shift  testor  k.  0  e.  1. 

J  l  jk 

efficiency  index  for  product  i  tested  on  test  sta 

tion  j  with  second  shift  testor  m.  0  e. .  1. 

J  13m 

a  variable  which  is  1  in  period  t  if  product  i  is 
scheduled  to  start  testing  on  test  station  j  with 
first  shift  testor  k  at  time  t;  0  otherwise, 
x.  . ..  .  need  not  be  treated  as  a  variable  in  all 
periods,  since  it  equals  0  for  t  a^. 

arrival  period  of  product  i. 

first  shift  processing  time  =  number  of  periods 
required  to  test  product  i  entirely  on  first 
shift . 
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DT\  =  due  period  of  product  i. 

u^,v^  0  and  integer. 

.  =  a  variable  which  is  1  if  a  product  would  be  late 

x  if  processed  entirely  on  first  shift  and  0  other¬ 

wise  . 

Ri(  k)t  =  a  variable  which  is  1  in  period  t  if  product  i  is 
J  being  processed  in  station  j  by  first  shift  testor 

k  at  time  t,  0  otherwise.  Defined  only  for  prod¬ 
ucts  not  requiring  second  shift  testing. 

R .  / ■ k »  =  a  vector  of  dimension  T  of  0's  and  I's,  repre- 
-1  senting  the  first  shift  resource  occupancy  in  the 

horizon  by  product  i  and  pair  (jk)  when  second 
shift  is  not  used. 

Yi('k)t  =  a  variable  which  is  1  in  period  t  if  product  i  is 
-*  being  processed  in  station  j  by  first  shift  testor 

k  at  time  t,  0  otherwise.  Defined  only  for 
products  requiring  second  shift  testing. 

Yi(ik)  =  a  vector  of  dimension  T  of  0's  and  l's,  repre- 

-1  senting  the  first  shift  resource  occupancy  in  the 

horizon  by  product  i  and  pair  (jk)  when  second 
shift  is  used. 

Zi(  'm)t  =  a  vari-akle  which  is  1  in  period  t  if  product  i  is 
being  processed  in  station  j  by  second  shift 
testor  m  at  time  t,  0  otherwise.  Defined  only  for 
products  requiring  second  shift  testing. 

Z.  ...  =  a  vector  of  dimension  T  of  0's  and  l's,  repre- 

J  senting  the  second  shift  occupancy  in  the  horizon 

by  product  i  and  pair  (jm). 

P..  ,.  =  fFSP./2~)  =  first  shift  processing  time  for  prod- 

^  ucts  Requiring  second  shift. 

Pf(jm)2  =  second  shift  processing  time. 

r.  =  amount  of  resource  of  type  c  required  on  product 


RSc  (i  1  '  2) =arnount  resource  of  type  c  available  between  the 
C  -1  -1  test  stations  j  ^  and  j  2 . 


Since  each  product  must  be  assigned  to  a  test  station 
and  a  testor,  this  assignment  will  be  treated  as  a  pair, 
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"station,  operator",  either  on  first  shift  or  second  shift. 
If  an  operator,  or  a  test  station,  is  not  eligible  to  work 
on  a  product  or  class  of  products  the  productivity  is  zero 
for  that  product  and  "station,  operator"  pair.  This  infor¬ 
mation,  along  with  the  operator's  efficiency  rating  is 
included  in  the  efficiency  indexes,  e. and  e. .  .  In  this 
way  varying  processing  times  is  possible  in  that  the  stan¬ 
dard  allotted  test  time  when  divided  by  the  efficiency  index 
yields  the  processing  time  as  a  function  of  the  "station, 
operator"  pair,  i.e.,  FSP^  =  FSP^(jk). 


Constraints 

Product  Assignment 

Testing  on  each  product  must  be  started  sometime  within 
the  planning  horizon,  [1,  Tj ,  in  a  compatible  test  station 
with  an  eligible  testor.  The  interval  for  start  of  test  is 
further  reduced  for  each  product  by  its  arrival  time  and 
Pi(jk)l'  t^ie  Processi-n9  time  requirement  on  first  shift  for 
those  products  tested  on  second  shift.  It  should  also  be 
noted  that  second  shift  testing  on  a  product  will  not  occur 
until  a  first  shift  testor  has  been  assigned.  An  assignment 
constraint  is  required  to  ensure  that  a  product  is  assigned 
to  only  one  "pair"  in  the  horizon. 


T-P 

'  ( J  k )  1  V  _  , 

4-  i(ik]t 

t  =  a . 

j 


(1) 


,  for  each  product  i.  It  should  also  be  noted  that  the 
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"pair"  (jk)  will  not  be  defined  if  either  the  product  is  not 
compatible  with  the  test  station  or  the  operator  is  not 
qualified  to  test  the  product.  In  addition,  the  constraint 
automatically  requires  that  a  product  is  assigned  to  only 
one  pair  (jk)  for  all  t.  The  constraint  is  actually  a 
"multiple  choice"  constraint  over  all  t  and  all  pairs  (jk). 

Additionally,  a  product  assignment  or  capacity  con¬ 
straint  is  required  to  ensure  that  test  stations  and  testors 
are  assigned  at  most  one  job  at  a  time.  Obviously,  there 
are  two  ways  in  which  to  proceed.  First  is  to  bound  the 
"station,  operator"  pairs,  (j,k)  and  (j,m).  The  second  is 
to  bound  the  test  station,  first  shift  testors,  and  second 
shift  testors  individually.  As  stated  previously,  the 
"pair"  model  will  be  presented. 

Let  y^  ( ,  z i  ( jm)  '  an<^  Ri(jk)  vectors  of  dimension 
T  representing  the  resource  occupancy  in  the  horizon  by¬ 


product  i  and  "pair"  (jk)  and  (jm) .  Then, 

?  (yi (jk)  +  Ri (jk) 1  -  lT  (2) 

,  for  each  "pair"  (jk)  (1st  Shift)  and 

^  Zi(jm)  -  1 T  (3) 

m  for  each  "pair"  (jl)  (2nd  Shift).  1  represents  a  vector 

of  ones  of  dimension  T.  It  should  be  noted  that  constraints 


(9),  (10),  (11),  (12),  (13),  and  (14)  will  further  restrict 

the  variables  R.  ,  ,  Y.  ,  and  Z  .  .  .  .  so  that  if  R.  ,  . 

l  (]k)  i (3k)  l (3m)  1 (3k) 

exists  then  Y...,.  and  Z.,.  .  will  not  and  if  Y.,  . .  exists 
1 (3k)  1  (3m)  1  (3^) 

then  Z ^  ( jm)  will  exist  but  not  ( j )  * 
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Limited  Resources 

In  some  instances  resources,  such  as  batteries,  power, 
or  voltage  set  ups,  may  be  shared  between  stations.  A 
"pairwise"  station  constraint  then  exists  when  the  avail¬ 
ability  of  these  resources  is  reduced  or  denied  to  one 
station  when  in  use  in  another.  This  consequently  reduces 
the  compatibility  of  the  station.  These  types  of  con¬ 
straints  take  two  forms. 

First,  the  use  of  a  specific  resource  in  one  station 
denies  its  use  in  another.  For  example,  only  one  50  Hz  unit 
may  be  tested  in  stations  j^,  and  j  ^  at  any  time. 


E  E  (V.,.  +  Y.,.  ,,  +  R.,.  ,,  +  R...  ,,  <  1.,.  (4) 

10,10  i(J2k)  iCjjk)  ifj-,10  -  1 


This  is  to  be  written  for  every  pair  and  j  0  thus  speci¬ 
fied  where  i  ranges  over  those  products  requiring  the  spe¬ 
cific  resource. 

Secondly,  the  use  of  a  specific  resource  in  one  station 
reduces  its  use  in  another  so  that  the  total  use  may  not 
exceed  a  certain  value.  For  example,  the  maximum  power 
requirement  in  test  stations  and  j ^  is  limited  to 


RS 


C  ( j  1  j  2 ) 


EE 


^  ric  ^  i  ( j  j k ) 


+  R. 


ifjjM 


+  Y. 


U,U 


R . , .  ,,11 

i  0,kl 


<  RS 


(5) 


Again,  this  is  to  be  written  for  every  pair  j  ^  and  j2,  thus 
specified  where  i  ranges  over  those  products  requiring  the 
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specific  resource.  Both  of  the  forms  of  the  resource  con¬ 
straint  may  be  expanded  to  any  subset  of  stations  j  ^ , 
j  2  »  -  -  instead  of  a  single  pair  j  ^  and  j  ^  • 

Inclusion  of  the  terms  Z.  .  and  Z.  ,  .  .  are  not 

i  l  (]2m) 

necessary  in  constraints  (4)  and  (5).  This  is  due  to  the 
fact  that  products  must  be  assigned  on  first  shift  prior  to 
second  and  that  the  test  station  assignment  remains  constant 
throughout  the  processing  period.  Therefore,  any  conflicts 
would  have  been  resolved  prior  to  second  shift  processing. 


Precedence  Relations 

Precedence  relations  exist  for  the  processing  of  paral¬ 
lel  systems.  A  parallel  system  is  composed  of  a  number  of 
components  offering  redundancy  of  protection  to  the  user. 
Therefore,  each  component  of  the  system  must  be  tested  as  a 
single  unit  prior  to  testing  as  a  system.  If  the  number  of 
periods  required  to  process  a  component  n  of  Product  i  is 

designated  by  PP .  ,  then 

n 

T 


and , 


t-?.  '"‘.OMt  *  \0Ht  *  :inU",t) 
1 

n 


T-P. 


i  (  j  k )  I 
n 

t — \ 

t=a. 

i 


T-P. 


tX.  r ♦  PP.  -  1  < 
xn(jk)t 


i  C  j  K )  1 

£  tx 

t  =  a. 


i (jklt 


(6) 


where  n  =  1,  2,...N;  N  =  number  of  components  in  parallel 


product  i.  This  would  be  written  for  each  component  of  each 
parallel  system  scheduled  to  be  tested  in  the  horizon. 
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Second  Shift 

Second  shift  operation  takes  place  only  if  time  is  net 
sufficient  to  complete  the  processing  of  the  product  by  its 
due  time  or  within  a  specified  period  of  its  due  time  on 
first  shift  only.  A  product  that  can  be  finished  by  the 
specified  period  of  its  due  time  is  not  permitted  to  have 
second  shift  testing  and  thus  finish  it  even  earlier. 


T-P 


z 

t  =  a . 


tX.rn  +  2  ( FSP .  -  1)  -  (DT .  +  C)  =  u.  -  v. 

i(ik)t  l  l  l  i 


(7) 


,  for  all  i,  where  u^,  v^  _>  0  and  integer  and  C  is  the 
specified  period  of  the  due  time  before  which  products  must 
be  completed  without  requiring  second  shift. 


u  .  <  6. K 

x  -  1 

vi  1  (l“«i)K 


(8) 


,  where  K  is  a  large  number  greater  than  zero  and  each 
constraint  is  written  for  all  i.  This  constraint  with  the 
proper  valuation  of  u^  and  v^  in  the  criterion  function  will 
result  in: 

1.  If  product  i  is  late,  completed  after  the 

specified  period  before  the  due  time,  when  processed 

only  on  first  shift,  u.  >  0,v.  =  0,  and  6.  =  1. 

i  i  l 

2.  If  product  i  is  early,  completed  prior  to  the 
specified  period  before  the  due  time,  when  processed 
onlv  on  first  shift,  u.  =  0,  v.>0,  and  f~.  =  0. 

'  l  l  l 

From  the  assignment  constraint,  x.,.,.^  -  1,  for  some 
^  l (gk) t 

t  =  t  ,  and  x.  .  , .  =  0,  for  all  t  ^  t  .  Time  must  now  be 

O  1  |  J  K )  t  O 
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blocked  to  complete  processing  on  first  shift  only  if 
product  i  is  "early",  or  on  first  and  second  shift  if  it  is 
"late"  . 


(i)  Early:  Let 
"regular  time"  (first 


R  •  ,  . ,  . 
i  (nk) 

shift 


t  =  0,  1  variable,  representing 
testing  only)  operation. 


(jk) 


2 

t  =  a. 


R.  ,  .  < 

l  ( j  k )  t  — 


(1 


<5.)FSP.  + 

1  1 


,  for  each  i,  and 


(9) 


t+Jd-SP. -1) 


V-V 

(Jk) 


2 

T=t 


R. 


i  (jk)x 


FSP .  X.  ...  . 
l  l  ( i  k )  t 


6.K 

l 


(10) 


,  for  each  i,  t. 

(ii)  Late:  Here  there  are  two  modes  of  operation  to 
consider.  The  first  is  to  split  the  job  evenly  between 
first  shift  and  second  shift.  This  results  in  an  Integer 
Linear  Programming  model.  The  other  is  to  do  as  much  as 
possible,  without  exceeding  the  product  due  time,  on  first 
shift  and  assign  the  remainder  to  second  shift.  This 
results  in  a  non-linear  model.  The  ILP  model  will  be  pre¬ 
sented  . 

Let  yi (  t  =  0,  1  variable  representing  first  shift 
operation , 

hum  -  rFSP(i(jk|/2l  and  bum  -  FSbuM  - 

Pi  ( jk)  1  '  where  Pi  ( jk)  1  and  Pi  ( jm)  2  differ  b>’  at  most  one- 


,  for  all  i  and  (jm)  pairs.  It  should  be  noted  that  if 


5 . =0 ,  constraint  (11)  and  (13)  force  all  Y .  .  . ,  .  and  Z  .  . 

1  1  \  ]  K )  t  llJiTiJi, 

equal  to  zero.  Then,  the  remaining  sets  of  constraints, 

(12)  and  (14),  are  automatically  satisfied. 

However,  if  5^=1  (late  job),  constraints  (11)  and  (13) 
are  automatically  satisfied,  but  the  remaining  pairs  are 


constrained  only  at  t  ,  for  which 
2  o 


t  +2  (P.  ......  - 1 ) 

o  i (ik)  1 

y1  Y .  , . ,  .  1\  ...  .  .  and  V'  Y .  , . ,  ,  <P.,.,,, 

i(jk)T-  i  (jkj  1  ^  i(jk)T-  j  (j  k  1 1 


to+2(Pi(j»0l"1) 


T=t 


T=t 


forcing 


t  +2(P.  ,.-l. 

o  i  (jk)  1  ) 

V  =  p 

T^t  i(jk)T  i  C i k ) 1 

o 


implying  that  Y  t  ^  =  i  for  every  first  shift  period 

between  t  <  T  <  (t  +2n';n;  ,-11). 

o—  —  v  O  1(1  111 


Similarly, 

t  +2P.  , .  .,-1 
o  i(jm)2 


V  Z .  - .  >  P  n  and 

Zj.  i(jm)T-  i(jm)2 


t  +2P.,.  ,,-1 
o  i  ( j  m)  2 


2 

T=t 


<  P 

'i  (jm)  t  —  i  (jm)2 


,  .  t  +2P.  ,.-1 

forcing  o  i(jnij- 

^  Z  =  P 

£->  *"i  ( jm)  t  i(jm)2 

o 


implying  that  Z.  1  for  every  second  shift  period  in 

t  <  t  <  (t  +2P .  . .  _-l). 
o  —  o  l  ( jm)2 


Objectives 

Determining  when  a  product  should  be  tested  depends  on 


the  desired  objective.  As  determined  in  the  problem 
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definition  phase  the  immediate  area  of  concern  is  lateness. 
Therefore,  formulations  for  the  following  objectives  are 
presented : 

1.  Minimize  number  of  late  jobs. 

2.  Minimize  total  lateness. 

3.  Minimize  total  cost  of  lateness. 


Criterion  1  (Minimize  number  of  late  jobs) : 

A  product  is  late  if  it  is  completed  after  its  due 
date.  Therefore,  if  at  t=DT^  +  1,  the  period  after  the  due 
date,  Ri(jk)t  +  ^  i  ( j  kl  t  +  “i ( j  k )  t  =  l’ 

the  product  is  late  because  it  is  still  occupying  the 
resource.  Let  0,1  variable,  then, 


'T'  n1' 

(jml  (jk)  t  =  ivr  +1 


^Ri(jk)t  +  ^i(jk)t  "i'.iinit  — 


for  all  i,  where  K  is  a  large  number  >  0. 

If  the  left  side  of  this  inequality  is  greater 

some  processing  must  be  taking  place  after  the 

a.  must  =  1:  otherwise  a.  =  0.  Therefore,  to 
l  i 

number  of  late  jobs 

Minimize  X-  a- 
i 

is  the  desired  function. 


n.K 

3 


than  zero 
due  time  and 
minimize  the 


Criterion  2  (Minimize  total  lateness) : 

Again,  a  product  is  late  if  it  is  completed  after  its 

due  period,  DT.  .  Therefore  a  product  is  late  if  ( R .  .  ,  .  + 

i  l  (  3  K }  t 

Y .  ...  .  .  +  Z...  . . )  >  0  in  the  periods  after  the  due  period, 
i(]k)t  i ( jm) t  r  t 
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indicating  that  testing  is  not  complete.  If  total  lateness 
is  to  be  minimized,  then 


Minimize  V 


?  (i)  (i)  t.DT.,1  CBi0k5t  *  Vifjk)t  * 


is  the  desired  function,  where  lateness  is  defined  as  the 
number  processing  periods  required  after  the  due  period. 


Criterion  3  (Minimize  total  cost  of  lateness) : 


If  a  penalty  or  cost  of  is  assessed  when  product  i 

is  not  completed  by  period  t,  the  total  cost  of  lateness  is 
minimized  by  the  function 

Minimize  ^  ^  „  T 

(jk)  (ji  t=DT.+l"itCRi^iklt  +  + 


This  expression  reduces  to  total  lateness,  Criterion  2,  if 
all  Wit  =  1. 

To  ensure  that  constraints  (7)  and  (8)  function  as 
intended,  the  variables  u^  and  must  also  be  a  part  of  the 
chosen  objective  function.  This  is  accomplished  by  adding 
the  following  term  to  the  chosen  Criterion: 


+  2  M  f  u .  +  v . 


M  ( u .  +  v  . ) 
i  l 


where  M  is  a  positive  number  sufficiently  large  to  ensure 

that  the  contribution  of  the  additional  term  is  less  than 

that  of  any  (R .  ...  .  4  y .  ,  ■ ,  .  .  +  Z .  .  .  .  ). 

7  i  (ilk)  t  l  (3k)  t  1  (3m) t 


Computational  Complexity 

The  computationa 1  complexity  of  this  formulation  is 
easily  demonstrated  by  a  realistic  example.  Using  the  data 
in  Table  4.1  the  necessary  number  of  constraints  for  each 
constraint  as  presented  in  the  formulation  is  tabulated  in 
Table  4.2. 


Table  4.1  Example  Planning  Data 


InPut  1  Input  2 


Horizon 

40 

shifts  20 

Number  of 

products  to  be  tested 

55 

22 

Number  of 

test  stations 

14 

14 

Number  of 

first  shift  testors 

13 

13 

Number  of 

second  shift  testors 

7 

7 

"Pairwise" 

stations 

3 

3 

Number  of 
Number  of 

parallel  systems  to  be  tested 
components  in  each  parallel 

4 

O 

t- 

system  2  2 


Table  4.2  Number  of  Scheduling  Constraints 


(Constraint  Number 


Data  1 


Number  Required 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 

(11) 

(12) 

(13) 

(14) 


55 

182 

98 

3 

3 

8 

55 

110 

55 

2,200 

10,010 

14,560 

5,390 

718,830 


Data  2 


2 

18 

9 


4 

22 
44 
22 
440 
4 ,004 
7 ,280 
2,156 
358,876 


TOTAL 


751,559 


443,156 


It  is  readily  apparent  by  the  total  number  of  constraints 
required  that,  even  though  the  scheduling  model  as  presented 


ro  ro  ro 


in  the  tormulation  is  an  integer  linear  programming  model, 
it  is  unsolvable  and,  therefore,  another  approach  must  be 
used.  However,  the  first  objective  of  this  study:  Develop 
a  model  to  describe  the  System  Test  assignment  functions 
with  regard  to  resource  constraints  has  been  met. 


Scheduling  Alternative 

Returning  to  the  scheduling  model  and  Table  4.2  it  is 
obvious  that  the  number  of  required  constraints  balloon 
after  constraint  (10).  The  remaining  constraints  all  deal 
with  second  shift  operation.  If  these  constraints  were 
eliminated  and  the  assumption  made  that  all  processing 
occurs  on  first  shift  the  scheduling  constraints  would 
reduce  to  the  following: 

Product  Assignment 

(1)  T-^Pi 


,  for  each  product  i. 


(2)  \  Ri  (jk)  -  1T 

,  for  each  pair  ( jk) . 


Limited  Resources 


(4)  £  ^  CRi  Cj  jk;  +  Ri(i2k)}-  !T 


,  for  each  pair  j^  and  j  "pairwise"  constrained  and 
where  i  ranges  over  the  products  requiring  the  specific 


resource . 
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(5) 


<0 


V 

k 


r .  (R. . . 

1CK  i(] 


k) 


R.  ,  ■  ,  J  <  RS 
i(j-,k)  - 


c  C  j  j  j 


2^T 


,  for  each  pair  and  "pairwise"  constrained  and 
where  i  ranges  over  the  products  requiring  the  specific 
resource . 


Precedence  Relation 


(6)  T^Pi 


T-FSP. 


?  «i_oot *  £  Rifjk)t'1-  ^  txi(jkn 


t=a.  n  t=a .  n  t  =  a. 

ii  i 


n  n 

,  for  each  component  of  each  parallel  system. 


Z 


(9)  L  L  R. 


* 


(jk)  t=a . 


i(jk)t 


=  FSP . 


,  for  all  i 


t+2(FSP.-l) 


z 

(jk) 


T=t 


R...,,  >  FSP. 

i  (jk)t  -  i  i  (jk)t 


,  for  each  i ,  t 


Again  using  the  data  of  Table  4.1  the  computational  com¬ 
plexity  of  the  reduced  formulation  is  represented  in  Table 
4.3  by  the  number  of  required  constraints. 
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Table  4.3  Number  of  Reduced  Scheduling  Constraints 


Number  Required 


Constraint  Number  Data  1  Data  2 


(1)  55  22 

(2)  182  18 

(4)  3 

(5)  3 

(6)  8 

(9)  55  2 


TOTAL  306  236 


The  totals  are  now  well  within  the  bounds  of  a  mixed  integer 
programming  package. 

Continuing  the  reduction  to  the  objectives  of  the 
scheduling  model,  the  reduced  formulation  yields: 


Criterion  1  (Minimize  number  of  late  jobs) : 

_  T 

Z  Z  RifiiOt-aiK 

(jk)  t=DT  +1  J 

,  for  all  i,  where  oo  =  0 ,  1  variable  and  K  is  a  large 
positive  number.  The  desired  criterion  function  remains 
Minimize  a. 


Criterion  2  (Minimize  total  lateness) : 


Minimize  Z  Z 


T 

i  (jk)  t=OT.+l 


A 


(N  m  fO  'T  CNJ 
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Criterion  3  (Minimize  Total 


Minimize  ^  2  X) 

i  (jk)  t=DT.+i 


cost  of  lateness): 


"it  Ri(jk)t 


An  iterative  approach  to  the  complete  scheduling  prob¬ 
lem  is  now  proposed  and  illustrated  by  a  simple  example  with 
only  one  test  station.  The  following  table  of  data  will 
describe  the  products  of  the  example. 


Step  1:  Solve  the  scheduling  problem  with  the  reduced  formu¬ 
lation  . 


i 

1 

2 

3 

mm 

1 

5 

mm 

7 

11 

m 

figfpP: 

3 

5 

m 

I 


Step  2:  Identify  those  products  which  are  late  and  whose 
lateness  is  less  than  the  completed  processing  time  prior  to 
the  period  DT^  +  1.  Then  reduce  the  processing  time  on 
these  products  by  their  lateness.  In  the  example,  the  total 
processing  time,  FSP^,  of  product  2  is  5  periods  but  it  was 
completed  2  periods  late.  Therefore,  the  new  processing 
requirement,  FSP^ ' ,  for  product  2  would  be  3  periods. 


1  3  5  7  9  11  13  15  17  19  21  23  25 

Step  3:  Solve  the  scheduling  problem  again  with  the  reduced 
formulation  and  the  new  processing  requirements,  FSPV,  for 
the  specified  products  of  step  2.  This  second  solution  will 
then  serve  to  compact  the  first  solution  after  the  reduction 
of  processing  times. 


Step  4:  Repeat  Step  2  for  any  products  in  the  category 
which  have  not  previously  been  reduced  or  whose  previous 
reduction  plus  current  lateness  is  less  than  current 
processing  time  prior  to  the  period  DT^  +1.  In  the 
example,  product  3's  total  processing  time,  FSP^,  is  4 
periods  but  it  was  completed  1  period  late.  Therefore,  the 
new  processing  requirement,  FSP^ ’ ,  for  product  3  would  be  3 
periods.  As  another  example,  if  a  product,  whose  total 
processing  requirement  is  8  periods,  had  been  reduced  2 
periods  in  Step  2  and  is  now  after  Step  4  one  period  late, 
the  product's  processing  time  would  again  be  reduced  by  one 
period . 


Step  5:  This  step,  in  essence,  schedules  the  second  shift. 
Using  the  reduced  formulation,  solve  a  reduced  problem  con¬ 
sidering  only  the  products  whose  processing  requirements 
have  been  curtailed  in  Step  2  or  Step  4.  Since  we  are  now 
concerned  with  second  shift,  the  test  station  for  each  of 
these  products  is  predetermined  by  the  first  shift  solution 
of  Step  3.  Therefore,  all  other  combinations  (ij)  will  be 
undefined.  The  arrival  times  of  the  products  are  set  by  the 
(first  shift  processing  start  time)  +  1 ,  also  determined  in 
the  solution  of  Step  3,  to  force  arrival  at  the  next  second 
shift  period.  The  due  periods,  are  also  set  by  the  solution 
of  Step  3,  in  that  the  minimum  of  the  (processing  stop  time) 

-  1  and  the  (due  period)  -  1  is  chosen.  This  tries  to  force 
second  shift  testing  to  occur  between  the  periods  that  the 
product  is  undergoing  first  shift  processing.  The  processing 
time  is  the  amount  of  the  reductions.  The  following  table 
displays  the  reduced  problem  for  the  example: 


i 

2  3 

a . 

DT1 

P1 

l 

8  14 

10  16 

2  1 

The  testors  to  be  assigned  are  now  the  second  shift  testors, 
i.e.,  k  is  replaced  by  m  in  the  formulation. 
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When  this  solution  is  combined  with  the  solution  of 
Step  3  and  the  reductions  of  Step  4,  a  complete  schedule  is 
obtained . 

2nd 


Shift 


1st 


Due  to  the  scarcity  of  second  shift  testors,  a  quali¬ 
fied  second  shift  testor  may  not  be  available  to  test  a 
product  between  the  periods  that  the  product  is  undergoing 
first  shift  testing.  With  the  preceding  iterative  approach, 
this  would  result  in  the  second  shift  testing  of  the  product 
extending  past  or  occurring  after  the  (first  shift 
processing  stop  time)  +  1.  This  would  be  a  violation  of  the 
established  second  shift  due  date  and  Systems  Test  operating 
procedure.  The  following  example,  with  two  test  stations, 
will  illustrate  this  situation  and  the  rectifying  procedure. 


. 

1 

1 

2 

3 

4 

5 

a  . 

1 

5 

7 

1  3 

13 

DT1 

7 

9 

1 1 

17 

17 

FSP1 

l 

3 

5 

4 

5 

5 
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Step  5:  (Assuming  only  one  qualified  second  shift 

testor  available) 


rm — rm — rrrr 


ttti — rm — rr-n — rrrr 

4  16 


t  Ep  t  l  t  r  n  i  r  i  t  rr 

20  22  24  26 


2  4  hb  ts  10  12  ^14 

(\  o  \  A  ,  \. 

a->  a,  DT  DT,  a.  r 

2  3  2  j  4,5 


18 


DT 


4,5 


Combining  the  solution  of  Step  5  with  the  solution  of  Step 
4,  it  is  apparent  that  the  second  shift  testing  of  product  4 
extends  past  its  first  shift  processing  stop  time  in  excess 
of  one  period. 


It  is  now  necessary  to  return  to  the  solution  of  Step  4  to 
determine  if  this  "lateness"  can  be  processed  on  first  shift 
in  a  "processing  gap"  produced  by  the  previous  reduction  of 
processing  time  in  this  step.  If  so,  a  feasible  solution  is 


obtained . 


1 


— rm rm rm — rm rm rm rm — m~i rm rm rm rm rm - 

1  5  5  7  9  11  13  15  17  19  21  23  2S 

*  V  V 

DT  DT'  DT,  DT,  r 

X  i—  *3  *4  j 

In  this  example,  it  was  possible  to  find  such  a  gap. 
However,  if  a  product  had  been  reduced  both  in  Step  4  and 
Step  2  the  "lateness"  may  not  fit  in  the  "processing  gap". 

In  this  situation,  it  would  be  necessary  to  return  to  the 
solution  of  Step  2  followed  by  Step  5  without  the  condensing 
step,  Step  3,  and  the  second  reduction  step,  Step  4.  This 
would  ensure  that  any  subsequent  "lateness"  would  fit  in  the 
"processing  gap." 

The  preceding  iterative  approach  to  the  complete 
scheduling  problem  terminates  with  a  feasible  solution  in  a 
maximum  of  seven  steps.  However,  the  solution  is  not  opti¬ 
mal  due  to  the  possible  requirement  of  a  return  to  Step  2 
and  the  subsequent  loss  of  the  condensing  step.  If,  how¬ 
ever,  the  number  and  qualifications  of  the  second  shift 
testors  mirrored  the  first  shift  testors,  an  optimal 
solution  could  be  obtained  through  this  procedure  by 
returning  to  Step  3  after  Step  4  until  no  reduction  could 
take  place  in  Step  4.  Since  reductions  are  not  allowed  to 
exceed  the  completed  processing  time  prior  to  the  period  DT. 


+  1,  it  is  ensured  that  the  second  shift  processing 
requirement  will  not  exceed  the  time  on  first  shift. 
Additionally,  since  testors  are  mirrored  on  first  and  second 
shifts,  it  is  also  ensured  that  a  qualified  second  sr.ift 
testor  will  be  available  during  the  second  shift  periods 
following  any  first  shift  testing. 

Capacity  Model 

If  instead  of  scheduling,  the  focus  is  placed  on 
capacity,  the  model  may  be  greatly  simplified  by  assumptions 
with  the  objective  of  obtaining  a  gross  capacity  estimate. 
The  assumptions  made  in  this  formulation  are  (1)  that 
testors  are  readily  available  and  (2)  that  products  are 
immediately  available. 

In  contrast  to  the  scheduling  model,  neither  testors 
nor  time  are  a  parameter  of  the  formulation.  Therefore, 
testor  qualifications  and  specific  product  arrival  times  are 
not  taken  into  consideration.  By  not  considering  testor 
qualifications,  all  testors  are  considered  capable  of 
testing  all  products  and  only  the  number  of  available 
testors  has  an  impact  on  the  solution.  By  not  considering 
arrival  times,  the  dynamic  structure  of  the  complete 
scheduling  model  is  converted  to  a  static  structure  by 
ignoring  the  nonsimultaneous  arrivals  of  products  to  Systems 
Test.  Additionally,  unlike  the  complete  scheduling  model,  a 
criterion  does  not  exist  for  second  shift  testing  and, 
consequently,  any  product  may  be  tested  on  second  shift. 


In  order  to  avoid  a  third  assumption  that  all  test 
stations  are  the  same,  which  would  be  questionable,  the 


assignment  of  a  product  to  a  test  station  will  be  treated  a 
a  "pair"  (ij) ,  as  in  the  scheduling  model  with  (jk)  and 
(jm) .  In  this  way  a  distinction  may  be  made  between 
stations.  This  is  accomplished  in  the  following  manner;  if 
product  i  is  not  compatible  with  station  j  then  the  "pair" 
(ij)  is  undefined  and  the  variables  x,. ..  and  y,.  .  are 

subsequently  reduced.  This  will  ensure  that  a  product  is 
assigned  only  to  a  compatible  test  station. 

Definitions 

W1  =  number  of  first  shift  testors. 

W2  =  number  of  second  shift  testors. 

i  =  product  number,  i  =  1,  2,  I;  1=  number  of 

systems  to  be  tested  in  planning  horizon. 

j  =  test  station  number,  j  =  1,  2,  J  =  number 

of  test  stations. 

=  processing  time  of  product  i. 

X, . ■  ,  =  duration  of  job  i  in  test  station  j  on  first 

'Vi  Shift. 

Y, . ..  =  duration  of  job  i  in  test  station  j  on  second 

1  shift. 

H j  =  number  of  periods  available  for  testing  on  first 
shift  =  planning  horizon  x  min  [J,  Wl], 

=  number  of  periods  available  for  testing  on  second 
shift  =  planning  horizon  x  min  [J,  W2], 

DT^=  due  period  for  product  i. 

6  .. ..  =  a  variable  which  equals  1  if  product  i  is 
Jrested  in  station  j  and  equals  0  otherwise. 
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Constraints 
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Product  Assignment 

Testing  for  each  product  must  be  started  sometime 
within  the  planning  horizon  in  a  compatible  test  station. 
An  assignment  constraint  is  required  to  ensure  that  each 
product  is  tested  in  only  one  test  station. 


J 

for  all  i,  where  6.  is  a  0,  1  variable 

l 


for  all 


X.. ..  <  6. . K 
(ij)  “  ij 

(ij)  pairs,  where  K  is  a  large  positive 


number . 


,  for  all  (ij)  pairs,  where  K  is  a  large  positive  number. 
This  set  of  constraints  ensures  that  x....  and  v....,  are 
equal  to  zero  if  product  i  is  not  assigned  to  test  station  j 
and  greater  than  zero  if  product  i  is  assigned  to  test 
station  j.  It  should  be  noted  that  the  "pair"  (i])  will  not 
be  defined  if  the  product  is  not  compatible  with  the  test 
station . 


Processing  Time 

In  this  formulation  any  product  may  be  tested  on  second 
shift.  The  constraint,  however,  is  that  the  duration  of 
processing  time  on  first  and  second  shift  must  equal  the 
processing  time  of  the  product.  Since  testors  are  not  a 
parameter  of  the  formulation,  consideration  of  tester  effi¬ 
ciency  in  calculating  processing  time,  as  in  the  complete 
scheduling  model,  is  not  possible.  Therefore,  variable 


processing  times  are  not  permitted  and  the  processing 
requirement  for  each  product  is  assumed  known. 
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x(ij)  +  Y(ii)  =  Pi'  for  each  ”Pair"  (i3>  (2) 

Processing  time  on  second  shift  must  also  be  less  than  the 

processing  time  on  first  shift  for  each  product.  This  is 
due  to  the  fact  that  second  shift  testors  are  more  scarce 
than  first  shift  testors  and  it  is  desirable  to  limit  their 
use.  Additionally,  it  is  a  result  of  the  constraint  that  a 
product  may  not  be  tested  on  second  shift  without  first 
being  tested  on  the  immediately  preceding  first  shift. 


Y(ij)lX(ij)'  for  each  "Pair"  (ij) 


(3) 


Capacity 

A  constraint  must  also  be  written  to  force  testing  to 
occur  within  the  planning  horizon.  This  is  achieved  as 
follows , 


for  every  j,  where  =  H 


1/J' 


and 


-H:j 


(4) 


,  for  every  j,  where  H_  .  =  H_,_. 

7  J  2j  2/ J 

This  constraint  pair  is  then  a  capacity  constraint  in  that 
the  sum  of  the  testing  durations  are  not  allowed  to  exceed  a 
given  value,  the  capacity  horizon. 
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Objective 

The  objective  continues  to  be  the  minimization  of  the 
number  of  late  jobs.  Since  time  is  not  a  parameter  of  the 
formulation,  however,  a  function  of  time  must  be  introduced 
in  the  objective.  The  following  minimum  slack  function  will 
be  used:  [ (DT . -a . ) -P . ] ,  where  DT .  is  the  "work  day"  that 

the  product  is  due,  a ^  is  the  "work  day"  that  the  product 
arrives  at  Systems  Test,  and  is  the  product's  processing 
requirement  in  shifts.  Work  days  are  numbered  sequentially 
through  the  planning  horizon.  For  example,  if  the  planning 
horizon  covered  four  weeks  and  the  plan  was  to  work  five 
days  a  week,  there  would  be  20  work  days,  numbered  one  to 
twenty.  In  this  way  the  minimum  slack  function  indicates  if 
a  product  will  be  late  when  tested  only  on  first  shift. 

This  indication  is  given  by  a  negative  value  for  the 
function.  For  a  number  of  examples  refer  to  Table  4.4. 

Table  4.4  Minimum  Slack  Function  Example 


H 

DT. 

i 

Comments 

■ 

5 

9 

3 

■ 

B 

5 

4 

3 

Bl: 

Product  arrived  after 
due  date 

B 

5 

9 

5 

mm 

Processing  time 
exceeds  slack 

By  combining  the  minimum  slack  function  with  the  second 
shift  testing  duration,  the  desired  time  link  is  achieved. 
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If  this  function  is  then  minimized, 

Minimize  ££  [(DT.-a.)  -  P.]Yr ... 

i  j  i  i  i  (ij) 

it  forces  the  jobs  with  the  minimum  slack  to  be  tested  on 
second  shift.  It  should  be  recalled  that  constraint  (3) 
requires  the  second  shift  processing  duration  to  be  less 
than  or  equal  to  the  first  shift  duration  for  each  product. 
Therefore,  it  would  be  expected  that  products  with  the  least 
slack  would  receive  equal  processing  on  first  and  second 
shift,  thereby  minimizing  the  chance  of  the  product  being 
late. 

Computational  Complexity 

Using  the  same  technique  to  demonstrate  complexity  as 
with  the  Scheduling  Model  and  with  the  input  data  of  Table 
4.1  the  number  of  required  constraints  are  presented  in 
Table  4.5. 


Table  4.5  Number  of  Capacity  Constraints 


Constraint  Number 

Number 

Data  1 

Required 

Data  2 

(1) 

1,595 

638 

(2) 

770 

308 

(3) 

770 

308 

(4) 

28 

28 

TOTAL 

3,163 

]  ,282 
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Again,  it  is  obvious  by  the  number  of  required  constraints 
that  the  formulation  is  not  within  the  bounds  of  a  mixed 
integer  programming  package. 

The  assumptions  made  in  the  capacity  model  were  very 

broad  and,  therefore,  the  solution,  as  stated  previously,  is 

at  best  a  gross  estimate  of  capacity.  However,  the  solution 

does  give  an  indication  of  feasibility  of  the  desired  load, 

i.e.,  whether  it  violates  the  capacity  constraint  or  not.  A 

preliminary  check  for  this  violation  would  be  to  sum  the 

processing  times  of  the  input  load  and  compare  the  value 

with  the  sum  of  the  capacity  horizons.  If  P.  >  (II  +  tl7) 

i  1 

then  you  may  immediately  assume  that  the  desired  load  is 
infeasible.  However,  if  52  P.  <  (H,  +  H  ')  ,  there  is 

l—l  2 

) 

a  possibility  that  the  desired  load  is  feasible. 

If  it  is  determined  that  the  desired  load  is  infea¬ 
sible,  either  by  the  preliminary  check  or  the  solution  of 
the  capacity  problem,  a  decision  must  be  made  as  to  which 
product  to  eliminate  from  the  desired  load.  This  heuristic 
could  be  based  on  a  number  of  factors,  such  as,  sold  versus 
stock  orders,  earliest  due  date,  shortest  processing  time, 
or  minimum  slack.  This  problem,  therefore,  demands  a  sepa¬ 
rate  treatment  and  will  not  be  investigated  here.  However, 
combining  an  infeasible  solution  with  an  appropriate 
heuristic  to  determine  which  products  to  eliminate  from  the 
desired  load  would  ultimately  yield  a  feasible  solution. 

A  feasible  solution  would  give  an  assignment  of 
products  to  est  stations.  However,  since  the  "pairwise" 
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test  station  constraints  and  testor  qualifications,  as  delin¬ 
eated  in  the  scheduling  model,  were  not  a  part  of  the  capa¬ 
city  formulation,  the  assignment  would  probably  be  infea¬ 
sible.  Additionally,  a  sequencing  link  does  not  exist 
between  first  and  second  shift.  Therefore,  a  product  may  be 
tested  on  second  shift  even  though  the  product  was  not 
tested  on  the  immediately  preceding  first  shift. 

As  stated  previously,  since  time  is  not  a  parameter  of 
the  formulation,  arrival  times  are  not  taken  into  considera¬ 
tion.  By  reducing  the  problem  to  subhorizons,  not  only 
could  more  consideration  be  given  to  arrival  times  but  also 
a  reduction  in  computational  complexity  would  occur  through 
the  consideration  of  less  products  at  one  time.  In  each 
subhorizon  the  product  input  would  be  those  products 
scheduled  to  arrive  in  that  subhorizon  plus  any  products 
eliminated  from  the  previous  subhorizon. 

The  combination  of  the  solutions  from  the  subhorizon 
problems  would  again  yield  an  estimate  of  capacity  for  the 
particular  product  mix  but  with  consideration  of  scheduled 
arrival  times.  Combining  this  measure  of  capacity  with  the 
previously  discussed  iterative  approach  to  scheduling  would 
result  in  a  testing  schedule  approximating  the  capacity  of 
the  system  for  the  specific  product  mix. 

V.  Simulation 


During  the  formulation  of  the  scheduling  model  it 
became  readily  apparent  that  the  model's  complexity  would 
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overshadow  its  usefulness.  Therefore,  in  order  to  meet  the 
second  of  the  project  objectives,  a  simulation  approach  was 
chosen.  Two  advantages  of  this  approach  were  that  it  per¬ 
mitted  controlled  experimentation,  especially  in  regard  to 
resource  constraints,  and  it  permitted  sensitivity  analysis 
by  allowing  manipulation  of  input  variables.  The  objectives 
of  the  simulation  were: 

1.  Develop  a  simulation  to  model  the  Systems  Test 
assignment  functions  with  regard  to  resource 
constraints . 

2.  Compare  and  evaluate  alternative  system  designs  in 
an  effort  to  make  scheduling  procedure  and  operation 
recommendations ; 

3.  Develop  a  sequential  simulation  program,  without 
randomization,  which  would  allow  the  manager  of  an  UPS 
test  facility  to  determine  whether  he  would  be  able  to 
meet  monthly  test  requirements. 

The  simulation  portion  of  the  study  was  accomplished  in  the 
following  four  phases: 

Phase  1.  Formulation  of  a  sequential  simulation  model 
and  program  and  its  subsequent  validation. 

Phase  2.  Collection  of  supporting  data  and  randomi¬ 
zation  of  the  sequential  simulation  model  and  program. 
Phase  3.  Design  of  experiment. 

Phase  4.  Analysis  of  simulation  output. 
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Phase  I.  Sequential  Simulation  Development 

Simulation  Model 

In  formulating  the  simulation  model  a  flow-approach  was 
used.  Using  the  observed  assignment  sequence,  Figure  2.3, 
the  information  from  the  problem  definition  phase,  and  the 
insight  gathered  from  the  scheduling  model  formulation,  a 
sequential  simulation  model  was  developed  which  incorporated 
time  and  organized  queues.  A  flow  diagram  of  the  model  is 
presented  in  Figure  5.1. 


i 


Lgure  5.1  Simulation  Flow  Chart 


(continued) 
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Specification  of  Sequential  Simulation  Components 

The  model  requires  three  inputs:  (1)  test  station 

versus  product  efficiency,  (2)  testor  qualifications  versus 
product  lines,  and  (3)  products  scheduled  to  be  tested.  The 
first  of  these  required  an  update  to  the  Test  Facility 
Capability  Study,  Appendix  B.  In  that  study  it  was  speci¬ 
fied  which  test  stations  could  handle  which  products.  It 
remained  to  determine  the  efficiency  with  which  a  giver, 
product  could  be  tested  at  a  given  test  station.  In  order 
to  arrive  at  this  information  expeditiously  the  supervisor 
and  a  number  of  senior  testors  were  requested  to  rate  their 
test  station  assignment  preferences  for  each  product.  These 
preferences  were  then  treated  as  efficiency  ratings  and  are 
presented  in  Appendix  J.  Testor  qualifications  were 
obtained  from  Systems  Test  personnel  records  and  are 
represented  by  a  zero-one  variable  with  one  indicating 
qualification  to  test  that  specific  product  line  as  shown  in 
Appendix  K.  The  final  required  input  is  the  list  of 
products  scheduled  to  be  tested.  This  is  handled  by  a 
simple  array  of  records  with  each  record  containing  all  per¬ 
tinent  technical  data  for  the  specific  product. 

Four  major  decisions  are  required  in  the  processing  of 
a  product:  (1)  is  a  test  station  available  to  test  the 
product,  (2)  is  a  first  shift  testor  available  to  begin 
testing  the  product,  (3)  is  second  shift  required,  and  (4) 
if  second  shift  is  required,  is  a  second  shift  testor  avail¬ 
able  to  test  the  product.  The  decision  on  whether  second 
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shift  is  required  or  not  is  arrived  at  by  determining  if  the 
number  of  first  shift  periods  from  the  first  shift  tester 
assignment  time  until  the  due  period  is  greater  than  the 
required  processing  time.  If  so,  then  second  shift  is 
required,  since  processing  only  on  first  shift  will  make  the 
product  late. 

Three  major  assignments  are  also  made  in  the  processing 
of  a  product:  (1)  assignment  of  a  test  station,  (2)  assign¬ 
ment  of  a  first  shift  testor,  and  (3)  assignment  of  a  second 
shift  testor,  if  required.  All  of  these  assignments  commit 
the  specific  resource  for  the  duration  of  the  requirement  in 
order  to  ensure  continuity  of  testing. 

Formulation  of  Program 

A  main  objective  in  formulating  the  program  was  to 
maintain  flexibility  for  user  defined  variables.  Examples 
of  used  defined  variables  are  number  of  hours  in  a  shift, 
number  of  days  in  a  work  week,  and  planning  horizon.  By 
allowing  these  elements  to  remain  flexible  it  was  felt  that 
the  utility  of  the  simulation  would  expand. 

In  following  the  flexibility  theme  it  was  also  a 
decided  advantage  to  employ  dynamic  data  structures  for 
queue  management.  It  then  became  relatively  easy  to  order 
and  modify  the  lists,  again  increasing  utility. 

The  program  was  coded  in  Pascal  to  facilitate  the  even¬ 
tual  use  of  the  simulation  by  a  manager  as  a  desk-side  refer¬ 
ence  through  the  use  of  a  personal  computer.  By  length 


alone,  the  program  is  complex  and  with  the  use  of  dynamic 
data  structures  the  complexity  increases  substantially. 

This  may  be  viewed  as  a  disadvantage,  as  an  uninformed  user 
may  have  some  hesitation  to  use  the  product. 

Additionally,  it  should  be  noted  that  the  "pair-wise" 
bay  constraints  as  noted  in  the  scheduling  model  formulation 
and  Appendix  B  were  coded  as  a  procedure  which  was  called 
prior  to  any  test  station  assignment. 

Two  queues  were  established  as  indicated  in  Figure  5.1. 
A  product  entered  the  assignment  queue  if  it  was  scheduled 
to  arrive  at  Systems  Test  at  time  t,  the  present  time 
period.  All  products  were  assumed  to  arrive  at  1200  hours 
on  the  arrival  date.  This  queue  was  processed  during  first 
shift  tine  periods,  with  a  product  leaving  only  after 
assignment  of  a  test  station  and  a  first  shift  testor.  The 
assignments  were  made  with  regard  to  availability,  maximum 
efficiency  and  qualification.  If  required,  the  product 
woulo  also  enter  the  second  shift  testor  queue.  As  with  the 
assignment  queue,  second  shift  testors  are  assigned  with 
regard  to  availability  and  qualification.  This  queue  was 
processed  only  during  second  shift  time  periods.  Time  in 
the  simulator  is  clock  oriented  with  a  fixed  increment  of 
advance.  This  increment  was  chosen  as  one  hour.  The  number 


of  tine  periods  in  a  day  is  set  by  the  user  and  represents 
the  number  of  hours  worked  by  first  and  second  shift 
together.  Work  days  are  treated  as  sequential  with  the 
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clock  running  continuously .  A  listing  of  the  program  is  at 
Appendix  L. 

Evaluation  of  the  Model 

The  evaluation  of  the  model  was  broken  into  two  phases, 
the  conceptual  phase  and  the  implementation  phase.  In  the 
conceptual  phase  the  logical  flow  of  the  model  as  presented 
in  Figure  5.1  was  reviewed  by  the  manager  of  Systems  Test. 
Once  his  confirmation  that  the  model  accurately  reflected 
the  assignment  process  was  received,  program  formulation 
began . 

In  the  implementation  phase,  the  validity  of  the  model 
was  tested  with  the  use  of  actual  schedules.  The  comparison 
involved  the  number  of  tested  hours  and  consistency  of 
assignment.  It  was  determined  that  the  simulation  model  was 
a  realistic  representation  of  the  system  as  number  of  tested 
hours  were  equal  and  there  were  no  breaks  in  assignments 
until  testing  of  a  product  was  complete. 

Phase  II.  Input  Data  Analysis 

The  following  four  steps  are  essential  in  the  develop¬ 
ment  of  a  valid  model  for  input  data  and  constituted  the 
basis  for  the  input  data  analysis: 

Step  1.  Collection  of  raw  data. 

Step  2.  Identification  of  the  underlyina 
statistical  distribution. 


A 
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Step  3.  Estimation  of  parameters  that 
characterize  the  distribution. 

Step  4.  Test  of  the  distributional  assumption 
and  the  associated  parameter  estimates  for 
goodness  of  fit. 

It  was  determined  through  analysis  of  the  system  that  data 
on  the  following  areas  needed  to  be  collected  and  analyzed: 
product  mix,  arrival  times,  interarrival  times,  processing 
times,  and  resource  availability  (testors  and  test 
stations) . 

Product  Mix 

Data  on  the  first  of  these  areas,  product  mix,  was 
obtained  from  the  planning  department  in  the  form  of  a  12 
month  rolling  sales  plan,  Appendix  M.  This  plan  identifies 
by  month  the  number  of  products  by  series  that  are 
forecasted  for  the  year.  The  specific  voltage  configuration 
of  the  products  was  not  addressed.  By  totaling  the  number 
of  each  product  type  through  the  12  month  projection  and 
dividing  by  the  total  number  of  projected  products  for  the 
period,  a  probability  was  associated  with  each  product  type. 
It  should  be  noted  that  projections  are  based  on  sales  and 
that  no  dependence  exists  between  products.  Taking 
advantage  of  these  probabilities  a  valid  product  mix  model 
is  arrived  at  through  the  use  of  a  uniform  distribution  ard 
the  cumulative  distribution  function  for  the  12  month 


projection,  reference  Table  5.1. 


As  stated  previously,  the  voltage  configuration  data 
was  not  a  part  of  the  12  month  projection.  This  information 
could  be  collected  from  historical  data  and  fitted  to  a 
suitable  distribution.  It  is  assumed,  however,  that  the 
voltage  configuration  follows  a  uniform  distribution  and, 
therefore,  each  feasible  voltage  configuration  has  an  equal 
probability  of  occurring. 


Table  5.1  12  Month  Projection 


Product  Type  For 

Projected 

12  Month  Period 

Probability 

Cumulative 

Probability 

2015 

60 

.08 

0.080 

2030 

134 

.  179 

0.259 

2045 

93 

.124 

0.380 

2715 

27 

.036 

0.419 

2730 

4 

.005 

0 .424 

Sub-Total 

2000 

318 

.424 

3060 

77 

.103 

0.526 

3100 

70 

.093 

0.620 

3180 

33 

.044 

0.664 

3250 

24 

.032 

0.696 

3330 

30 

.04 

0 .736 

3450 

76 

.101 

0.837 

3600 

12 

.016 

0.853 

Sub-Total 

3000 

322 

.429 

4080 

41 

.055 

0.908 

5060 

22 

.029 

0 .937 

5100 

12 

.  026 

0.953 

5160 

9 

.012 

0.965 

5200 

12 

.016 

0.981 

5300 

13 

.017 

0.998 

5400 

1 

.001 

1  .000 

Sub-Total 

5000 

69 

.092 

TOTAL 

750 

i 


It  should  be  noted  that  the  Hertz  description  of  the 
product  is  a  function  of  series  and  voltage  and,  in  essence, 
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has  already  been  accomplished .  In  addition  to  all  5000 
series  products,  all  2000  series  products  with  an  input 
voltage  of  400  or  415  volts  have  a  50  Hertz  rating.  All 
4000  series  products  have  a  415  Hertz  rating  and  all  other 
products  have  a  60  Hertz  rating.  This  completes  the  input 
model  for  the  product  mix  as  all  technical  data  required  tc 
completely  describe  the  product  has  been  formulated.  This 
section  is  summarized  by  the  partial  flow  diagram  presented 
in  Figure  5.2. 

Processing  Times 

In  the  development  of  a  valid  input  model  for 
processing  times,  data  was  collected  on  the  actual 
processing  times  of  149  products.  In  order  to  limit  the 
collection  period  it  was  decided  that  products  should  be 
grouped  in  accordance  with  the  standard  allotted  test  hour 
scheme,  reference  Appendix  F. 

2000  Series 

For  the  2000  series  UPS,  data  was  collected  for  60 
products,  presented  in  Table  5.2.  A  frequency  distribution 
or  histogram  was  then  constructed,  reference  Figure  5.3,  in 
order  to  infer  a  known  probability  distribution  function. 
Based  on  the  shape  of  the  histogram,  a  normal  distribution 


was  assumed. 


Partial  Flow  Chart  for  Product  Mix 


G  2 


If  the  observations  in  a  sample  of  size  n  are  x,  , 

X2''‘',Xn'  t^ie  samP^e  mean  (x)  is  defined  by 

n 

y  x. 

-  i=i  1 
X  =  - 

n 

2 

and  the  sample  variance  (S  )  is  defined  by 

n  __  -» 

Z  ex.-xr 

s2  =  - 

n- 1 

For  the  60  products  of  Table  5.2,  x=38.827  and  S2=2 59.92 
(S=16. 122) . 

In  applying  a  goodness-of-f it  test  the  Geary  Test  of 
Normality  was  used.  This  procedure  is  based  on  the  ratio  of 
the  average  absolute  deviation  to  the  square  root  of  the 
average  squared  deviation.  A  random  variable  U  is  defined 
by : 

vT/2  Z I x .  -  x  I  /n  1 . 2533  Z  I  x .  -  x  |  /n 

11  = - : -  =  ■  ■  ■■■  - — — 

^ ZOj  -  x)“/n  J  y(xi  -  x)“/n 

When  the  underlying  distribution  is  normal,  both  the 
numerator  and  denominator  of  U  estimate  C,  and  the  expected 
value  of  U  is  approximately  1.  A  departure  from  normality 
is  indicated  by  a  value  of  U  which  differs  considerably  from 
1.  The  resulting  hypotheses  and  corresponding  test 
statistic  are: 

Hq:  the  underlying  distribution  is  normal 

H^:  the  underlying  distribution  is  not  normal 

ll  -  1 
.  .'(-(.I  /  n 


Test  Statistic: 


Rejection  Region:  Hither  Z  >  2  or  :  <*2  .... 

J  3  —  a/2  -  a/2  , 

In  computing  U,  it  can  be  noted  that  £|x.  -x|  =  2(£x.-n  x), 
where  n1  denotes  the  number  of  X^'s  which  exceed  X  and  £' 
is  the  sum  of  those  n'  observations. 

For  the  specific  problem;  X=38.827,  n=60,  n'=29  and 

f 

£  xi=1509.8. 


x. 


xl/n  = 


:(£ 


x^  -  n  x) /n  = 


2(1509. 8  -  29(38. 82") 
60 


)  =  12. "959 


- 1 

M 

/ - N 

X 

H* 

1 

_ iil 

1/2 

1409.435 

L  n  _ 

60 

15.710651 


(1.2535) (12.7959) 
15.710631 


1.0210587 


1.0210587  -  1 
. 2661 / /oO 


.  6130026 


Using  a  confidence  level  of  95%,  a/2=.025,  Z  q0^=1.96. 
Since  Z  <zaj2'  t^le  assumPtion  that  the  underlying  distri¬ 
bution  is  normal  is  accepted. 

3000  Series  <.100  Kw 


The  procedure  for  arriving  at  the  input  model  for  the 
3000  Series  UPS  with  Kw <  100  was  the  same  as  the  procedure 
for  the  2000  Series  UPS.  A  sample  size  of  35  products  was 
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taken  with  actual  processing  durations  exhibited  in  Table 

5.3. 

The  resulting  histogram  is  Figure  5.4.  Based  on  the 
histogram  a  normal  distribution  was  assumed. 


Table  5.3  Actual  Processing  Times  for  3000  Series  UPS 
with  Kw  <  100 


34 

19.5 

49 

55 

41  .5 

41.5 

39 

79 

39.5 

5 

30 

5.5 

77.5 

41 

3 

46 

24 

5 . 5 

52.5 

22.5 

36 

41.5 

54 . 5 

14.5 

17.5 

64 

24 . 5 

13.5 

25 

9 

57 

16 

19.5 

16.2 

2  - 

The  estimates  of  the  parameters  u  and  a  are  X=32.918  and 
2 

S  =417.0,  respectively. 


Figure  5.4.  Histogram  for  3000  Series  K  ^_100 

Again,  using  Geary's  Test  of  Normality: 

H  :  the  underlying  distribution  is  normal 
o  J 

Ha:  the  underlying  distribution  is  not  normal 


For  the  specific  problem;  X=32.918,  n=3 4,  n'=17, 
Sx  .=848.5. 


£lxi  -  x|  2(£  xi  -  n  X1 


2  (848.5  -  1" (52.918) 
34 


=  16.994 


1/2 

13761. 009 

n 

_  _ 

34 

1/2 


=  20. 118 


_  ( 1 . 2553) ( 16.  994 )  _ 

20.118 

1.059  -  1 


1 . 349 


.2661/* 31 


Using  a  confidence  level  of  95%,  a/2=.025,  2  Q25=1.96. 

Since  Z<  za/2'  t^'e  hypothesis  that  the  underlying  distri¬ 
bution  is  normal  is  accepted. 


3000  Series,  180lKw£450 

A  sample  size  of  30  products  was  taken  with  actual 
processing  durations  exhibited  in  Table  5.4.  The  resultin 
histogram  is  Figure  5.5.  Based  on  the  shape  of  the  histo¬ 
gram  a  normal  distribution  was  assumed.  The  estimates  of 

2  —  2 

the  parameters  u  and  a  are  X=38.157  and  S  =246.505, 
respectively . 

Using  Geary's  Test  of  Normality: 

Hq:  the  underlying  distribution  is  normal 

the  underlying  distribution  is  not  normal. 


t 
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Figure  5.5  Histogram  for  3000  Series  180£Kw_<450 

3000  Series,  Kw=600 

A  sample  size  of  5  products  was  taken  with  actual 
processing  durations  exhibited  in  Table  5.5. 

Table  5.5  Actual  Processing  Times  for  3600  UPS 

21  56.6  88 

9  78.5 

The  sample  average,  X,  is  equal  to  50.6  The  standard 

allotted  test  hours  for  this  type  product  is  50.  Due  to  the 

lack  of  sufficient  data,  the  assumption  will  be  made  that 

the  underlying  distribution  is  normal  and  that  u=50.  The 
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standard  deviation,  ,  will  be  assumed  to  be  equal  to  17. 
This  figure  falls  between  the  maximum  and  minimum  standard 
deviations  previously  observed. 


Actual  test  data  was  aviiilable  for  only  two  4000  scries 
UPS.  As  with  the  3600  UPS,  the  assumption  will  be  made  that 
the  underlying  distribution  for  the  processing  durations  of 


the  4000  series  UPS  is  normal.  The  mean  is  assumed  to  be 
equal  to  the  standard  allotted  test  hours,  u=35,  and  the 
standard  deviation  is  assumed  to  equal  17. 

5000  Series 

A  sample  size  of  18  products  was  taken  with  actual 
processing  durations  exhibited  in  Table  5.6.  Due  to  the 
lack  of  sufficient  data,  the  assumption  will  be  made  that 
the  underlying  distribution  is  normal.  Since  the  sample 
mean,  X=36.294,  differs  significantly  from  the  standard 
allotted  test  hours,  50,  the  sample  mean  and  sample  standard 
deviation,  S=13.909,  will  be  used  as  estimators  for  u  and 
a“  ,  respectively. 

Table  5.6  Actual  Processing  Times  for  5000  Series  UPS 


6 

50.5 

55 

26  .  3 

32 

36.5 

34.5 

31 

20 

62.5 

35 

57 

3  6 

32.8 

45.5 

23 

37 

3  2.7 

Table  5.7  presents  a  summary  of  the  parameters 
resulting  from  the  processing  time  input  analysis. 
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Table  5.7  Parameters  for  Processing  Times 


Product 

U 

0 

2000 

38.827 

16.122 

3000  KwA  100 

32.918 

20.421 

3000  1  80  <_  Kw  £  450 

38.15'’ 

15.701 

3000  Kw= 600 

50 . 0 

17.0 

4000 

35.0 

17.0 

5  00  0 

36.294 

13.909 

All  processing  times  were  found  or  assumed  to  fit  a  normal 
distribution.  The  distributions  will  be  trunctated  to  the1 
left  of  the  value  3,  the  minimum  observed  test  time. 
Processing  time  generation  is  illustrated  in  Figure  5.6. 


Figure  5.6  Processing  Time  Generation 

Arrival  and  Interar rival  Times 

All  attempts  to  collect  data  on  arrival  and  inter¬ 
arrival  times  failed.  It  was  observed,  however,  that  the 
maximum  number  of  UPS  completed  in  a  month  was  93  and  that 
on  the  average  three  to  four  products  arrived  each  day.  It 
has  also  been  shown  lhat  in  numerous  situatiors,  such  as  the 
arrival  of  jobs  at  a  job  shop,  arrivals  occur  according  to  a 
Poisson  process.  Although  Systems  Test  is  not  a  job  shop, 
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numerous  similarities  exist.  Checking  the  assumptions  asso¬ 
ciated  with  a  Poisson  process  the  following  was  observed: 

1.  Products  are  completed  by  the  production  line  one 
at  a  time.  Since  this  is  the  source  of  products  for 
systems  Test,  products  can  be  assumed  to  arrive  at 
Systems  Test  one  at  a  time. 

2.  Arrivals  are  completely  at  random  throughout  the 
first  shift  without  rush  or  slack  periods. 

3.  A  large  or  small  number  of  arrivals  on  one  shift 
had  no  effect  on  the  number  of  arrivals  in  subsequent 
shifts . 

Arrivals  to  Systems  Test,  therefore,  met  the  assumptions 
associated  with  a  Poisson  process.  Accordingly,  a  Poisson 
distribution  with  X  =4,  will  be  assumed  to  model  arrivals. 

The  relationship  between  arrival  date  and  due  date  was 
then  investigated  with  the  use  of  the  PM01  Report,  Appendix 
D.  Two  dates  were  of  interest;  the  Original  Operations 
Promise  Date  (Appendix  D,  #8),  and  the  Original  Scheduled 
Start  for  Systems  Test  (Appendix  D,  #1) .  By  taking  the 
difference  of  these  two  dates  and  subtracting  nonworking 
days  (weekends,  holidays,  etc.)  and  the  standard  lead  time 
for  packing  and  shipping,  a  distribution  for  planned  test 
window  times  was  developed  for  each  product  category  as  with 
processing  times. 
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3000  Series,  Kw  £  1 0 0 

A  sample  size  of  53  products  was  taken  with  planned 
test  window  times,  calculated  as  described  previously, 
exhibited  in  Table  5.8.  The  resulting  histogram  is  Figure 
5.7.  An  attempt  was  made  to  fit  the  data  to  a  number  of 
distributions,  all  of  which  failed  the  goodness-of-f it  test. 
Based  on  the  shape  of  the  histogram  and  the  fact  that  lead 
times  often  follow  a  gamma  distribution  the  determination 
was  made  to  fit  a  gamma  distribution  to  the  data. 


Table  5.8  Test  Window  Times  for  3000  Series 
UPS  with  Kw  <  100 
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Figure  5.7  Test  Window  Times  Histogram  for 
3000  Series  UPS  with  Kw  ^100 

The  gamma  distribution  has  the  following  density 

function : 


rn r  r-1  -0 

i  0  y  e 
fy(y)  =  J  r(r) 

L° 


if  y  >  0 


if  v  <  0 


Where  0  and  r  are  positive  numbers  and  r (r)  is  the  gamma 
function . 

To  find  values  of  0  and  r  from  the  data  in  Table  5.8 
the  mean  and  sample  variance  must  first  be  calculated. 


n 

V  X. 
i= i  1 


11  9 

£(xrGYr 


°Y  = 


i=  1 


10. 925 


?8. 447 


The  parameters  ©  and  r  are  now  determined: 


0  =  —  = 


10. 925 
28.447 


=  0.584 


r  = 


(uYr 


Y 


(10.925)' 

28.447 


=  4.196 


3000  Series,  180  1  Kw  1450 

A  sample  size  of  41  products  was  taken  with  planned 
test  window  times  exhibited  in  Table  5.9.  The  resulting 
histogram  is  Figure  5.8.  Again  the  data  was  not  found  to 
fit  a  standard  distribution  and  it  was  determined  that  a 
gamma  distribution  would  be  fit  to  the  data. 

To  find  values  of  ©  and  r  from  the  data  in  Table  5.9, 
the  mean  and  sample  variance  were  first  calculated  with 
u  =11.935  and  ^  =24.206.  The  values  of  ©  and  r  were  then 

y  y 

determined  to  be  0.493  and  5.885,  respectively. 


Table  5.9  Test  Window  Times  for  3000  Series  UPS 
with  180  <  Kw  <  450 
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Figure  5.8  Test  Window  Times  Histogram  for  3000  Series 
UPS  with  180  <  Kw  <  450 


3000  Series,  Kw=600 

A  sample  of  seven  products  was  taken  with  planned  test 
window  times  exhibited  in  Table  5.20.  The  resulting 
histogram  is  Figure  5.9.  The  mean  and  sample  variance  for 
the  data  are  X-10.643  and  S^=25.807.  Due  to  the  lack  of 
sufficient  data  a  normal  distribution  will  be  assumed. 
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Table  5.10  Test  Window  Times  for  3600  UPS 


4.5 

11.5 

6 . 5 

17 . 5 

16.6 

11.5 

6.5 

4000  Series 

A  sample  of  15  products  was  taken  with  planned  test 
window  times  exhibited  in  Table  5.11.  The  resulting 
histogram  is  Figure  5.10.  The  mean  and  sample  variance  for 
the  data  are  X=12.0  and  S  =24.711.  Based  on  the  shape  of 
the  histogram  and  the  lack  of  sufficient  data  a  normal 


distribution  is  assumed. 
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Table  5.11  Test  Window  Times  for  4000  Series  UPS 
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Figure  5.10  Test  Window  Times  Histogram  for  4000  Series  UPS 

5000  Series 

A  sample  size  of  30  products  was  taken  with  planned 
test  window  times  exhibited  in  Table  5.12.  The  resulting 
histogram  is  Figure  5.11.  Based  on  the  shape  of  the 

histogram  a  normal  distribution  was  assumed.  The  estimates 

2  _ 

of  the  parameters  u  and  o~  are  X=11.833  and  S2=19.395, 
respectively .  Using  Geary's  Test  of  Normality,  as  described 
in  the  processing  time  section: 
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H  :  the  underlying  distribution  is  normal 
o 

H  :  the  underlying  distribution  is  not  normal 

cl 

For  the  specific  problem;  X=11.833,  n=30,  n  '  =  1 1  ,  £x  .  =1  7' 
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Table  5.12 
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Figure  5.11  Test  Window  Times  Histogram  lor  5000  Series  IT 


Using  a  confidence  level  of  95%,  a  /2=.025,  Z 

Since  Z  >-za//2'  t*le  hypothesis  that  the  distribution  is 
normal  is  accepted. 

2000  Series 

Only  two  products  were  available  for  investigation. 

Each  of  these  products  had  a  planned  test  window  time  of 

eight  days.  Due  to  the  lack  of  sufficient  data  a  normal 

2 

distribution  is  assumed  with  u=8  and  o  =25.  The  value  for 
the  variance  was  chosen  to  fall  within  the  observations 
taken  from  the  other  product  categories. 

A  link  between  planned  test  window  times  and  actual 
test  window  times  must  now  be  developed.  During  the  problem 
definition  phase,  it  was  determined  that  approximately  one 
half  of  the  products  arrive  two  to  four  days  late  and  that 
the  other  half  arrive  on  time  or  one  to  two  days  early.  A 
uniform  distribution  will  be  used  to  determine  if  the 
arriving  UPS  is  early  or  late  with  0.5  being  the  decision 
point.  The  number  of  days  that  an  early  product  arrives 
early  will  be  drawn  from  a  Poisson  distribution  with  A  =1. 
The  number  of  days  that  a  late  product  arrives  late  will  be 
drawn  from  a  Poisson  distribution  with  >  =3. 

It  was  also  determined  in  the  problem  definition  phase 
that  a  carryover  of  products  that  were  not  tested  in  the 


previous  month  contributes  significantly  to  monthly  commit¬ 
ment.  In  order  to  simulate  this  carryover,  products  will  be 
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generated  for  two  months  with  statistics  being  taken  on  the 
scheduling  that  occurs  in  the  second  month. 

Table  5.13  presents  a  summary  of  the  distributions  and 
parameters  resulting  from  the  arrival  time  input  analysis. 
Figure  5.12  illustrates  the  generation  of  arrival  dates  and 
due  dates. 

Table  5.13  Distributions  and  Parameters  for  Planned  Test 
Window  Times 


Product 

Distribution 

u 

2 

a 

e 

r 

2000 

Normal 

8.0 

25.0 

3000, 

Kw  100 

Gamma 

10.925 

28.447 

0.384 

4.  196 

3000, 

180  Kw 

450  Gamma 

11.935 

24.206 

0.493 

5.885 

3000, 

Kw=600 

Normal 

10.643 

25.806 

4000 

Normal 

12.0 

24.711 

5000 

Normal 

11.833 

19.393 
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Resource  Availability 

At  the  time  of  this  study,  test  station  unavailability 
was  negligible.  Future  plans,  however,  call  for  a  periodic 
maintenance  schedule  based  on  six  month  periods  which  should 
coincide  with  plan  shutdowns.  Again,  this  unavailability 
period  would  be  statistically  negligible  since  testing  would 
not  be  scheduled  during  the  period.  Specific  data  was  net 
available  for  testor  unavailability.  However,  it  was  deter¬ 
mined  that  major  periods,  such  as  vacations,  occur  only 
during  plant  shutdown.  Employee  unavailability,  therefore, 
will  be  assumed  negligible. 

Phase  III.  Experiment  Design 

Prior  to  the  discussion  of  the  experiment  design  a 
brief  overview  of  experimentation  and  a  number  of  pertinent 
definitions  should  be  covered.  Of  course,  experimentation 
begins  by  formulating  a  number  of  research  hypotheses.  The 
next  stage  is  the  selection  of  an  appropriate  experimental 
design  within  which  to  maximize  the  information  gained,  with 
respect  to  the  hypotheses,  from  the  experiment.  The  basic 
requirements  of  an  experiment  are  simple:  Differential 
treatments  are  administered  to  different  groups  of  subjects 
(or  to  the  same  subjects  in  different  orders)  and  perform¬ 
ance  on  some  response  measure  is  observed  and  recorded  fol¬ 
lowing  the  administration  of  the  treatments.  The  input 
variables  (independent  variables)  to  the  system,  such  as  the 
decision  variables,  the  structural  assumptions,  and  the 
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parameters  of  the  random  variables,  are  called  factors. 
Factors  may  be  classified  as  qualitative  or  quantitative. 
Quantitative  independent  variables  are  variables  that  repre¬ 
sent  variation  in  amount,  such  as  the  number  of  parallel 
servers.  Qualitative  independent  variables,  on  the  oth-  r 
hand,  represent  variations  in  kind  or  type,  such  as  queue 
discipline.  Each  possible  value  of  a  factor  is  called  a 
level  of  the  factor.  A  combination  of  factors  all  at  a 
specified  level  is  called  a  treatment.  When  a  simulation  is 
run  with  the  same  treatment  but  an  independent  stream  of 
random  numbers  is  used,  it  is  said  that  an  independent  rep¬ 
lication  of  the  experiment  has  been  made.  The  response 
variables  (dependent  variables)  are  those  measures  that  seem 
to  "capture"  the  phenomenon  being  studied  most  accurately. 

Two  factors  of  the  Systems  Test  scheduling  problem  have 
been  chosen  for  investigation;  Queue  discipline  (dispatch 
rules)  and  criterion  for  second  shift  testing  consideration. 
The  first  of  these  factors,  queue  discipline,  was  chosen 
with  the  expectation  that  the  sequential  simulation  will  be 
used  as  a  manager's  desk-side  reference  in  the  future. 
Therefore,  it  is  necessary  that  ar.  exhaustive  test  and 
evaluation  of  the  dispatch  rules  be  conducted  in  order  that 
the  sequential  simulation  package  contain  the  most  suitable 
dispatch  procedure.  It  was  determined  that  this  factci 
would  have  the  following  four  levels: 

Level  1.  First-Come-Fi rst-Cons idt re  ; :  I! 

orders  the  queue  with  respect  t  <  • 
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the  queue  by  placing  the  new  arrival  at  the  end  of 
the  queue. 

Level  2.  Earliest  Due  Date:  Products  are  placed 
in  the  queue  by  due  dates  with  the  product  having 
the  earlier  due  date  being  placed  ahead  of  one 
with  a  later  date. 

Level  3.  Shortest  Processing  Time:  Products  are 
ordered  with  respect  to  the  minimum  total 
processing  time  remaining  at  the  time  of  arrival 
at  the  queue. 

Level  4.  Minimum  Slack  Time:  Minimum  slack  time 
is  defined  as  the  time  remaining  between  the  due 
date  and  the  time  at  which  the  job  could  be 
completed  without  delay.  Priority  jobs  are  those 
jobs  whose  slack  is  negative.  Slack  is  defined 
for  a  job  as  the  due  date  minus  the  current  time 
minus  the  processing  time  remaining. 

The  objective  of  this  experiment  will  be  to  select  that 
dispatch  rule  which  has  the  most  positive  effect  in  reducing 
lateness  and  the  number  of  late  jobs.  It  is  recognized  that 
the  analysis  of  these  two  response  variables  may  yield  two 
different  optimizing  rules.  The  experiment  will  be  a 
single-factor  experiment  that  will  be  completely  randomized. 
This  implies  that  the  replications  of  the  model  at  each 
level  and  for  different  levels  of  the  factor  are  based  on 
independent  streams  of  random  numbers. 
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The  statistical  model  for  the  analysis  of  this  experi¬ 
ment  is 


Yrj=u+Tj+£rj  , r=l , 2 , . . . ,Rj ;  j=l,2,...,k, 

where  ^  is  observation  r  of  the  response  variable  for 

level  j  of  the  factor,  u  is  the  overall  mean  effect,  t_.  is 

the  effect  due  to  level  j  of  the  factor,  e  .  is  a  "random 

J  n  j 

error"  in  observation  r  at  level  j,  and  R.  is  the  number  of 

D 

observations  made  at  level  j.  This  model  is  called  a  fixed 
effect  model.  In  this  experiment  R^=4  and  k=4. 

The  initial  analysis  of  a  single-factor  fixed-effects 
completely  randomized  experiment  consists  of  a  statistical 
test  of  the  hypothesis 

Ha 

that  the  levels  of  the  factor  have  no  effect  on  the 
response.  The  applicable  statistical  test  is  a  one-way 
analysis  of  variance  (ANOVA) .  If  the  computed  value  of  the 
F  statistic  is  not  significant  the  analysis  will  be 


Tj— o  (j— l,2...,k) 


terminated  because  no  differences  between  the  u.'s  have  been 

3 

identified.  But  if  Hq  is  rejected  further  analysis  will 
follow  a  multiple  comparisons  procedure  such  as  the  range 
test. 


If  an  improvement  of  10%  in  lateness  or  the  number  of 
late  jobs  is  observed,  when  comparison  is  made  with  the 
present  system,  a  recommendation  for  employment  will  be 
made.  The  present  system  employs  Earliest  Due  Date  with 
ties  broken  by  kilowatt  rating  as  a  dispatch  rule. 
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The  last  factor  to  be  investigated  is  the  criterion  for 
second  shift  testing.  It  is  felt  that  the  current  criterion 
is  too  restrictive  and  that  by  expanding  this  criterion 
second  shift  testor  utilization  will  increase  and  lateness 
and  the  number  of  late  jobs  will  decrease.  It  was 
determined  that  the  factor  would  have  the  following  levels: 

Level  1:  All  products  that  will  be  late  if  tested 
only  on  first  shift. 

Level  2:  All  products  that  will  be  completed 
within  eight  yours  of  their  due  date. 

Level  3:  All  products  that  will  be  completed 
within  16  hours  of  their  due  date. 

Level  4:  All  products  that  will  be  completed 
within  24  hours  of  their  due  date. 

Level  5:  All  products  that  will  be  completed 
within  32  hours  of  their  due  date. 

Level  6:  All  products  that  will  be  completed 
within  40  hours  of  their  due  date. 

The  objective  of  this  experiment  will  be  to  select  that 
criterion  which  has  the  most  positive  effect  in  reducing 
lateness  and  the  number  of  late  jobs  and  increasing  second 
shift  testor  utilization.  The  experiment  will  be  a  single¬ 
factor  fixed-effects  completely  randomized  experiment  and 
will  follow  the  procedure  of  the  experiment  with  queue  dis¬ 
cipline  with  Rj=4  and  k=6. 

If  an  improvement  of  10%  in  the  lateness  or  the  number 
of  late  jobs  is  observed  a  recommendation  will  be  made  for 


86 


employment  of  that  specific  criterion.  The  present  system 
employs  the  criterion  that  if  the  product  will  be  late  when 
processed  only  on  first  shift,  it  will  be  considered  for 
second  shift  testing. 


Phase  IV.  Output  Analysis 

It  should  be  recalled  that  the  null  hypothesis  in 
single  factor  ANOVA  is  that  the  different  levels  of  the  fac¬ 
tor  have  no  effect  on  the  response: 

Hq:  t j  =  0  ( j  =  l ,  2  ,  .  .  .  ,k)  or  Hq:  u^  =  o>= .  .  .  =u^ 

The  alternative  claim  is  that  there  are  differences  among 
the  responses  for  the  different  levels;  this  can  be  stated 
as 


H  ;  T.^0  (j=l,2...,k)  or  H  :  at  least  two  of  the  u.'s 

are  unequal  1 

The  mean  of  each  level  j  is  computed  by  the  following 
formula: 


R 


Z  Y 

r=l 


rj 


for  j  =  1,2,. . .  ,k 


,  where  Yrj  is  the  response  variable,  r  is  the  replication, 
and  Rj  is  the  total  number  of  replications.  The  dot  in 
place  of  the  second  subscript  in  Y.^  is  used  to  indicate 
sums  over  the  indicated  subscript,  in  this  case  r.  The  ave¬ 
rage  of  all  observed  values  in  the  experiment  is  given  by 


Y 


k 

s 


r=l  j=l 


rj 


R 
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,  where  j  is  the  level  and  R  is  the  total  number  of  replica¬ 
tions  in  the  experiment.  This  is  often  referred  to  as  the 
grand  mean  and  is  equal  to  43.25  in  Table  5.14.  The  ANOVA 
test  is  based  on  partitioning  the  variability  of  the 
observed  responses  Yr j #  into  two  components,  one  component 
due  to  the  level  of  each  factor  and  one  due  to  the  inherent 
variability  of  the  process  being  simulated.  The  first  of 
these  components  is  measured  by  the  treatment  sum  of 
squares : 


SS 


treat 


A  -7  -7 

£  R  Y2  .  -  RY2 
j=l  J 


and  the  second  by  the  error  sum  of  squares: 

SSE  =  £  (Y rj-Y  )2 

b  r=l  j=l  rj  r 

This  leads  to  the  fundamental  identify  of  single-factor 
ANOVA : 


SS  .  =  SS,  ,  +  ss_ 

total  treat  E 

where  SStota^  is  the  total  sum  of  squares: 


SS 


total 


k  R 

=  Z  Z 

j=l  r=l 


—2 
-  RY 


If  the  assumption  of  a  common  variance  is  correct,  the  mean 
square  MS£  =  SSE/  (R-k)  is  an  unbiased  estimate  of  the  vari¬ 
ance,  a2,  of  the  response  variable  Y.  If,  in  addition,  the 
null  hypothesis  is  correct,  the  treatment  mean  square 
MStreat  =  sstreat/ (k-1)  is  also  an  unbiased  estimate  of  o2. 
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In  any  case,  MStreat  anc*  MSE  are  statistically  independent. 
When  H  is  true,  SS,  /  a2  and  SS_/  o2  have  Chi-square 

O  tTcat  ij 

distributions  with  (k-1)  and  (R-k)  degrees  of  freedom, 
respectively.  The  test  statistic  for  testing  the  hypothesis 
is  computed  by: 


F 


MS 


treat 


MS. 


SS 


treat 


/ (k-1) 


SSE/ (R-k) 


When  H  is  true,  this  test  statistic  has  an  F  distribution 
o 

with  k-1  and  R-k  degrees  of  freedom.  The  ANOVA  test  of  the 

null  hypothesis  is  to  reject  H  if  F  >  F  ,  .  ,  D  ,. 

o  —  a  K— 1 r  K— K 


Queue  Discipline 

As  stated  earlier,  R=16  replications  were  made,  R^  =  4  at 
each  of  the  k=4  levels,  for  the  queue  discipline.  The 
results,  Yr j ,  for  the  number  of  late  products  of  replication 
r  and  level  j,  are  given  in  Table  5.14. 


Table  5.14  Number  of  Late  Products  for  Queue  Discipline 


Number  of  Late 
Level  j  of  the 

Jobs ,  Yr j ,  at 

Queue  Discipline 

Replication , 

j  =  l 

j  =  2 

j  =  3 

j  =  4 

Total 

r 

FCFC 

EDD 

SPT 

MST 

Mean 

1 

47 

47 

45 

51 

2 

49 

45 

37 

54 

3 

47 

34 

36 

38 

4 

47 

42 

33 

40 

Totals 

190 

168 

151 

183 

692 

Means 

47.5 

42 

37.75 

45.75 

43.25 
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Table  5.15  ANOVA  for  Number  of  Late  Products  with  Queue 
Discipline 


Source  of 

Sum  of 

Degrees  of 

Mean 

F 

Variation 

Squares 

Freedom 

Squares 

Treatment 

SS.  =224 . 5 

treat 

k-l=3 

MStreat=74*833 

2.437 

Error 

SS  =368.5 

E 

R-k=l 2 

MS  =30.708 

E 

Total 

SStotal=593 

R-l  =  l  5 

The  ANOVA  for  the  number  of  late  products  is  shown  in  Table 

5.15.  With  a  95%  confidence  level,  a=.05, 

F  3.49.  Since  the  test  statistic,  F=2.437  ,  is 

a,k-l,R-k 

less  than  the  critical  value,  the  null  hypothesis  of  no 
treatment  effect  is  accepted. 

The  results  for  the  mean  lateness  for  replications  r 
and  level  j  are  given  in  Table  5.16.  The  ANOVA  for  mean 
lateness  is  shown  in  Table  5.17. 


Table  5.16  Mean  Lateness  for  Queue  Discipline 


Number 

of  Late 

Jobs,  Yr j , 

at 

Level  j 

of  the 

Queue  Discipline 

Replication , 

j  =  l 

j  =  2 

j  =  3 

j  =  4 

Total 

r 

FCFC 

EDD 

SPT 

MST 

Mean 

1 

105 

46 

210 

43 

2 

148 

90 

171 

86 

3 

61 

50 

141 

38 

4 

91 

49 

155 

42 

Totals 

405 

235 

677 

209 

1526 

Means 

101.25 

58.78 

169.25 

52.25 

95.375 

Table  5.17  ANOVA  for  Mean  Lateness  with  Queue  Discipline 


Source  of 

Sum  of 

Degrees  of 

Mean 

Variation 

Squares 

Freedom 

Squares 

Treatment 

SSk  =34772. 75 

treat 

k-l=3 

MS  =11590.917 

treat 

Error 

SS  =9433 

E 

R-k=12 

MS  =786.083 

h 

Total 

SS^  „  =44205.75 

total 

R— 1= 1 5 

Since  the  test  statistic,  F=14.745,  is  greater  than 

F  a, k-1 , R-k  =  F. 05, 3,12  =  3*49'  the  nul1  hyP°thesis  of  no 
treatment  effect  is  rejected  at  the  a =.05  level  of 
significance . 

It  now  remains  to  determine  which  of  the  u.’s  are  dif- 

l 

ferent  from  one  another.  A  method  for  carrying  out  this 
further  analysis  is  called  a  multiple  comparisons  procedure. 
The  procedure  which  will  be  used  in  this  case  is  Tukey's 
Procedure  (the  T  Method).  Tukey's  procedure  involves  the 
use  of  the  statistic  Q  called  a  studentized  range  statistic. 
The  T  method  for  identifying  significantly  different  u^'s 
follows  the  following  sequence  of  steps: 

1.  Select  a  and  find  Q  .  _  . 

X  f  K  f  K  ~  X 

2.  Determine  w  =  QXfkfR_k  "V  MSe/Rj 

3.  List  the  sample  means  in  increasing  order  and 
underline  those  pairs  which  differ  by  less  than  w.  Any 
pair  of  sample  means  not  underscored  by  the  same  line 
corresponds  to  a  pair  of  true  treatment  means  which  are 
judged  significantly  different. 


For  the  specific  problem,  mean  lateness,  a  is  chosen  as  .05, 


°a,k,R-k  Q. 05, 4, 12  “  4'20'  and 
w  =  4.20  ~\/7  86.083/4  =  58.878. 

Arranging  the  four  sample  means  in  increasing  order,  every 
pair  differing  by  less  than  58.878  is  underscored: 

MST  EDD  FCFC  SPT 

y.4  y.2  y.j  y.3 

52.25  55.75  101.25  169.25 

Thus  queue  disciplines  MST,  EDD,  and  FCFC  are  not  signifi¬ 
cantly  different  from  one  another  in  respect  to  mean  late¬ 
ness,  but  are  significantly  lower  than  SPT. 


Preliminary  Conclusions 

The  following  conclusions  are  made  in  respect  to  queue 
discipline : 

1.  The  true  average  number  of  late  products  does  not 
depend  on  queue  discipline. 

2.  Shortest  Processing  Time  increases  the  true 
average  mean  lateness  of  products. 

3.  A  recommendation  will  not  be  made  to  change  queue 
discipline . 

Second  Shift  Criterion 

As  stated  earlier,  R=24  replications  were  made,  R^=4  at 
each  of  the  k=6  levels,  for  the  second  shift  criterion.  The 
results,  Yr j ,  for  the  number  of  late  products  of  replication 
r  and  level  j,  arc  given  in  Table  5.18.  The  ANOVA  for  the 


i 
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number  of  late  products  is  shown  in  Table  5.19.  With  a  95* 
confidence  level,  a  =.05,  F  ^  r_)<=f  05  5  19=2,74,  since 

the  F  statistic,  F=2.621,  is  less  than  the  critical  value, 
the  null  hypothesis  of  no  treatment  effect  is  accepted. 

The  results  for  the  mean  lateness  in  respect  to  the 
criterion  are  shown  in  Table  5.20.  The  corresponding  ANOVA 
is  in  Table  5.21. 


Table  5.18  Number  of  Late  Products  with  Criterion 


Level,  j 

1 

Ref 

jlication 

Totals 

Means 

2 

3 

4 

1 

47 

45 

34 

42 

168 

42 

2 

39 

40 

29 

36 

144 

36 

3 

30 

40 

27 

35 

132 

33 

4 

36 

36 

24 

30 

126 

31.5 

5 

21 

39 

17 

36 

113 

28.25 

6 

30 

31 

23 

31 

115 

28.75 

Total 

798 

Grand 

Mean 

33.25 

Table  5.19  ANOVA  for  Number  of  Late  Products  with  Criterion 


Source  of 
Variation 

Sum  of  Degrees  of 

Squares  Freedom 

Mean 

Squares 

F 

Treatment 

SStreat=53°  k'1=5 

MS.  =106 

treat 

2.621 

Error 

SS  =76  8 . 5  R-k=l 9 

E 

MS  =40.447 

E 

Total 

SStotal=1298-5 
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Table  5.20  Mean  Lateness  with  Criterion 


Replication , 

r 

Level,  i 

1 

2 

3 

4 

Totals 

Means 

1 

46 

90 

50 

49 

235 

58 . 75 

2 

53 

91 

36 

71 

251 

62.75 

3 

45 

88 

48 

67 

248 

62 

4 

48 

99 

48 

69 

264 

66 

5 

60 

63 

51 

47 

221 

55.25 

6 

53 

88 

40 

52 

233 

58.28 

Total 

1452 

Grand 

Mean 

60.5 

Table  5.21 

ANOVA 

for  Mean 

Lateness 

with 

Criterion 

Source  of 

Sum  of 

Degrees  of 

Mean 

F 

Variation 

Squares 

Freedom 

Squares 

Treatment 

SS 

=  293 

treat 

k-1 

=  5 

MStreat=58* 

'  151 

Error 

SS  =7373 

E 

R-k= 

19 

MS£=388 .053 

Total 

SS,  ,  . =7666 
total 

The  critical 

value 

'  F . 05 , 5 , 

19=2.74, 

remains  the  same. 

Since 

the  F  statistic  is 

0.151,  which  is 

less 

than  the  critical 

value,  the  null  hypothesis  of  no  treatment  effect  is 
accepted . 

The  results  for  the  second  shift  testor  utilization 
with  respect  to  the  criterion  are  exhibited  in  Table  5.22. 
These  measures  were  taken  for  a  one  month  period  where  the 
maximum  possible  utilization  figure  would  be  1120.  The 
ANOVA  for  mean  lateness  is  shown  in  Table  5.23.  The  criti¬ 
cal  value  remains  2.74.  Since  the  F  statistic,  1.088,  is 
less  than  the  critical  value  of  the  null  hypothesis  of  no 
treatment  effect  accepted. 
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Table  5.22  Testor  Utilization  with  Criterion 


Ref 

Plication 

,  r 

Level , 

_1 

_2 

_3 

_ 4 

Tota 1 s 

Means 

1 

856 

830 

775 

806 

3267 

816.75 

2 

802 

879 

787 

862 

3330 

832.5 

3 

937 

872 

725 

860 

3394 

848 . 5 

4 

862 

976 

767 

890 

3495 

873.75 

5 

951 

961 

783 

891 

3586 

896 . 5 

6 

977 

918 

805 

933 

3633 

908.25 

Total 

20705 

Grand 

Mean 

862.708 

Table  5.23  ANOVA  for  Testor  Utilization  with  Criterion 

Source  of  Sum  of  Degrees  of  Mean 

Variation  Squares  Freedom  Squares 

Treatment  SS.  =26271 . 512  k-l=5  MS.  ,=5254.302 

treat  treat 

Error  SSE=91767.25  R-k=19  MSE=4B29.855 

Total  SS^.  ^  .=118038.76 

total 

Preliminary  Conclusions 

The  following  conclusions  are  made  in  respect  to  the 
second  shift  criterion. 

1.  The  true  average  number  of  late  products  does  not 
depend  on  the  second  shift  criterion. 

2.  The  true  average  mean  lateness  of  products  does 
not  depend  on  the  second  shift  criterion. 

3.  The  true  average  second  shift  testor  utilization 
does  not  depend  on  the  second  shift  criterion. 

4.  A  recommendation  will  not  be  made  to  change  the 
second  shift  criterion. 


Conclusions 


All  of  the  results  appear  to  be  counter-intuitive.  In 
the  queue  discipline  experiment,  since  there  are  many 
favorable  comparisons  that  can  be  made  between  Systems  Test 
and  the  single  machine  sequencing  problem  with  independent 
jobs,  one  would  expect  that  the  dispatch  rule  of  shortest 
processing  time  would  yield  the  most  favorable  result.  The 
fact  that  the  dispatch  rules  were  not  significantly  differ¬ 
ent  from  one  another  leads  one  to  believe  that  the  reason 
for  the  large  number  of  late  jobs  does  not  rest  with  Systems 
Test.  This  conception  is  strengthened  by  the  results  of  the 
experiment  investigating  the  second  shift  processing  crite¬ 
rion  . 

In  this  experiment,  one  would  expect  the  number  of  late 
jobs  to  decrease  and  second  shift  testor  utilization  to 
increase  as  more  jobs  were  considered  for  second  shift 
processing.  The  fact  that  the  criteria  were  shown  to 
produce  results  not  significantly  different  from  one  another 
implies  that  the  second  shift  testor  resource  is  already 
committed  to  products  which  would  be  late  without  a  safety 
factor.  The  combined  results  suggest  an  investigation  of 
why  products  arrive  late  to  Systems  Test  and  whether  lead 
times  for  product  manufacture  are  sufficient. 

As  a  result  of  the  nonincrease  of  second  shift  testor 
utilization  in  the  experiment  involving  the  second  shift 
processing  criterion,  an  additional  observation  can  be  made. 
The  maximum  possible  second  shift  testor  utilization  for  the 
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period  under  consideration  is  1120.  The  grand  mean  for  this 
experiment  was  only  862.7.  This  result  suggests  a  further 
investigation  into  the  testor  qualifications  of  second  shift 
personnel  as  qualifications  do  not  appear  to  meet  the  needs 
of  the  system.  An  initial  verification  of  this  rationale  is 
obtained  by  checking  the  qualification  matrix  for  second 
shift  operators  shown  below. 


0  0  1 
110 
1  1  0 
1  1  0 


1111 
0  0  0  0 
0  0  0  0 
0  0  0  0 


The  sparseness  of  the  matrix  supports  the  conclusion. 


VI .  Summary  and  Recommendations 

Subsequent  to  a  complete  problem  definition,  an 
analytical  scheduling  model  was  developed  to  completely 
describe  the  Systems  Test  assignment  procedures  with  regard 
to  limited  resources  and  constraints.  In  determining  the 
computational  complexity  of  the  model,  it  was  determined, 
however,  that  the  scheduling  model  could  not  be  solved  for 
an  optimum  schedule.  This  was  due,  primarily,  to  the  unique 
second  shift  organization  and  the  desire  not  to  process  all 
products  on  second  shift.  By  eliminating  this  consideration 
and  dealing  only  with  a  single  shift,  a  reduced  formulation 
was  developed  which  fell  within  the  constraints  of  a  mixed 
integer  programming  package.  Combining  this  formulation 
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with  an  iterative  algorithm,  a  solution,  although  not  opti¬ 
mum,  was  obtained. 

The  capacity  of  Systems  Test  is  a  function  of  the  spe¬ 
cific  product  mix.  Therefore,  an  analytical  capacity  model 
was  also  developed.  In  an  effort  to  reduce  the  computa¬ 
tional  complexity  of  the  formulation,  the  assumptions  that 
products  were  immediately  available  and  that  testors  were 
readily  available  were  made.  In  this  manner,  time  and 
testors,  unlike  the  scheduling  formulation,  were  not 
parameters  of  the  model.  Due  to  these  assumptions,  however, 
the  solution  is  at  best  an  estimate  of  capacity  for  the 
specific  product  mix.  A  related  heuristic  remains  to  be 
developed  which  will  select  products  for  deletion  from  an 
infeasible  solution.  Computational  complexity  again 
overshadowed  the  usefulness  of  the  formulation  and  it  was 
determined  that  a  simulation  approach  would  be  used. 

The  concept  was  to  employ  simulation  both  as  a  tool  of 
design  and  a  tool  of  operation.  As  an  operational  tool,  a 
sequential  simulation  was  developed  which  is  suitable  for 
use  by  a  manager  as  a  desk-side  reference.  It  assists  the 
manager  in  determining  whether  he  or  she  will  be  able  to 
meet  monthly  test  requirements,  and  in  planning  shift  organi¬ 
zation,  overtime  requirements,  et  cetera.  It  fully  con¬ 
siders  the  facility  resource  constraints,  both  personnel  and 
equipment . 

As  a  design  tool,  the  sequential  simulation  was 
combined  with  an  input  simulation  to  yield  a  randomized 
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simulator.  The  input  simulation  was  the  result  of  the 
collection  and  analysis  of  data  from  the  test  facility  and 
accurately  reflected  the  technical  description  of  the 
products  and  their  arrival  to  Systems  Test.  It  was  assumed 
that  product  arrival  followed  a  Poisson  distribution  with 
X=4  and  that  half  of  these  products  were  late  and  half  early. 
The  "lateness"  was  assumed  to  follow  a  Poisson  distribution 
with  X  =3  and  the  "earliness"  a  Poisson  distribution  with 
X  =1.  Probabilities  were  associated  with  product  types 
based  on  an  actual  12-month  production  projection.  Due 
dates  were  assigned  to  products  based  on  the  analysis  of 
planned  test  window  times  and  whether  the  product  was  deter¬ 
mined  to  be  late  or  early.  In  two  cases  a  gamma  distribu¬ 
tion  was  fit  to  the  collected  data  of  the  planned  test 
window  times,  as  the  data  did  not  fit  a  standard  known  dis¬ 
tribution.  In  all  other  instances,  a  normal  distribution 
was  shown  to  be  appropriate  or  assumed.  A  normal  distri¬ 
bution  was  also  shown  to  be  appropriate  for  product  pro¬ 
cessing  requirements. 

Two  factors  of  the  scheduling  procedure  were  chosen  for 
investigation  with  the  use  of  the  simulator.  The  first, 
queue  discipline,  was  tested  at  four  levels,  and  the  second, 
the  second  shift  processing  criterion,  was  tested  at  six 
levels.  It  was  statistically  determined  that  the  influence 
of  these  factors  on  the  number  of  late  products  was  not 
significant.  It  was  also  determined  that  the  influence  of 
queue  discipline  on  mean  lateness  was  statistically  signifi- 
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cant,  in  that  the  mean  lateness  associated  with  shortest 
processing  time  was  greater  than  all  other  tested  queue 
disciplines.  In  the  case  of  the  second  shift  processing 
criterion,  it  was  shown  that  the  influence  of  the  criterion 
on  mean  lateness  and  second  shift  testor  utilization  was  not 
significant . 

It  is  recommended  that  the  reduced  scheduling  model  and 
the  related  algorithm  be  employed  by  the  company,  both  as  a 
planning  tool  and  an  operational  tool.  As  a  planning  tool, 
it  is  recommended  for  use  as  a  capacity  model  and  the  first 
step  in  a  backward  planning  process.  By  not  considering  the 
planned  arrival  times  of  the  products  and  instead  assuming 
that  products  are  immediately  available  for  processing,  the 
reduced  formulation  and  related  algorithm  could  be  used  to 
determine  capacity  for  specific  product  mixes. 

Additionally,  this  method  could  be  used  to  determine  best 
arrival  times  for  products  to  Systems  Test  and  then,  through 
a  backward  planning  process,  related  manufacturing  times. 

As  an  operational  tool,  it  is  recommended  that  the  reduced 
formulation  and  related  algorithm  be  employed  by  the  manager 
of  Systems  Test  on  receipt  of  the  Monthly  Sales  Plan  in 
order  to  develop  the  manager's  Monthly  Plan  for  Test.  In 
this  way,  the  manager  will  be  better  able  to  allocate  his 
resources,  such  as  second  shift  personnel  and  overtime 
hours . 

The  manager  should  also  employ  the  sequential  simulator 
as  a  desk-side  reference.  Since  it  was  designed  to  be  flex- 
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ible,  he  or  she  can  experiment  on  an  as  needed  basis,  per¬ 
haps  weekly,  with  changes  in  user  defined  variables,  such  as 
shift  durations  and  testor  shift  assignments.  The  simulator 
would  also  provide  an  immediate  scheduling  update  as  infor¬ 
mation,  such  as  new  arrival  times,  becomes  available. 

Finally,  it  is  recommended  that  current  procedures 
with  regard  to  dispatch  rules  and  second  shift  processing 
considerations  do  not  change.  However,  as  other  changes  are 
made  outside  Systems  Test  that  may  affect  arrival  times,  in 
relation  to  due  dates,  it  is  recommended  that  the  simulation 
experiments  be  repeated  with  the  related  changes. 

All  of  the  objectives  of  the  project,  as  detailed  in 
Chapter  III,  and  the  simulation  objectives  of  Chapter  V, 
have  been  accomplished.  Many  aspects  of  the  Systems  Test 
facility  still  deserve  study,  however.  Recommendations  for 
continued  work  are  as  follows: 

1.  As  indicated  in  the  simulation  experiment,  second 
shift  testor  qualifications  may  not  match  the  future 
needs  of  the  system.  A  development  of  an  improved  mix 
of  testor  qualifications  and  numbers  is  warranted. 

This  could  be  accomplished  by  a  single-factor 
fixed-effects  completely  randomized  experiment,  with 
the  developed  simulator,  following  the  procedure 
presented  in  Chapter  V  of  this  study. 

2.  The  simulation  experiment  also  indicated  that  the 
source  of  late  products  lies  outside  Systems  Test.  An 


investigation  of  production  lead  times  would  further 
refine  the  search. 


3.  As  indicated  previously,  further  development  of 
the  capacity  model  and  the  related  heuristic  demands  a 
separate  treatment. 

4.  Plans  are  in  the  offing  at  the  plant  for  expansion 
of  Systems  Test  facilities.  A  complete  investigation 
of  requirements  prior  to  final  determination  would 
maximize  return  of  the  expansion. 

5.  Input  data  collection  should  continue  in  order  to 
increase  the  validity  of  the  simulation,  especially  in 
the  area  of  arrival  and  interarrival  times. 


f 
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Appendix  A 


UPS  Products 


Administrative 

Number 

Series  and 

Kw 

Voltage 

In/Out 

3000  Series 

1 

3025 

204/204 

3 

3030 

208/208 

4 

3030 

480/208 

5 

3030 

480/480 

6 

3050 

204/204 

8 

3060 

208/208 

9 

3060 

480/208 

10 

3060 

480/480 

54 

3060 

220/208 

14 

3100 

208/208 

15 

3100 

480/208 

16 

3100 

480/480 

17 

3180 

208/208 

18 

3180 

480/208 

19 

3180 

480/480 

63 

3200 

208/208 

20 

3250 

208/208 

21 

3250 

480/208 

22 

3250 

480/480 

23 

3330 

208/208 

24 

3330 

480/208 

25 

3330 

480/480 

26 

3400 

208/208 

27 

3400 

480/208 

28 

3400 

480/480 

29 

3450 

208/208 

30 

3450 

480/208 

31 

3450 

480/480 

32 

3600 

480/480 

2000  Series 

33 

2015 

208/208 

34 

2015 

480/208 

35 

201 5 

480/480 

36 

2015 

380/380 

56 

2015 

400/400 

37 

2030 

208/208 

38 

2030 

480/208 

39 

2030 

480/480 

40 

2030 

380/380 

66 

2030 

400/400 

41 

2045 

208/208 

42 

2045 

480/208 

43 

2045 

480/480 

44 

2045 

380/380 

62 

2045 

600/280 

Administrative 

Number 

Series  and 

Kw 

Voltage 
In /Out 

Single  Phase 

45 

2715 

208/120 

46 

2715 

480/208 

60 

2715 

480/120 

61 

2715 

415/120 

64 

2715 

400/230 

47 

2730 

208/120 

48 

2730 

480/120 

49 

2730 

400/400 

57 

2730 

400/230 

50 

2745 

400/400 

4000  Series 

11 

4080 

208/208 

12 

4080 

380/208 

13 

4080 

480/208 

5000  Series 


2 

5025 

380/380 

7 

5050 

380/380 

53 

5060 

380/380 

65 

5060 

400/400 

68 

5060 

415/415 

55 

5100 

380/380 

58 

5160 

380/380 

67 

5300 

380/380 

Parallel  Products 

51 

4080 

480/208 

52 

3450 

480/480 

69 

3400 

480/480 

70 

34  50 

480/480 

71 

3300 

480/480 

72 

3330 

208/208 

t  .. 
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Test  Facility  Capability  Study 


12/30/80 
updated  5/26/83 


NOTES: 

STA 

A2, 

A3,  Bl 

— 

No  load 

capability  for  over  208V 

Limited 

to  60Kw  or  less 

STA 

Bl, 

B2 

- 

No  50Hz 

or  gen.  test  capability 

- 

E .  G .  nu 

1st  article  Navy 

— 

Limited 
one  6  0 

to  one  250Kvv  or  less  plu 

STA 

C,  D 

Limited 

to  330Kw  or  less 

STA 

G 

— 

Limited 

to  480  in  and  out  units, 

except  for  514  Hz=  output  models 
Limited  to  total  of  250  Kw 
Limited  to  400A  generator  (50Hz) 

STA  E ,  F  ,  G  ,  H  ,  I ,  J  _  No  60  or  120  cell  DC  available 

STA  A,  B  -  100,  120  cell  DC  available 

STA  E,  F  -  Limited  to  total  of  one  50Hz=  unit 

at  a  time 

Limited  in  inductive  load  capacity 
(300  kvars  total  in  both  stations) 
Limited  to  one  208V  output  system 
in  either  station 

Limited  to  one  600Kw  unit  in  either 
station 

No  208V  inductive  load  available 


STA  H,  I  No  permanent  208V  supply 

No  inductive  Load 
No  generator  capability 

GENERATOR  -  Limited  to  425  KVA  ?  208  (60Hz)  354 

KVA 

Limited  to  550  KVA  @  480  (60Hz)  450 
KVA 


STA  J 


No  inductive  load 
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Facility  Capability  Study 


5/26/83 


Individual  Station  Notes: 


STA  A1 


STA  A2,  A3 
STA  B 1 

STA  B2 


STA  C 


STA  D 


Maximum  size  unit  -  60  Kw 

Load  #1  -  119.5  Kw  @  208  (max.  V=208) 

Kvar  @  208  (max.  v=208) 
Load  #2  -  300  Kw  @  240  or  480 
(Max.  v=208) 

750  Kvar  @  240  or  480  (miax.  v=480) 
(Notice:  this  Load  #2  is  same  load  as 

B2  and  G  load.  E.g.,  only  one  unit  or 
load  at  a  time . ) 

Voltages  available:  208,  480,  380  (or 
nom.  gen.  voltage.) 

Battery:  60  cell  Ex-13  (Notice:  this 

battery  is  common  with  A2 ,  A3,  B1 ,  B2.6C 
cell  battery,  hence  one  system  in  these 
stations  on  battery  at  one  time.) 

Same  as  A1  except  no  Load  #2 

Same  as  A1  except  no  Load  #2  and  no  380v 
50  Hz  voltage  available 

Max.  size  unit  250Kw 
Load:  same  as  Load  #2  (station  Al) 

Voltages  Available:  208,  480  (No  50  Hz) 
Battery  1:  60  cell  Ex-13  (same  battery 

as  Station  Al) 

Battery  2:  192  cell  Ex-27  (same  battery 

as  C,  D,  E,  F,  G.192  cell,  hence 
only  one  unit  in  these  stations  on 
battery  at  one  time) 

Max.  size  unit:  330Kw 
Load:  475  Kw  @  240  or  480 

350  kvar  @  240  or  480 
(Shared  with  Station  D) 

Voltages  available:  208,  480,  380(50Hz) 
(by  conversion  of  normal  480V  60Hz 
feeder ) 

Battery  1:  60  or  120  cell  FTC-21 

(Shared  with  Station  D  and  30 
charger  test  area,  hence  one  unit 
in  these  station  on  this  battery  at 
one  time . ) 

Battery  2:  same  as  battery  2,  Station 
B2 

Same  as  Station  C 
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STA  E  -  Max.  size  unit  60GKw  (max.  208  unit: 

480  Kw) .  (However,  can  be  used  for 
two  module  parallel  sys.  testing  up 
to  800  Kw . ) 


Load 

#1  ; 

500 

Kw 

480 

V 

only* 

300 

Kvar 

@  480 

v  only* 

Load 

#2: 

400 

Kw 

0 

480 

V 

only* 

Load 

#3: 

800 

Kw 

@ 

240 

or 

■  480* 

♦Shared  with  Station  F 
Voltages  available:  208,  480,  380  50Hz 
(max.  of  375  Kva) 

Battery:  Same  as  Battery  2,  Station  B2 

STA  F  -  Same  as  Station  E 


STA  G1 ,  G2 


STA  H 


STA  I 
STA  J 


Max.  size  unit:  250  Kw  (total  cap.  of 
Gl+G2=250Kw) 

Load:  same  as  Load  #2,  Station  A1 

Voltage:  480,  380  (50  Hz) 

(Only  1  50Hz  unit  at  a  time.) 
Battery:  same  as  Battery  2,  Station  B2 

Max.  size  unit:  600Kw  (can  be  used  for 
up  to  4  modules  in  parallel  at  1200  kw) 
Load  #1:  1200  Kw  @  240  or  480* 

Load  #2:  1200  Kw  @  240  or  480* 

♦Shared  with  Station  I 
Voltage  available:  480  only 
Battery:  192  cell  Ex-33B  (shared  with 

Stations  I,  J) 

Same  as  Station  H 

Max.  size  unit:  600  Kw 

(Up  to  3  units  at  1000  KW  total  in 
parallel  system) 

Load  #1:  500  Kw  @  240  or  480 

Load  #2:  500  Kw  A  240  or  480 

Voltages:  208,  480,  max.  of  2  units  on 

208;  max.  of  3  units  on  480 
Battery:  Same  as  Station  H 


no 


O  O  iiO  O  O 


o  o  o  o  o  o  o  o  oo  x  x  no  oo  oo  o  o  o  o  x  o  o  c* 


Appendix  C.  Project  Management  System  PM  21  Report 
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Appendix  F 


Allotted  Test  Times 


Hrs . 


30 

-  100  kw 

3000 

Series 

30 

180 

-  450  kw 

3000 

Series 

40 

6  00  kw 

3000 

Series 

50 

80 

kw  415Hz 

4000 

Series 

35 

15 

-  4  5  kw 

2000 

Series 

15 

30 

UPS  50Hz 

5000 

Series 

50 

Standard  Factory  Witness  Test 

8 

Custom  Content 

25 

Multimodule , 

Swg 

90 

Voltage  Change  Retest  8 


Appendix  G.  Monthly  Sales  Plan 
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Appendix  H.  Systems  Test  Manager's  Plan 


Appendix  I.  Ups  Production  Plan 
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Appendix  J. 

Test  Station  Efficiencies 


T  i  S  T  5  T  U  T  I  D  N 

*  1  2  3  4  5  6  7  b  91,11121314 

i  r-;jtCT****‘* . . 

1  *11  10*2  "'*C.  *0*0 

2  *  i  :  t  :•  o  2  c  o  o  .  :  o  ?  * 

3  *1112343  0  2,000 

4  *  i  i  i  2  :■  5  o  ;■  o  .  *  z  o 

b  *  2  o  ■’  ■'  1  7  9  P  r  „  9  o  " 

<'  •lii2''>'  i''.  °  " 

7  *  1  0  C  O  0  2  C  U  0  .  3  0  0 

«  *  1  1  1  2  2  3  ,  U  0  ,  0  0  0 

a  «  i  i  i  2  2  .3  :•  0  o  .  '  9  n 

1  '  «  5  i  *  ^  j  ^  '  n  r*  n  j 

11  O 

12  *0  3  3  00  1200,  300 

13  *0*0101*9456,^32 

14  •  0  n  •'  2  2  1  1  „•  *  ,  0  0  1 

lb  *  j  i  "  2  2  1  1  0  r  .  ”  0  0 

j(}  2  i1  "  r>  n 

17  *03  0  441  1  2  2-000 

19  *  *  *  •••  4  4  1  1  2  2  .  j  0  ? 

19  441122  -  635 

2  1  *  ‘  4  -■'«=*  1  2  2  -  0  0  0 

21  *330  0  53122.  300 

2?  *2  30  0  fall  11«  522 

2  3  *  0  *  *  *  0  5  4  l  1  „  *  0  3 

2  4  *  3  ?  "  0  *  5  4  i  i  „  0  3 

9  b  *  Till  nop 

2  b  *  0  3  f  G  0  0  3  1  1  .  C  0  0 

27  *0  0C  00*111,000 

2H  *  3  O  ‘  ?  *  3  3  1  1  .  7  2  2 

op  *ortr*ri'iii."n-ii 

30  *000000011-0001 

31  *000003011.  0112 

3  2  *  ‘  0  *  *  0  I  *  2  2  -  3113 

33  *  3  1  1  0  2  -  *  3  r  ^  4  0  3 

fit  *71irp-'-i-^43'i>' 

35  *2000030  o'  0.  1000 

3  6  *  3  1  1  0  2  3  0  O  0  t  4  0  0  0 

*  7  *  3  1  1  0  2  '  0  *  •.  4  0  0  * 

*8  «  2  *  2  0  "  -  j  "  *  1  3  0  * 

3  p  *11102  3  CO*  -  0  *  *  * 

40  *111023000.0000 


*>  f\i  r*j  i>  f\ j r,j  -*>  &  ^ Oj  t>i  j  j  j  ^  uo.j  j  jcj 
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Appendix  L. 


Sequential  Simulation  Program  Listing 


PROGRAM 


( INPUT .CUTPui  )  5 


(•PURPOSE;  TO  SChEOULE  mu. T IPROOUCTS  IN  AN  UPS  TEST  FACILITY  WITH  RE6ARC 
(•TO  LIMITED  AVAILABILITY  .F  RESOURCES. 

prigpammer:  Charles  j.  .enable 

input:  INPjT  IS  DIVIDED  INTO  the  FOlLOwI N &  FOUR  sections: 

1.  TEST  STATION  EFFICIENCY  DATA 

2.  SERIES  VERSUS  TESTOR  Q UAL  I F I C A T I ONS 

3.  PRODUCTS  SCHEDULED  TO  BE  TESTED 
A.  PARRALLEL  SYSTEM  COMPONENT  ARRAY 

THE  FIRST  OF  THESE,  T. S T  STATION  EFFICIENCY  DATA,  REPRESENTS  THE 
COMPATIBILITY  AND  PREFER. NCE  R  A T I N  G  FOR  THE  TEST  STATION  AND  THE 
SPECIFIC  PRODUCT  TYPE.  *  2ER0  REPRESENTS  A  NONCOMPATIBLE  SITUATION  AND 
THE  NUMBERS  1  TO  10  ARE  .SED  TO  INDICATE  PR EF ER E NC E /E F F I C I ENC Y  FOR  THE 
TEST  S T  A  T 1 0  N  AND  THE  SPE-TFIC  PRODUCT  TYPE.  THE  NUMBER  1  IS  ThE  MOST 
PREFERED/HIGHEST  EFFICIENCY  AND  THE  NUMBER  10  IS  THE  LEAST  PREFERED/ 
LOWEST  EFFICIENCY.  the  format  of  this  data  is  AS  follows: 

ADMINISTRATIV.  LINE 
Ell  E 1 2  .  .  .  E1B 
£21  E2B 

•  • 

•  • 

EU1  EU2  •  •  .  EU8 

WHE»E  EIU  IS  the  EFFICIENCY  RATING  FOR  PRODUCT  I  IN  TEST  STATION  U, 

P  IS  ThE  NUMBER  OF  THE  l.ST  STATIONS,  AND  U  IS  TNE  NUMBER  OF  PRCCLCT 
TYPES. 

THE  SECOND  IS  SERIES  .ERSUS  TESTOR  0 U AL IF  I C A T I  ON S .  THIS  IS  A  0,1 
MATRIX  WITH  0  INDICATING  THAT  THE  TESTOR  IS  UNQUALIFIED  TO  TEST  THE 
SPECIFIC  PROCUCT  LINE  ANl  1  INDICATING  THAT  THE  TESTOR  IS  QUALIFIED. 

ThE  DATA  FOlLOwS  THAT  OF  EFFICIENCY  IN  THE  FOLLOWING  FORMAT: 

ADMINISTRATIV.  LINE 
Q 2  1  022  .  .  .  02W 

Q31  Q3W 


Q51  052 


05  W 


AND  TESTOR  U  AND 


WHERE  QIu  IS  THE  QUALIFI.ATION  RATING  FOR  SERIES  1 
W  IS  THE  TOTAL  NUMBER  OF  TESTORS. 

THE  THIRD  SECTION,  PRlDUCTS  TO  BE  TESTED,  CONSISTS  OF  AN  ARRAY  OF 
ALL  PERTINENT  data  FOR  L«CH  PRODUCT  SCHEDULED  TO  BE  TESTED  IN  THE 
PLANNING  HORIZON.  IT  CO..SISTS  OF  13  COLUMNS  ORGANIZEC  IN  THE  FOLLOWING 

FQOWAT ; 


COLUMN 

1. 


2. 

3. 

A. 
5. 
6  • 
7. 

B. 


THIS  IS  A  SE0UENT1AL  NUMBER 
NUMBER  OF  PRODUCTS  TO  PE 


IS 


10. 


11. 

1  2. 


13. 


DESCRIPTION 

ADMINISTRATION  NUMBER. 

FROM  1  TC  N,  WHERE  N  = 

TESTED  IN  ThE  H0RI70N. 

SYSTEM  NUMBER 
SER IES 

INPUT  .OUTAGE 
KW 

OUTPUT  VOLTAGE 
HERTZ  RATING 
PROOUC.  TYPE  NUMBER 

PARALLlL  SYSTEM  INDICATOR.  IF  r  0  THEN  PRODUCT 
NOT  PA*T  OF  A  PARALLEL  SYSTEM.  0THFRW1SE,  T HE 
ADMINI.TATION  NJMBER  OF  The  PARAL.EL  UNIT  IS  ENTERED. 
FOR  EXAMPLE,  IF  THE  PRODUCT  IS  A  PART  OF  PARALLEL 
SYSTEM  2  A  2  IS  ENTERED. 

COMPON. NT  NUMBER.  IF  =  0  THEN  NOT  A  PARALLEL  SYSTEM 
COMPON.NT.  OTHERWISE,  THE  COMPONENT  NJMBER  IS 
ENTEREs.  FOR  EXAMPLE,  IF  THE  PRODUCT  IS  THE  SECOND 
COMPON.NT  OF  A  PARALLEL  SYSTEM  THEN  2  IS  ENTERED. 
NUMBER  OF  REQUIRED  TEST  HOURS 

CURREN.  SCHEDULED  ARRIVAL  DATE.  THIS  DATE  IS  GIVEN 
IN  NUM.ER  OF  WORM  DATS  FROM  START  OF  PLANNING 
HORIZO..  . 

CURRENi  DUE  DATE.  THIS  DATE  IS  GIVEN  IN  NUMBER  OF 
WORK  O.YS  FRO**  START  OF  PLANNING  HORIZON. 


THIS  MATRIX  IS  IN  THE  F0-L0WIN6  FORMAT 


<  * 

ADMINISTATI,  E  LINE 

•  ) 

(  * 

PI,  1 

PI, 2  .  .  .  PI, 13 

*  ) 

(  * 

P2 , 1 

P2 , 1  3 

•  > 

(  • 

• 

♦  ) 

(  * 

. 

• 

•  ) 

(  * 

• 

• 

•  ) 

(  * 

PN,  1 

PN , 2  .  .  .  PN , 1 3 

•  ) 

WHERE  PIJ  IS  DATA  OF  PRO-UCT  I  IN  COLUMN  J,  AND  N  IS  THE  NUMB  ER  OF 
PRODJCTS  TO  BE  TESTED  IN  ThE  PLANNING  HORIZON.  PARALLEL  PRODUCTS  MUST 
BE  FIRST  in  THE  SEOUENCE  WITH  ARRIVAL  TIME  E C U AL  TO  ZERO.  ALL  OTHER 
PRODUCTS  MUST  BE  CRCEREu  EY  ARRIVAL  DATES,  EARLIEST  FIRST. 

THE  LAST  SECTION  OF  I„PUT  CONCERNS  THE  PARALLEL  SYSTEMS  SCHEDULED 
TO  BE  TESTED.  THIS  IS  A  0,1  MATRIX.  IF  r  C  THEN  THIS  COMPONENT  IS 
PRESENT.  IF  -  1  THEN  THE  COMPONENT  IS  NOT  PRESENT.  ThE  M A  X I M UM  • 

NuMpFR  OF  COMPONENTS  FOR  ONE  SYSTEM  IS  FOUR.  THE  COMFONENT  ARRAY  IS 
IN  ThE  FClLOwINF  FORMAT: 

ADMINISTRATIVE  line 
CU  C12  .  .  .  CIA 
C21  C2“ 


CNFSYS1 


.  CNPSVSA 


1  1 


LATE  OR  EXCEED  THE  HORIZON.  THIS 
IN  THE  FOLLOWING  7  COLUMNS  OF 


W u  E D  E  CIU  IS  THE  0,1  VAR. ABLE  REPRESENTATION  Of  COMPONENT  J  OF 
PARALLEL  SYSTE-  I,  ANC  Nr S y S  IS  THE  NUMBER  OF  PARALLEL  SySTE“S 
SCHEDULED  TO  EE  TESTED  I..  ThE  HORIZON.  FOP  EXAMPLE.  THE  LINE:  C  0 
INDICATES  THAT  yhe  PARALLEL  SYSTEM  IS  COMPOSED  OF  TWO  COMPONENTS. 

outf-ut:  THE  FIRST  THREE  SECTIONS  OF  THE  INPUT  ARE  ECHOED  FOLLOWED  BY 

A  LIST  OF  THE  PRODUCTS  W.hICh  WILL  BE  "  ‘ 

IS  FOLLOWED  BY  AN  ASSIGNMENT  SUMMARY 
DATA  FOR  EACH  PRODUCT: 

COLUMN  DESCRIPTION 

1.  ADMINISTRATIVE  NUMBER 

2.  SYSTEM  NUM3ER 

3.  TE. T  STATION  ASSIGNMENT 

A.  FImST  SHIFT  TESTOR  ASSIGNMENT 

5.  SEvONC  SHIFT  TESTOR  ASSIGNMENT 

6.  TE-TING  START  TIME 

7.  TESTING  COMPLETION  TIME 

IT  SHOULD  PE  NOTED  THAT  .F  ALL  ASSIGNMENT  NUMBERS  ARE  ZERO  THAT  THE 
PRODUCT  NEVER  BEGAN  TEST. NG  IN  THE  HORIZON  AND  WILL  NCT  BE  REFLECTED 
IN  The  STATISTICS  FOLuOw. NG  THE  ASSIGNMENT  SUMMARY. 

the  following  statist. cs  are  now  outputed: 

1.  NUMBER  of  late  products 

2.  ME. N  LATENESS 

3.  PROPORTION  OF  JOBS  TARDY 
A.  MEmN  JOB  FLOW  TIME 

THIS  IS  FQ..OWED  BY  A  STATEMENT  SPECKING  WHICH  DISPATCH  RULE  WAS 


USED. 


IF  DESIRED  A  COMP. ET  E  SCHEDULE  FOLLOWS  BT  DAT  AND  HOUR  OF  TEST 


STATION  AND  TESTOR  ASSIGNMENTS. 

ASSUMPTIONS:  ALL  PRODUCIS  ARE  ASSUMED  TO  ARRIVE  AT  1200HRS  ON  THE 

ARRIVAL  DAY.  ALL  PRODUCiS  ARE  ASSUMED  DJE  AT  1200HRS  ON  THE  DUE  DATE. 
TEST  STATIONS  ARE  ASSUME.,  EMPTY  AT  START  OF  HORIZON.  THIS  CAN  BE 
CHANGED,  HOWEVER,  BY  INSERTING  ASSIGNMENT  STATEMENTS  AFTER  INITIAL¬ 
IZATION  OF  ASSIGNMENT  MAiRICES.  PRODUCTS  ARE  ASSIGNED  TO  SECOND  SHIFT 
ONLY  IF  They  WILL  BE  LAT.  IF  PROCESSED  COMPLETELY  ON  FIRST  SHIFT. 

FIRST  SHIFT  TESTOR  ASSIGNMENT  OCCURS  IN  CONJUNCTION  WITH  TEST  STATION 
ASSIGNMENT  ANC  PRIOR  TO  -ECONO  SHIFT  TESTOR  ASSIGNMENT. 

LIMITATIONS:  the  ARRIVA.  AND  DUE  DATES  MUST  BE  IMPUTED  as  The  number 

OF  WORK  DAYS  UNTIL  THE  AlTUAL  DATE.  THEREFORE,  A  DUE  DATE  "UST  BE 
CHANGED  IF,  FCR  INSTANCE,  AN  ADCITICNAL  DAY  IS  ADDED  10  THE  WORK  WEEK 
PRIOR  TO  THE  ACTUAL  DUE  -ATE.  PRODUCTS  WHICH  ARE  NOT  ASSIGNED  OR 
WORKED  ON  DJRINS  ThE  PLA.NING  HORIZON  ARE  NOT  A  PART  CF  THE  TERMINAL 
STATI STICS. 

ERROR  CHECKING  AND  RESPO„SES:  IF  A  PRODUCT  CAN  NOT  COMPLETE  TESTING 

in  the  p.anning  horizon  «  statement  to  that  effect  is  printed  along 
with  the  message  that  it  will  be  late. 


124 


ALGORITHM/STRATEGY: 

00  FOR  Tzl  TO  T 

IF  PRODUCT  ARi\  I V  AL  T IME  zT  THEN 

NO V E  RRODUwT  TO  ASSIGNMENT  CllEUE 
IF  FIRST  S  H I F I  HOUR  THEN 

00  FOR  C A0 1  PROOJCT  IN  ASSIGNMENT  QUEUE 
IF  COMP»TIBLE  BAY  AVAILABLE  THEN 

if  qualified  testor  available  then 

ASSIGN  BAY  WITH  HA*  EFFICIENCY 
ASSIGN  FIRST  SHIFT  TESTOR 
IF  SECO.O  SHIFT  REQUIRED 

ASSIGN  FIRST  SHIFT  TESTOR  FOR  REMAINDER 
Of  SHIFT 

DECR.HENT  PROCESSING  TIME 
HOVE  PRODUCT  TO  SECOND  SHIFT  TESTOR 
Oj  CUE 

ELSE 

ASSIuN  FIRST  SHIFT  TESTOR  UNTIL  TESTING 
C-hplETE 

ASSIGN  BAY  UNTIL  TESTING  COMPETE 
FREE  BAY  AND  TESTOR  FOR  ASSIGNMENT  AT 
CwMpLETION  TIME  ♦  1 

ELSE 

00  FOR  E A C i  PROOJCT  IN  SECOND  SHIFT  TESTOR 
QUEUE 

IF  QUALIFIED  testor  available  then 
ASSIuN  SECOND  SHIFT  TESTOR 
ASSIuN  FIRST  AND  SECOND  SHIFT  TESTORS 
U.TIL  TESTING  COMPLETE 
ASSIGN  BAY  UNTIL'  TESTING  COMPLETE 
FREE  BAY  AND  TESTORS  FOR  ASSIGNMENT  AT 

completion  time  ♦  1 

ELSE 

ASSIuN  bay  for  t 

IF  L~ST  HOUR  OF  SECOND  SHIFT  THEN 

AuSIGN  BAY  AND  FIRST  SHIFT  TESTOR  FOR 
UPCOMING  FIRST  SHIFT 
Is  TESTING  COMPLETE  THEN 

FREE  BAY  AND  TESTORS  AT  COMPLETION 
TIME  .  1 

Ei.  SE 

DECREMENT  PROCESSING  TIME 

*  C-tS  200000.X  O.P  0  •> 


CONST 

N 

_ 

iqt; 

(‘NUMBER  OF  MACHINES  TO  BE  TESTED 

* 

W1 

= 

u; 

(•NUMBER  „F  TESTORS  ON  FIRST  SHIFT 

•  ) 

w2 

= 

7; 

(•NUMBER  j F  TESTORS  ON  SECOND  SHIFT 

•  ) 

B 

= 

lR ; 

(‘NUMBER  u  F  TEST  STATIONS 

«  ) 

U 

= 

72; 

(•NUMBER  w  f  product  types 

*  ) 

W 

z 

2  0 ; 

(•NUMBER  »F  TOTAL  TESTORS 

•  ) 

HOR 

r 

5; 

(•PLANNINu  HORIZON  IN  WEEKS 

•  ) 

OAY 

3 

16; 

(•NUMBER  uF  HOURS  IN  WORK  DAY 

•  ) 

W  M 

z 

5; 

(•NUMBER  uF  DAYS  IN  WORK  WEEK 

*  ) 

TIME 

Z 

400 ; 

(•NUMBER  uF  HOURS  IN  PLANNING  HOR I ZON zHOR *D A Y* WK 

«  ) 

SHIFT 

Z 

b; 

(•NUMBER  j  F  HOURS  IN  SHIFT 

•  ) 

NPSYS 

r 

4 ; 

(•NUMBER  v  F  PARALLEL  SYSTEMS 

•  ) 

PRINT 

- 

0; 

(•PRINTED  SCHEDULE  SPECIFIER:  IF  z  1  PRINT 

*  ) 

TYPE 

SPOINTER  z  3SYSTEMNO0E!  (.QUEUE  POINTER 
PR  CD  z  RECORD 

SERIES  :  INTEGER;  (‘SERIES 

voltin  :  INTEGER;  (.INPUT  voltage 

kw  :  integer;  (.power  size 

VOLTOUTI  INTEGER!  (.OUTPUT  VOLTAG 

hertz  :  integer;  (.hertz 

proonum:  integer  -  - - ‘ 

end;  (.prod*) 


(  ‘SERIES 
(.INPUT  VOLTAGE 
(•POWER  SIZE 
(•OUTPUT  voltage 
( ‘HERTZ 

l ‘ADMIN  NUMBER  ASSIGNED  TO  PRODUCT  TYPE 


SYSTEM  :  RECORD 

ADMIN  :  INTEGER!  k  *  S£  OUE  NCE  NUMBER  * 

STSNJM  :  INTEGER!  k  ‘SYSTEM  NUMBER  • 

sys  :  prod:  i  *ups  description  • 

PL  :  INTEGER!  k  *  NUMBER  OF  PARALLEL  SYSTEM  TC  WHICH  THIS  • 

( ‘COMPONENT  BELONGS  * 

COMP  :  INTEGER!  i ‘NUMBER  OP  COMPONENT  IN  PARALLEL  SYSTEM  • 

TESThR  :  INTEGER!  k ‘NUMBER  OF  HOURS  REQUIRED  TO  TEST  UPS  • 

cscbs  :  integer:  ‘‘Current  arrival  date  to  test  * 

OUEOT  :  INTEGER:  ( ‘CURRENT  PROMISE  date  • 

START  :  INTEGER;  i‘TIM£  PRODUCT  PIRST  ASSIGNED  TP  TEST  STATION  * 

FINISH  :  INTEGER;  i • T I  ME  TESTING  COMPLETED  ON  PRODUCT  • 

FLAGP  :  INTEGER:  k ‘PARALLEL  SYSTEM  ALREADY  IN  BAT  OUEuE  IF  =  I  • 

FLAG’  :  INTEGER:  k‘P»OOuCT  ALREADT  IN  SECOND  SHIFT  QUEUE  IF  a  1  • 

bat  :  integer;  ^number  of  assigned  bay  •> 

T 1  :  INTEGER;  k ‘NUMBER  of  assigned  first  shift  testor  *> 

Y2  ;  integer  (-number  of  assigned  second  shift  testor  *> 

ENC;  {‘SYSTEM.) 

STSTE"N00E  =  RECORD 
stsinfo:  system; 

LINK  :  qpotnter 

end;  ( -systemnode* ) 


JPS  =  ARRAYS1..N3  OF  ST^TEM!  (‘ARRAY  OF  PRODUCTS  TO  BE  TESTED  •) 

V  AP 

0  l^EAD  |Q2H£A0  .NEYT.P  ,0  ,(k  .S  :  QP01NTER! 

Eii«j«KtL«M«x  :  integer* 
product  :  ups; 

EFF  ;  ARRAY81..U.1.. B9  OF  INTEGERS  (‘BAY  EFFICIENCY  *) 

QUAL  :  ARR ArS2. .N , 1 ..  va  OF  INTEGER;  (‘TESTOR  QUALIFICATIONS  •> 

3ASSI6N  :  ARRAYS1..9 .0..  times  of  INTEGER!  (‘BAY  ASSIGNMENT  •) 

tassign  :  ARR&Yai..w,o.. times  of  integer;  ctestor  assignment  •> 

T  :  INTEGER;  (‘Tint  •) 

flag  :  integer;  (.ch.ck  flag:  if  =  1  violation  •) 

TEST0F1  :  INTEGER!  (‘FUST  SHIFT  TESTOR  NUMBER  *) 

TEST0R2  :  INTEGER!  (‘SE.OND  SHIFT  TESTOR  NUMBER  •> 

MAX  ;  INTEGER!  (‘MAaIMUM  NUMBER  IN  SEARCH  ‘> 

TEST1.TEST2  :  INTEGER!  k ‘TESTOR  ASSIGNMENT  FLAG;  IF  :  1  TESTOR  ASSIGNED  ») 
PCOMP  ;  ARRAT31..NPSYjtl..A8  OF  INTEGER!  (‘PARALLEL  SYS  COMPONENT  array* 

LATE  :  INTEGER!  (‘NUiBER  OF  LATE  PRODUCTS  •) 

LATENESS;  INTEGER:  (‘TOiAL  number  of  LATE  PERIODS  ») 

flCu  :  integer;  (.amount  of  time  in  system  •> 


PROCEDURE  OCREATE  (VAR  QME a  0 : 0 p 0  IN T ER  )  ; 

(‘PURPOSE:  CREATES  A  QUEU_  WITH  QHEAD  POINT  ING  TO  The  FIRST  (DUMMY)  ELEMENT.) 


BEGIN 

NEw(QHEAO)! 

WITH  GHEADS,  SYSIVPO  DO 
BEGIN 

admin  :  =  o; 
stsnum  :=  o; 

STS. SERIES  ;r  0! 

SYS. VOLTIN  :=  o; 

sys.kw  :-  o; 
sys. volt  out  :=  o; 
SYS. HERTZ  0! 
sys.proonum  :=  o; 
pl  o; 
comp  :=  o; 
testhp  :  =  o; 
cschs  ;=  o; 
duedt  : r  o; 

START  :=  o; 

FINISH  ;z  0; 

flagp  :=  o; 

FLAG2  l-  05 
BAY  :=  05 

t i  :=  o; 
t 2  :=  o; 
link  :=  NIL 
END 

end;  (»qcpeate*) 


PPOCEDURE  CREATENODE  (IrlNIFGER!  VAR  NEXT : Q PO IN  TER > ! 

(♦purpose:  allocates  a  ne.  nooe  for  storage  of  the  next  product  in  the  bueueo 


BEGIN 

new(ne»T); 

NEX  T  3. $Y  S  INFO  :=  PRODUCTS!*? 
NEXT3.LINK  :=  NIL 

end;  (.createnooe*) 
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PROCEDURE  CHECK  <J,T:INTEG.RS  GIOPOINTERS  VAR  FI  AG : I NTE *ER ) 5 

(•»u*pose:  checks  to  see  .f  sat  constraints  are  violated  bt  possible  assign-  •  > 

(•KENT  •) 

VAR 

k, i  :  integers 


BEGIN 

IF  (J  =  A)  OR  (J  =  5>  then 
BEGIN 

(*S0h2  UNITS  NOT  ALLOVEj  IN  B1  OR  B2*) 

IF  G3. SYSINFO. SYS. HERTZ  =  50  THEN 
FLAG  :=  1! 

(•IF  25  OK  V  UNIT  IN  B1  On  B2,  HA*  SIZE  IN  OTHER  IS  60KV»> 

IF  G3.SYSINFO.SrS.KV  =  ..5  0  THEN 
BEGIN 

IF  8ASSIGN3A.T3  O  0  THEN 

IF  PRODUCTSBASS  IG<«  3  5  t  T S3.  SYS.KU  >  60  THEN 
FLAG  :=  IS 

IF  BASSIGNS5*T3  <>  0  THEN 

IF  PR0DUCT3BASSIGaS5iT38.SYS.Ktf  >  60  THEN 

flag  :=  1 

ENDS 

IF  G3.  SYSINFO.  SYS.KU  =  o0  THEN 
BEGIN 

IF  BASSIGN3A.T3  O  0  THEN 

IF  PROD UCT3BASSlGia8A«T33. SYS.KU  >  250  THEN 
PLAG  :=  IS 

IF  BASSIGN35.T8  O  0  THEN 

IF  PR0DUCT3BASSIG.35 «T33. SYS.KU  >  250  THEN 
FLAG  S=  1 

END 

ENDS 

IF  ( J  =  6 )  OR  ( J  =  7 )  THEN 
BEGIN 

(*C*D<=330KV*> 

IF  BASSIGN36.T3  <>  8  TH_N 

IF  (PRD0UCT3BASSIGN3jtT33.SYS.KU  ♦  3 3 . SY S I NF 0 . S Y S .K U )  >  330  THEN 
FLAG  :=  IS 

IF  BASSIGN37 1 T3  <>  0  T(C  N 

IF  ( PR 00U C T3B ASS  I GN3(  .T33.SYS.KV  ♦  G3  .  SYS  I  NF 0 . S VS  .K u )  >  330  THEN 

flag  :=  1 


IF  (  J  =  8  )  OR  (J=9)  THEN 
BEGIN 

IF  BASSIGN3S.  T3  <>  0  TH.N 
BEGIN 

(*f.F<=600ku*> 

IF  (PR00UCT3BASSIGN3O .T33. SYS.KU  ♦  G3 . SYS  I NF 0. S YS .K U >  >  600  THEN 
FLAG  :=  15 
(*E*F<=1(50HZ)*) 

IF  PR0DUCT3BASSIGN38. T33. SYS. HERTZ  =  50  THEN 
IF  G3. SYSINFO. SYS. HERTZ  =  50  THEN 
FLAG  :=  IS 
(*E*F<=1(2090UT)*) 

IF  PR0DUCT3BASSIGN38. T33.STS .VOLTOUT  :  208  THEN 
IF  G3.STSINF0. STS. VOLTOUT  =  208  THEN 
FLAG  S-  1 

END 

ENOS 

IF  (JrlO)  OR  ( J  =  1 1 >  THEN 
BEGIN 

IF  BASSIGN310.T3  <>  0  TnEN 
BEGIN 

(»G1*G2<=250KVO 

IF  (PR00UCT3BASSIGN3. 0,T3i. SYS.KU  ♦  G3. SYSINFO. SYS .KU >  >  250  THEN 

flag  :=  is 
( *G1»G?<  =  I ( 50hZ> •  I 

IF  PR0DUCT3BASSIGN31. ,T33. SYS. HERTZ  =  50  THEN 
IF  G3.STSINF0. STS. HERTZ  =  50  THEN 

flag  :=  l 

ENOS 

IF  BASSIGN3U.T3  <>  0  TnEN 
BEGIN 

IF  (PR0DUCT3BASSIGN3. 1,T3S. SYS.KU  ♦  G3. STS  INFO .SYS  .KU  )  >  250  THEN 
flag  :=  is 

IF  PR0DUCT33ASSIGN31. .T33. STS. HERTZ  s  50  THEN 
IF  G3. SYSINFO. SYS. HERTZ  s  50  THEN 
FLAG  1 

END 

END 

ENOS  ( 'CHECK • ) 
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PROCEDURE  QINSERT  ( NEXT , QH_ AD IQPOINTER )  ! 

(•PURPOSE:  INSERT  PRODUCT  IN  OUEUE  WI TH  DESIRED  DISPATCH  RULE.  CODE  •) 

(•REQUIRES  ALTERATION  IF  Oi HER  THAN  PRESENT  DISPATCH  RULE.  EARLIEST  DUE*) 
(•DATE  WITH  TIES  BROKEN  BY  GREATEST  Kb.  IS  DESIRED.  •) 

(•INSERT  PRODUCT  IN  LIST  W.TH  EARLIEST  DUE  DATE  FIRST*) 

VAR 

a.b.c  :opointer; 

I.K.L.H  :  INTEGER! 

PROCEDURE  INSERTNODE  ( NEX T.  G P OINTE R  i  VAR  D:Q POINTER!  EIOPOINTER)! 

BEGIN 

OS.LINK  :=  NEXT! 

NEXT3.LINK  :=  E 
END!  (* INSERTNODE* ) 


BEGIN 

(•SEARCH  FOR  CORRECT  POSITION*) 

B  :=  OHEAD! 

while  ( NEXTS.SYSIVFO.DUEDT  >  B 8 . S Y SI NF 0 . DUE D T )  AND  (38.LINX  <>  NIL)  DO 
BEGIN 

A  !=  B8.LINK! 

3  CS.LINK! 

A  BS.LINK 

END!  (*UHILE*) 

(•INSERT  PRODUCT  IN  PROPER  POSITION*) 

IF  NEXT8.SYSINFO.OUEDTrB8.-YSINFO.DUEDT  TtN 
BEGIN 

IF  NEXT8.SYSINFO.SYS.KW  >  B8. S Y SI NFO . S YS .K W  THEN 
(•INSERT  BEFORE  B*> 

INSERTNODE(NEXT.C.B) 

ELSE  IF  NEXT8.SYSINF0.Si S.KU  =  B8 . S YS I NF 0 .$ Y$ .K U  THEN 
BEGIN 

I  : s  NExTS.SYSINFC.DoECT! 
d  :=  Ba.SYSINFO.DLEDi ! 

L  :=  NExT8.SYSINF0.Si  s.kw; 

n  :=  ea.SYSiNFo.sYs.vw; 

WHILE  (L=H>  AND  (I=J i  AND  (B8.LINK  <>  NIL)  DO 
BEGIN 

A  :=  B8.L INK  5 
C  :=  Bi 
B  :=  C8.LINK! 
a  :=  Ba.LlNx; 

J  :r  Ba.SYSlNFC.D^ECT ; 

M  Ba.SYSINPO. Sr  s.x  W 
END!  (‘WHILE*) 

IF  (IOJ)  OR  (LOP)  j  HEN 

(•INSERT  NODE  BEFORE  B«) 

INSERTNOCE(  NEXT.  C.  B) 

ELSE 

(•ATTACH  AT  END  Or  QUEUE*) 

INSERTNODC(NEXT,B, NIL  ) 

ENO  (‘ELSE  IP*) 

ELSE 

(•INSERT  AFTER  B*> 

INSERTNODE (NEXT. B. A) 

ENC 

ELSE  IF  NEXTS.SYSINFO.DUEDi  <  B 8 . S Y S  I  NF 0  .DUEDT  THEN 
(•INSERT  NODE  BETWEEN  C  AND  B  * ) 

INSERTNODE(NEXT.C.B) 

ELSE 

(•ATTACH  TO  END  OF  QL'EU-*) 

INSERTNODE (NE  XT .3 .NIL ) 

END!  (»QINSERT*) 
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BEGIN  (*EXI0E*> 

(•READ  BAY  EFFICIENCY  DATA*) 

WRITELNS 

WRITELN  (»PR0DUCT  VERSUS  T.ST  STATION  EFF 1 C I ENC Y : • ) { 

WRITELNS 

WRITELN  ( ’ TEST  STATION* :35<  • 

write  (• 

FOR  I  :=  1  TO  B  DO 
write  <i:3>; 

WRITELNS 

WRITELN  ( ’PRODUCT””””” > ; 

RE40LNS 

FOR  I  :=  1  TO  U  00 
BEGIN 

WRITE  (IMt*  *•  )S 
FOR  J  :=  1  TO  e  DO 
BEGIN 

READ  UNPUT. EFFil.ja*  S 
WRITE  (EFFal.oa:3) 

END  S 
reaoln; 

WRITELN 
ENDS 

reaoln; 

(•READ  TESTOR  OUAUFICATIO.S’) 

WRITELNS 

WRITELN  (’SERIES  VERSUS  TEo  T0e  OUAL I F I C A TI ON S : • ) 5 
WRITELNS 

WRITELN  (’TE$TOR’:AO)S 

write  (•  •  •>; 

FOR  I  :=  I  TO  W  DO 
WRITE  (1:3)5 
WRITELNS 

WRITE  (  ’SERIES* ) S 
FOR  I  :=  1  TO  67  DO 
WRITE  (•”); 

WRITELNS 

FOR  I  S=  2  TO  5  DO 
BEGIN 

WRITE  (ISA.*  •*>; 

FOR  J  : -  1  TO  W  DO 
BEGIN 

READ  (  INPUT, QUAL8I.J. )  S 
WRITE  (QUALSI, JS:3) 

ENOS 
REA0LNS 
WRITELN 
ENDS 
READLNS 

(’INITIALIZE  ASSIGNMENT  MAiRICES  TO  EMPTY’) 

(•INSERT  READ  STATEMENTS  M.RE  IF  STATIONS  ARE  NOT  EMPTY  AT  TIME  ZERO’) 

FOR  I  :=  1  TO  B  DO 

FOR  U  lr  0  TO  TIME  DO 
BASsiGNai.ja  :=  o; 

FOR  I  *.=  I  TO  W  DO 

FOR  J  :z  0  TO  TIME  DO 
TAssiGNai,ja  :=  o; 

(•READ  INPUT  DATA*) 

WRITELNS 

WRITELN  (’THE  FOLLOWING  AR_  THE  PRODUCTS  SCHEDULED  TO  BE  TESTED:*); 

WR I  TEL  N  5 

WRITE  (•  ADMIN  STSNUM  SLUES  VOLTIN  K W  VOLTOUT  HERTZ  PROONUM’) 

WRITELN  (♦  PL  COMP  TESTHR  CSCHS  DUEDT’IS 


1  TO  N  DO 


FOR  I  := 

BEGIN 

WITH  PR0DUCT8I8  DO 
BEGIN 

READ  ( INPUT* ADMIN,ST.,NUH,SYS.SERIES,SYS.V0LTIN,SYS.KW ) ; 

READ  (INPUT, SYS. VOLT. UT, SYS. HERTZ *SYS.PRODNUM, PL *COMP,TESTHR,CSCHS) ; 

READ  (INPUT, DUEDTIS 

READLN*. 

start  :=  o; 

finish  :=  o; 

flagp  o; 

FLAG2  :=  0! 

bay  :=  o; 
ti  :=  o; 

T2  :=  0! 

write  <admin:a,sysnuh :g,sys. series: 7, sys.voltin:9,sys.kw: 8) ; 

WRITE  (SYS. VOLT OUT:b,  SYS. hertz : B ,SYS.PRODNUH :8 ,PL :*,:OHP: 8.TE  SThR:7) 
WRITELN  (CSCHS:9,DUE-T:8> 

END 

end; 

readln; 

(.READ  PARALLEL  SYSTEM  COMPONENT  ARRAY*) 

FOR  I  :=  1  TO  NPSYS  00 
BEGIN 

FOR  J  :=  1  TO  R  DO 

READ  (INPUT, PC0MP81**8)« 

READLN 

END; 

writeln; 

(•INITIALIZE  QUEUES* ) 

OCREATE( Q1HEAD)  ; 

8CREATE(02hfad> ; 

(•INITIALIZE  ACCUMULATORS*, 

late  :=  o; 
lateness  :=  o; 
flow  :=  o; 

(•SINCE  TEST  BAYS  ASSUMED  -MPTY  AT  START, TIME  INITIALIZED  AT  1200  HOURS*) 
(•FIRST  DAY*) 

t  :=  J? 

E  :  =  npsys  *  1; 

(•THE  FOLLOWING  FOR  LOOP  I*  USED  IF  ALL  PROOUCTS  ARRIVE  ON  THE  SAME  DAY.) 
(•AND  SHOULC  BE  USED  IN  PL.CE  OF  THE  PREVIOUS  ASSIGNMENT  STATEMENT:  •> 

(•  •  ) 

(•  FOR  I  :r  (NPSYS**)  TO  N  DO  •> 

(•  BEGIN  •) 

(•  CREATENODE(I»,EXT)  •) 

(*  QINSERT(NEXT,. 1HEAO)  •) 

(•  END  •) 

(*  •> 

(•ADDITIONALLY,  THE  FOLLOW, NG  STATEMENTS  WHICH  APPEAR  AFTER  THE  BE6IN  •) 

(*DF  THE  TIME  LOOP  SHOULD  jE  DELETED:  •) 

(  *  •  > 

(•  FLAG  0  •> 

(♦  WHILE  (FLAG=0)  «ND  (I<=N)  do  *) 

(•  BEGIN  •) 

(•  IF  P( PRODUCT. 18. CSCHS-1)*DAY»3>  =  T  THEN  •) 

(•  BEGIN  •) 

(•  CREATENOO,.(I,NEXT)  •) 

(•  SINSERT(  N_  XT,  Q1HEAD  )  •) 

c  I  :=  i*i  •) 

(•  END  •) 

(•  ELSE  «) 

(•  r  LA  G  I =  1  •) 

(•  END  •) 

WHILE  T  <=  TIME  DO 

(•RUN  TIME  LOOP  UNTIL  E..0  OF  PLANNING  HORIZON*) 

BEGIN 

flag  :=  o; 

while  (flag=o  and  (E<=.)  do 

BEGIN 

IF  ( (PR0DUCT8E8,CSCH*-1 )»DAY*3)  r  T  THEN 
BEGIN 

CREATENODE(E,NEXT) ; 

0  INSERT! NEXT, 01HC»D) ; 

E  :=  E* l 
END 
ELSE 

flag  :=  l 

ENDS 
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(•DETERMINE  If  FIRST  SM.FT  HOUR:  IF  SO  RUN  ASSIGNMENT  QUEUE*) 

IF  (NOT  ODD  (  T  DIV  SHIFT/  >  OR  (T  OIV  SHIFT  =  0)  THEN 
BEGIN 

(•ESTABLISH  ASSIGNMENT  QUEUE  POINTERS*) 

o  :=  oihe«d; 

P  09. LINK  5 

(•MOVE  THROUGH  ASS IG . RENT  QUEUE*) 

WHILE  P  <>  NIL  00 
BEGIN 

max  :=  u ; 
testi  :=  o; 

TEST0R1  :=  i; 

FOR  U  ;=  1  TO  B  D- 
BEGIN 

flag  :=  o> 

(•IDENTIFY  BAY.  AVAILABLE  AND  DETERMINE  3 A Y  WITH  MAX  ErFICIENCY* 
IF  EFFSPi.SYSI..  F0.SYS.PRODNUM»ja  <>  0  THEN 
BEG  T  N 

IF  BASSIGNaw.TS  =  fl  THEN 
BEGIN 

check  <  ji i ,p,flag> ; 

IF  flag  =  0  THEN 
BEGIN 

IF  EFr  8P 8 • SY S INFO • S YS .PRODN UM . JS  <  MAX  THEN 
BE/  IN 

M  Aa  :=  EFF8P8.SYSINF0.S YS. PR ODNUM , J 8 5 

X  .  =  J 

ENw 

END 

END 

END 
END  ; 


(•IF  COMPATIBLE  B„ Y  FOUND  CHECK  FIRST  SHIFT  TESTOR  AVAILABILITY*) 
IF  MAX  <>  11  THEN 
BEGIN 

WHILE  (TESTlsOz  AND  ( T EST0R1 <  =  W 1 )  DO 
BEGIN 

(•CHECK  FIR/T  SHIFT  QUALIFICATIONS  FOR  PRODUCT*) 

IF  3UAL8P8./YSINF0. SYS. SERIES. TEST0R18  C>  0  THEN 
BEGIN 

(•IF  TEST  OR  AVAILABLE  ASSIGN  PRODUCT*) 

IF  TASSIxN8TEST0Rl ,T8  =  0  THEN 
BEGIN 

P2.SYwlNF0.Tl  :=  TEST  OR  1 ! 

PRODU.  T8PS.SYSINF0.  ADMINS.  T1  TESTORli 
TASSI- NSTESTOR 1 *Ti  l-  P 3 . S Y S I NF 0 . A C M IN i 

ps. SYw info. bay  :=  x; 

BASSIwN8PS.SYSINF0.BAY.TB  1  =  P8 . S Y S I NF 0 . AD M I N ; 

TESTI  :=  1$ 

PRCDu. T8PS.SYSINF0. ADMINS. BAY  ::  x; 

PRCDU. TSPS. SYSINFO. ADMINS. START  :=  T) 

PS. SY/INfO. START  :  =  T 
END 
END! 

TEST  OR  1  :=  iESTORl  ♦  1 

END  (‘WHILE  TC  S  T  1  •  ) 

End;  (*if  max*/ 


(•CHECK  IF  SECOND  SHIFT  REQUIRED*) 

(•CURRENT  PREREQU.SITE  IS  IF  TESTING  ON  FIRST  SHIFT  ONLY  WILL*) 
(•CAUSE  PRODUCT  Tw  BE  LATE.  IF  IT  IS  DESIRED  TO  EXPAND  THIS  •) 
(•CATEGORY  THE  DEwIREO  NUMBER  OF  HOURS  OF  THE  EXPANSION  •) 

(•SHOULD  BE  SUBTRACTED  FROM  THE  FIRST  ELEMENT  CF  THE  NEXT  IF  •) 
(•STATEMENT.  FOR  EXAMPLE.  IF  IT  IS  DESIRED  THAT  ALL  PRODUCTS*) 
(•THAT  ARE  SCHEDU. ED  TO  BE  COMPLETED  WITHIN  J  HOURS  Oc  THEIR  •) 
(•DUE  DATE  WITH  F.RST  SHIFT  TESTING  ONLY  BE  CONSIDERED  FOR  •) 
(•SECOND  shIFY  TEwTING  YHEN  thE  NUMBER  THREE  WCULD  BE  •) 

(•subtracted  from  the  first  element  of  the  next  if  statement.*) 

(♦PRODUCTS  ARE  AS.UMED  DUE  AT  1200HRS  ON  DJE  DATE*) 

L  :=  ( P8. SYSINFO. w UEDT-1 ) «SH1FT*J ; 


IF  (L-(T  DIV  DAY). SHIF T* < T  MOD  DAY))  <  ( P8. St S I NF 0 . TES 1  HR )  THEN 
BEGIN 

(•SECOND  SHIFT  REaUIAEO*  ) 

IF  TEST1  =  1  T-tFN 
BEGIN 

(•ASSIGN  FUST  SHIFT  testor  for  remainder  OF  SHIFT*) 

FOR  K  :s  T  iO  (  T*  SH IFT- 1  -  (  T  HOD  SHIFT))  DO 
BEGIN 

M  •  z  K  « 

If’m  >  T.ME  THEN 

m  : =  i inf ; 

TASSIGNJi  EST0R1-1,M8  :  =  PS  .S  Y  S  I  NF  0  .  AD  M  I N  ! 
BASSlGNSi-8.STSINF0.BATtR8  :=  P  8.  S  Y  S  I  NF  0  .  A  DM  I  N 

end; 

(•calculate  test  HOURS  REGAINING  TO  BE  ASSIGNED*) 
PS.SYSINFO.I  FSTHR  PS  .  SY-SI  NFO.  TE  S  THR  -  SH  I  F  T  ♦  (  T  MCC  SHIFT) 
»  :=  T ♦ SH I F i -<T  ROD  SHIFT); 

if‘m  o’time  THEN 

h  :  =  Tin.; 

BASS  IGNSPS.^  TSINF0.9A  TtNS  Ir  Pi  .  S  T  S  I  NFO.  ADM  I  N  ; 

(*«DVE  PRODUCT  TO  SECOND  SHIFT  CUEuE»> 
as. l ink  :=  r a. link ; 

G INSERT (  P  ,Q. HEAD  )  ; 
p  :=  qs.lINa 
END 
ELSE 

(•FIRST  SHIrT  TESTOR  NOT  AVAILABLE-PRODUCT  MUST  REMAIN*) 
BEGIN 
o  :=  p; 

P  PS.LlN\ 

ENC 

END  ( *  I F  SECON-  SHIFT  REOUIREC*) 


ELSE 

(•SECOND  SHIFT  NOT  RE3UIRED*) 

BEGIN 

IF  TEST1  r  1  TnEN 

(•ASSIGN  FUST  SHIFT  TESTOR*) 

BEGIN 

FOR  K  :=  T  iO  (T*SHIFT-1-(T  MOD  SHIFT))  DO 
BEGIN 
L  !  -  K  ; 

if'l  >  t* me  then 
l  ;=  i  ime ; 

TASSIGNS,  ESTORl-ltLS  :=  P  8  .S  TS  INF  0  •  »  C  M  I  N  *, 
BASSIGNat-S.SYSINFO.BAY.LS  P8.SYSINF0. ADMIN 

pnd  ; 

(•CALCULATE  TEST  HOJRS  REMAINING*) 

L  SHIFT-iT  mod  SHIFT); 

P8.SYSINF0.1 ESTHR  I:  P8.SYSINF0.TESTHR-L ( 

»  : r  P3.SYS*NFO.TESTHP; 

FOR  J  :=  0  iO  ((X  OIV  SHIFT)-!)  DO 
FOR  K  0  TO  7  DO 
BEGIN 

m: =  t  ♦. *8* (day • j )  *k  ; 

IF  M  >  TIME  Then 
M  .  r  TIME*. 


TASSIwN8TCST0Rl-l,M8  :  =  PS . SYS  I NFO . A  DM  I N 

end; 

m  :=  t«l*b*»x  div  shift).day; 

L  :=  M- 1 ♦ ( X  MOD  SHIFT); 

Pa.SYSINFO.r  INISH  :=  L ! 

PR0DUCT8F8.O YS INFO. ADMIN8. FINISH  :  =  L • 

IF  Pi.SYSlNr O.COMP  <>  0  THEN 

PC0MP8P*.  SYSINF0.PL»P8. SYSI NF O.C OMPi  :=  L*. 

IF  L  >  TIME  THEN 
BEGIN 

L  I  r  TIM.; 

NRITELN  (‘PRODUCT  * , P8 . S Y S  INF  0 . ADM  IN : 3 , •  EXCEEDS  HOR« ) 

end; 

FOR  K  ;r  M  I O  L  DO 

TASSIGNS! ESTOR1-1, *8  :  =  P8 .S Y S I NFO . ADM  I N ; 

FOR  K  :=  T  I  0  L  DO 

BASS IGNir 8. SYS  INFO. BAY, Ki  1=  P8  .SY  S  I  NF  0  .  A  0“  I N  ; 

L  :=  M- 1  * ( X  MOD  SHIFT )-( (P8. SYSINFO. DJEDT-1 > «DA Y ♦?) ; 
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IF  L  >  3  TH.N 
BEGIN 

WRITE  (  *r  R  0  DU  C  T  • . P 8 . S Y S I N FO . ADM  IN  I J , •  WILL  BE  • «  L : 3  >  * 

WRITELN  (•  TINE  PERIODS  LATEMS 
LATENESS  1=  LATENESS  »  LS 
LATE  :=  uATE  ♦  1 

pnd; 

08. L  INK  :=  I-  8. LINK  5 
DIS’OSEIP)  ! 
p  :=  oa. link, 

END 

ELSE 

(•FIRST  S  H I  -  T  TESTOR  NOT  A VA I L A  3 L E-PR  00 UC T  MUST  REMAIN  IN  OuEJE»> 

BEGIN 

G  i  =  PS 

P  :=  PS  .  L  I  Nr 

END 

END  (.else*) 

END  <  *  WHI L  £  P  O  ..IL*> 

fnd:  « » i f  not  odoo 

(•Check  Ic  SECOND  Shift  hojr-if  so  RUN  SECOND  SHIFT  QUEUE*) 

IF  00  0  (  T  DIV  SHIFT)  THE.. 

0  EG  IN 

(•initialize  oueue  pointers*) 
s  Q2hEao; 

=  : =  s  a. l ink  : 

(•MOVE  through  secon.  shift  QUEUE*) 

.mIlf  R  <>  NIL  00 
°EG  I  N 

TEST2  :=  05 
TESTJR2  : i  W1  *  1. 

WHILE  (TEST2  =  0)  A.,0  (TEST0»2<  =  W>  DO 
BEGIN 

(•CHECK  SECOND  SHIFT  0 J AL I F I ; A T I  ON S  FOR  PRODUCT*) 

IF  QUAL8R8. SYS. NFO. SYS. SERIES, TEST0R28  <>  0  THEN 
BEGIN 

(•IF  TESTOR  AVAILABLE  ASSIGN  PRODUCT*) 

IF  TA  SSIGN8  .  EST0R2*  T8  =  0  THEN 
BEGIN 

TASSIGNSi  EST0R2.T8  R8  .  S  YS  I  NFO  .  ADM  I  N  5 
8ASSIGN8r8.SYSINF0.BAY.T8  5=  R 8 . SY S I NF 0 . A  DM  I N 5 
T  EST  2  15 

RS.SYSIN--0.T2  TEST0R25 
PROOUCT8k8.SYSINFO.ADMIN8.T2  1=  TE  S  TOR  2 
END 

E  NO  5  (  *  I  F  Qo  A  L  •  ) 

TEST0R2  5=  TESi  0R2  ♦  1 
ENOS  (‘WHILE  T_  ST  2  * ) 

IF  TEST2  r  I  THEN 

(•BOTH  SHIFTS  r  SSI6NED • ) 

3E  G  I  N 

L  T*R8.SYSI,F0.TESTHR-15 
R8. SYS  INFO. F IN. SH  I r  L  5 

P0CCUCT8R8. SYS. NFO. ADMINS. FINISH  5=  L! 

IF  R8.  SYSI  NFO..  OMP  <>  0  THEN 

PC0"PSR8. SY. INFO. PL, R8. SYSINFO.  COMPS  5-  L5 
IF  L  >  TIME  TH.N 
BEGIN 

L  5=  TIMES 

URITELN  CPkODUCT  Y.Ra.SYSINFO. admin: 5*'  EXCEEDS  HOR»> 

ENDS 

FOR  K  : r  T  TO  .  00 
BEGIN 

IF  00D<  K  DI.  SHIFT)  THEN 

(•SECOND  SHIFT  HOUR-ASSIGN  SECOND  SHIFT  TESTOR*) 

BEGIN 

TASSIGN8I EST0R2- 1, K8  :=  R 8 . S Y S I NFO . A  0 P I N S 
BASSIGN8n  8  .S  YSINFO  .BA  Y,K8  :=  R8  .  S  Y  S  I  NF  0  .  A  DM  IN 
END 
ELSE 

(•FIRST  „HIFT  HOUR-ASSIGN  FIRST  SHIFT  TESTOR*) 

PFGIN 

T  A  SS I GN8n B.SYSINF0.T1 . K  8  R 8 .S YS I NF C . AD M I N 5 
RASSIGN8i8.SYSINF0.BAY,KS  :  =  R8.SYSINF0. AOMIN 
END 

ENDS  ( •FOR  r • ) 

L  :=  T*R8.SYSI..FO.TESTHR-l-((R8.SrSINFO.DUECT-l  )*DAY*3)  5 


i 

L 
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If  L  >  0  THEN 
BEGIN 

WRITE  t'PPO^UCT  '.Ra.SYSlNFO.ADHINlS. •  WILL  BE  »»L:J)5 
WRITELN  (•  i I ME  PERIODS  LATE'); 

LATENESS  :=  LATENESS  ♦  LI 

LATE  :=  LAT.  ♦  1 

ENDS 

(•PROOUCT  IS  ASSIGNED  ON  SECOND  SHIFT  AND  MOST  BE  REHOvED*) 
ss.link  :=  R3._ink; 

DISPOSE (R ) ; 
r  :=  ss.link 

ENO  t » I F  TE  S  T2*  ) 

ELSE 

<  .PRODUCT  MJST  REGAIN  IN  SECDNO  SHIFT  TESTOR  GJEUE') 

BEGIN 

BASS  IGN3R3. SYS. NF0.6AY , T8  :=  RS . STS  I Nf 0. ADH IN S 
BASSIGN3Ra.SVS.NF0. BAY, T*18  S=  P  8.  S  YS  I  NF  0  .  A  DH  IN  ; 

if  (T  moo  ship;  i  =  7  then 

(•ASSIGN  FIiST  Shift  TESTOR  FOR  UPCOMING  FIRST  SHIFT*) 

BEGIN 

IF  R3.SYSlNr  O.TESTHR  >  SHIFT  THEN 
BEGIN 

l  :  =  t  ♦  p  ; 

IF  L  >  T.M£  THEN 

begin 

l  :  =  1 1 m e ; 

WRITE. N( 'PRODUCT  • , R3 . ST S I NFO. A DH I N: J , •  EXCEEDS  HOR') 

end; 

FOR  K  :=  (T»l>  TO  L  DO 
BEGIN 

TASSIj N8R8. SYS  INFO .T 1 ,K8  1=  R 8 . S YS I NF 0 . ADM  I N S 
BASSIjN3Ra.SvSlNF0.BAY.Ka  :=  Ra.SYSIN-O. ADMIN 

end; 

Ra.SYSlNr O.TESTHR  R3 . S Y S I NFO . TE STHR - B ; 

L  5=  T  ♦“), 

IF  L  >  T. ME  THEN 

l  :=  i i he  ; 

BASSIGNdr^a.SYSINFO.BAY  ,  L  3  1=  R  3  .  SY  S  I  NF  0  .  A  DH  I N 
END 
ELSE 

BEGIN 

l  :  =  t  *  rs.sysinfq.testhr; 

Ra.SYSlNr  o. finish  ;=  l; 

ppoduct3(\  a. SYsiNFo.ACHiNa. finish  :r  l; 

IF  R3. SYj INFO. COMP  <>  0  then 

PCOMP.P3.SVSINFO.PLiR3.SySINFO.C3MP8  :  =  l; 

IF  L  >  T. HE  THEN 
BEGIN 

l  :=  r imf ; 

WR ITE.N< 'PRODUCT  • , R8 . S YS I NF 0. A DH I N:  3 ,  •  EXCEEDS  HOR') 

end; 

FOR  K  I;  <  T* 1 )  TO  L  DO 
BEGIN 

TASSIjN3Ra.SYSINF0.Tl.K3  l-  R a . S Y S I NFO . ADH 1 N ; 
BASSIjN3Ra.SYSINF0.3AY.K8  :=  R i . S V S I N * 0 . A D M I N 
ENO  ; 

L  :=  T*R« .SYSINFO.TESTHR-C (RS.SYSINFO.DUEDT-I >«DAY*3> ; 

IF  L  )  3  THEN 
BEGIN 

WRITE  ('PRODUCT  • , R 3 . S Y S I NF 0. AD M I N : 3 . »  WILL  BE  *.L:3) 
WR  I  TE.  N  (•  TIME  PERIODS  L  A  T  E*  > ; 

LATEN.SS  ;=  LATENESS  ♦  l; 

LATE  .  =  LATE  ♦  1 
end; 

S3. LINK  ;  =  R3.LINK ; 

DISPOSEU ) ; 

R  :=  s 
END 
end; 
s  :  =  r  ; 

R  :=  R3.LINK 
ENO 

END  ( ‘WHILE  R  C>  . IL*) 

ENO;  ( • I F  ODD*) 
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(•DETERMINE  IF  ALL  PARALLEL  SYSTEM  COMPONENTS  HAVE  COMPLETED  TESTIN6*) 

FOR  I  1  TO  N  3  S  T  S  DO 
BEGIN 
max  :=  i; 
flag  o; 

FOR  J  1  TO  »  DO 
BEGIN 

IF  PC0MP3I.J3  >  Mm  X  then 
MAX  :=  PC0MP3I.JS 

else  if  pcompsi.j.  =  o  then 
f.AG  H  1 
end;  <*f o»  j‘> 

IF  (FLAG  <>  1J  AND  Cr  R  0DUCT3  13  .  FL  AGP  =  0)  AND  (T  =  MAX)  THEN 

(•ALL  COMPONENTS  - OMP L E T E 0 -E N TE R  PARALLEL  SYSTEM  IN  ASSIGNMENT  OUEUE*) 
BEGIN 

P50DJCT3I3.CSCHS  .=  (T  DIV  DAY)  ♦  1} 

PRCDUC  T3I3 . FLAGP  .  =  1 i 
CREATENOOE(  I.NExT.  S 
g;nsert(next,gihEmD) 

END  ( *  I F  c  L  AG  *  > 

end;  (*for  i*) 

T  t  *  l 
eno; 

writeln: 

(•DETERMINE  IF  ANY  PPOOLCT.,  DID  NOT  eEGIN  PROCESSING  IN  HORIZON*) 
for  i  ;i  i  to  n  do 

IF  PR00OCT3I3. START  p  I  THEN 

WRITELN  (•PRODUCT  •  ,  *  R  0  D  UC  T  3  I  J  .  A  DM  I N  !  J  ,  *  GAS  NOT  SCHEDULED*); 

writeln* 

w  D  I  TEL  N  («THE  FOLLOWING  IS  the  ASSIGNMENT  SUMMARY:*)! 
wRIteln; 

WRITELN  (•  ADMIN  SYSNUM  j  A  Y  TE  S  T  OR  1  TE  ST  OR  2  START  FINISH*)! 

FOR  I  :p  1  TO  N  DO 

WITH  PRO0UCT3I3  DO 

WRITELN  (ADMIN:5*SYS.<UM:8tBAY:6tTi:7tT2:9tSTART:6*FINISH:8)* 

«r iteln; 

WRITE  (‘THERE  WERE  '.LATEIj**  LATE  PRODUCTS  WITH  A  MEAN  LATENESS  OF  •)! 

W  R I  *  £  L  N  ((LAT£MFSS  CIV  LAT.KJ,*  TIME  PERIODS*)! 

W  P ITELN  ( • THE  PROPORTION  Or  JOBS  TAROY  WAS  * . ( I  LA T E/N ) • I  00 )  M , •  PERCENT*)! 

FOR  I  :=  1  TO  N  00 

FLOW  :i  FLOw*PRODUCT 313. FINISH- ( (PR0DUCT3I3.CSCHS-1 > *DAY*3 ) ! 

WRITELN  ( • THE  MEAN  uOE  FLC.TIME  WAS  '.(FLOW  DIV  NKA.*  TIME  PERIODS*)! 

WRITELN! 

.RITELN  (• "HESE  ASSIGNMENT^  WERE  made  WITH  SPT*)! 

IF  PRINT  —  1  THEN 
BEGIN 
WRITELN! 

WRITELN  ('TEST  STATION  ASSIGNMENTS  FOR  DAY  1*)! 

WRITELN! 

WRITE  (•  BAY  •)! 

FOR  I  :=  1  TO  UA  DO 
WRITE  (*••)! 

WRITELN! 

FOR  I  ::  1  TO  8  DO 
BEGIN 

UR  ITF  (  I  :  A  ,*  •*  )  ; 

FOR  J  :=  0  TO  OAY-l)  DO 
WRITE  (BASSIGN3I*w3:A)| 

WRITELN! 

WRITELN 
E  NO! 

WRITE  (•  »)! 

FOR  I  Z-  1  TO  UA  DO 
WRITE  (»••)! 

UR  ITELN ; 

wR ITE  ( *  *  )  5 

FOR  I  0  TO  (DAY-1)  U., 

WRITE  ( I :  a  ) ; 

WRITELN! 

WR  ITELN  ( • T IME* : 55  >  5 
wR ITELN ! 

WRITCLN  (‘TESTOR  ASSIGNMENTS  FOR  DAT  1*)! 

WRITELN! 

WRITE  ("TESTOR  •)! 

FOR  I  :=  I  TO  11A  DO 
WRITE  (***)! 

WRITELN! 


1  TO  W  DO 


FOR  I  := 

BES  IN 

write  <  I : r •  •  »•>; 

FOR  J  0  TO  (DAY-U  do 
WRITE  <TASSIGNai,u8:R>; 
written; 

WR ITTL  N 

end; 

WRITE  < •  •  )  ; 

FOR  I  1  TO  114  DO 
WRITE  <•**); 

W  R  I  T  E  L  N ; 

write  (•  •>; 

FOR  I  0  TC  (OAT-1)  0- 
writc  (i:o; 
writeln; 

WRITELN  ( •TI«E* :55) 5 

writeln; 

FOR  K  2  TC  <HOR*WK>  „0 
BEGIN 

WRITELN  ('TEST  STATIjN  ASSIGNMENTS  FOR  DAT  •,<:?) 
writeln; 

WP ITE  ( •  BAT  • >  ; 

FOR  I  :=  1  TO  114  DO 
WRITE  (•••),* 
writeln; 

FOR  I  l =  1  TO  B  DO 
BEGIN 

WRITE  ( I  1  »  .  •  *  •  )i 

FGR  J  <(K-l)«j,Y)  TO  ( *• (DAT-1 ) tK-1 )  DO 

write  ( bassign- i« js:r)  ; 
writeln; 

writer 

end; 

wr ;te  < .  • ) ; 

F  OF  I  :=  1  TO  HR  DJ 
WRITE  <••»>; 

written; 

WR ITE  (  •  » ) ; 

FOR  I  (  (<-l)*DAT)  TC  (K*(DAT-1)*K-1  )  DC 

write  «i:«i; 
written; 

writeln  ( *  T I  WE*  :55)  ; 
writeln; 

WRITELN  ( *  TE  STOP  ASSIGNMENTS  FOR  DAY  * , K : 2  > ! 

written; 

WRITE  (TTESTOR  »>; 

FOR  I  1-  1  TO  114  DO 
WRITE  (••’>; 
written; 

FOP  I  :=  1  TO  W  DO 
BEGIN 

WRITE  <  I : 4 ,  •  •  •  ), 

FOP  J  :=  f  «K-1 > • J„ V  >  TO  (K*(DAT-1)*K-1)  DO 
WRITE  ( TASS IGN. I f js: 4 ) ; 

written; 

WRITTEN 
END  ; 

WRITE  (*  »>: 

FOR  I  :=  1  to  114  DC 

write  <•••>; 

WR  I  T  E  L  N  ; 

WR I TE  <  •  » ) ; 

FOR  I  :=  ((K-I)»D<T>  TO  IK*IDAY-1)»K-1>  DO 

WRITE  (  i:*l  ; 
written; 

WRITELN  < *TI«C':55) ; 

WRITELN 

END 

ENL 
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Goodness-of-Fit  Tests  for  Random  Number  Generators 


Chi-Squared  Goodness  of  Fit  Test  for  Uniform  Distribution: 
The  following  procedure  was  used  to  generate  uniform 

random  numbers,  where  IX  is  a  uniform  random  number  seed: 

*  -  a  t  ;  v  ;  v  :  :  i  *  -  ,  •  a  v  ytl  t  j  >  \  ; 

( ■  i  ,  j  .  ;  ■  .  ;  *  *  .  i .  ■  > 

: *  .>i  : 

I  y  :  I <  ■*  -  ■  > T  «  i 

I »  .  r  <  r  1  “  \ 

.  :  =  :  y  ♦  : :  ♦  ? 

yi  i  :  r  i  y  ; 

Y  L  ;  =  y  L  *  l  t  .  ' 

'  .  » 

The  procedure  was  called  200  times  with  the  results  as 
indicated  in  the  following  table: 


Interval 

0.0-.  I 

.  I-.2 

.2-. 3 

.3-. 4 

.4-. 5 

.5-. 6 

.6-.  7 

.  7 -.8 

.8-. 9 

.9-1.0 

Obs  Freq 

28 

22 

17 

17 

19 

21 

22 

14 

20 

20 

Exp  Freq 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

The  hypotheses  and  corresponding  test  statistic  are: 

H  :  the  underlying  distribution  is  uniform 
o  -1 

H^:  the  underlying  distribution  is  not  uniform 

2  =  V  (OBS  -  EXP)2 
X  Lj  r.XP 

2  -> 

Rejection  Region:  If  x  x~  j.  ^  Reject  Hq 

2  -> 

Acceptance  Region:  If  x  ±  Accept 


'  f  Xoc,k- 1 -m 


X“ 


Xa,k-1 


< 


< 


Withhold  judgment 
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where  k=  number  of  cells  and  m=  number  of  estimated 
parameters-  In  this  case  k=  10  and  m=  1  (p) . 

With  a  95%  confidence  interval  a  =  .05  and 


Xa,k-1  =  \05,9  =  16'919 

Xa,k-l-m  =  X. 05 , 8  =  15-50' 
For  the  specific  problem. 


2  =  y  lOBS-EXP)2  = 

X  EXp 

2  2 

Since  X  <  Xa  x-i-m  the  hypothesis  is  accepted  and  the 

uniform  distribution  provides  a  good  fit  to  the  data. 
Therefore,  the  procedure  is  an  acceptable  generator  of 
uniform  random  numbers. 


Chi  Squared  Goodness  of  Fit  Test  for  Normal  Distribution: 

The  following  procedure  was  used  to  generate  random 
numbers  from  a  normal  distribution,  where  J  is  a  uniform 
random  number  seed  and  R  is  a  uniform  random  number: 


t'X-.Oi  >J;-  .  "•  *  <  \!  '>■  i :  >  i_  )  ; 

LA  •  L 


</  >'■ 

I  it  :  ’ll.,:'; 

T  1  T  i  f  i  •  t  :  i  i  ; 

m  , 

t  -  .  ; 

II  <  f  T  ' 

1  '  ■  I  * 

~  •  i  :  ■  ;  :  >  i 

!■_  ,1  . 

IJ  .  r  :  ’  < 

)  :  :  _ 

T  .>  r  :  - 


r  it  ♦ 


The  procedure  was  called  210  times  with  the  results  as 
indicated  in  the  following  table: 
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Cell 

(-00,-1.07) 

(-1.07 ,-.57) 

(-.57 ,-.18) 

(-.18,. 18) 

(.18, .57) 

(.57,1.07) 

(1.07,00) 

Obs  Freq 

35 

25 

24 

34 

34 

30 

28 

Exp  Freq 

30 

30 

30 

30 

30 

30 

30 

The  hypotheses  and  corresponding  test  statistic  are: 


Hq:  the  underlying  distribution  is  normal 

H  :  the  underlying  distribution  is  not  normal 

d 

2  y*  (OBS  -  EXP)2 
X  “  ^  EXP 

2  2 

Rejection  Region:  If  x  >  x  ,  ,  Reject  H 
J  ^  A  -  a,k  - 1  J  o 


2  ° 

Acceptance  Region:  If  x  li  k  1  m  AccePt 


If  < y2  <  y2  ,  Withhold  judgment 

xa,k-l,m  *  xa,k-l 

,  where  k=  number  of  cells  and  m=  number  of  estimated 
parameters.  In  this  case  k=  7  and  m=  2  (u,  o) • 

With  a  95%  confidence  interval  a=  .05  and 

xl  ,  i  =  X2„  .  =  12-592 


\*,k-l  A. 05 , 6 


xa,k-l-m  X.05,4 


9.488 


For  the  specific  problem 


4.0667 


Since  <  X?  the  hypothesis  is  accepted  and  the 

®,k-l-m 

normal  distribution  provides  a  good  fit  to  the  data. 
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Therefore,  the  procedure  is  an  acceptable  generator  of 
normal  random  numbers. 


Chi-Squared  Goodness-of-Fit  Test  for  Poisson  Distribution: 

The  following  procedure  was  used  to  generate  random 
numbers  from  a  poisson  distribution,  where  J  is  a  uniform 
random  number  seed  and  Lambda  is  the  mean,  X  ,  for  the  dis¬ 


tribution: 


p«  ; c j  1  j:  '■  (  m  \  ->  j:  l -  *-l  *  v;  *  ■» : :  \  rr*  z  '  ° i  ; 

l  a  i  :  l 
r>  ! 

V  A"' 

i  1 y 

l  :  I  :  t  -  .  •  >.  ; 

!tf  ", ! 

y  :  -  i  ; 
f-  :  -  : 

o  :i)  ir  t  m  ( j  ,l  »  '  •  > 
j  :  -  l  ; 
y  :  -  y  ; 

I  r  y  >  (.<’(•.  i‘i  i  )  t  i ■ 

i:,!  j 

■>  ;  r  I  >  ♦  ’  : 

■j ;  ■  T  / 

:  'i  > 

1  )  i  >•  i 


The  procedure  was  called  200  times  with  X~  3.  The  results 
are  as  indicated  in  the  following  table. 


143 


Xi 

Obs  Freq 

■HHH 

Exd  Frea 

0 

15 

.  0498 

9.96 

1 

30 

.  1494 

29.88 

2 

39 

.  2240 

44 . 8 

3 

43 

.  2240 

44.8 

4 

36 

.  1680 

33.6 

5 

22 

.  1008 

20.16 

6 

7 

lS 

.0504 

10.08 

8 

0 

9 

0 

>  6 

.  0334 

6.68 

10 

0 

1 1 

The  probabilities  associated  with  various  values  of  x 
are  obtained  using  the  probability  mass  function  for  the 
poisson  distribution 


p  (X)  = 


< 


-  X,  x 

c _ X_ 

x.' 


X=  0 ,  1 ,  2, 


0  otherwise 


The  hypotheses  and  corresponding  test  statistic  are: 

H  :  the  random  variable  is  Poisson  distributed 
o 

H  :  the  random  variable  is  not  Poisson  distributed 
a 

2  ^  (OBS  -  EXP)2 

X  ^  EXP 


2  2 

Rejection  Region:  If  y  >  y  Reject  H 

x  -  xa,k-l  o 


1  "> 

Acceptance  Region:  If  x*-  <  X*"  Accept  H 

—  a,k-l-m  o 


->  ~y 

I  f  X"  <  y-  <  Y“  withhold  judgement 

a,k-l-m  x  xa,k-l 

where  k=  number  of  cells  and  m=  number  of  estimated  parame¬ 
ters.  In  this  case  k=7  and  m=l  (a  ).  With  a  95%  confidence 
interval 
a  =  .05  and 


xa,k-l 

2 

Xa,k-l-m 
For  the  specific  problem 

-  3.898 

.  2  2 

Since  x  <  X  a  m  the  hypothesis  is  accepted  and  the 
Poisson  distribution  provides  a  good  fit  to  the  data. 
Therefore,  the  procedure  is  an  acceptable  generator  of 
Poisson  random  numbers. 

Chi  Squared  Goodness-of-Fit  Test  for  Gamma  Distribution : 
The  following  procedure  was  used  to  generate  random 


*:05,6 =  i2-592 


*:0s)5  =  n-o7° 


numbers  from  a  Gamma  distribution,  where  IX  is  a  uniform 
random  number  seed  and  Lam  and  Neta  are  the  parameters  of 
the  distribution: 
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The  procedure  was  called  200  times  with  the  parameters 
Lam=  1.0  and  Neta=  4.0  with  the  results  as  indicated  in  the 
following  table: 


m 

Obs  Freq  P(Xi) 

Exp  Freq 

0 

3 

.  019 

3.8 

1 

24 

.124 

24.8 

2 

43 

.21 

42 

3 

43 

.214 

42.8 

4 

36 

.168 

33.6 

5 

23 

.114 

22.8 

6 

13 

.069 

13.8 

7 

7 

.04 

8 

8 

2\ 

9 

3 

10 

1 

>8  .042 

8.4 

11 

0 

12 

1 

12 

The  probabilities  associated  with  various  values  of  X 
are  obtained  using  the  probability  mass  function  for  the 
gamma  distribution. 


P  (X)  =  N 


nr  r-1  -0x 
Ox  e 

T(r) 


,  if  X  >  0 


,  if  X  <  0 


The  hypothesis  and  corresponding  test  statistic  are: 

H  :  the  random  variable  is  gamma  distributed 
o  3 

H  ;  the  random  variable  is  not  gamma  distributed 

2  =  v  (obs  -  r.xr)2 
x  ^  nxr 


]  47 


'i  2 

Rejection  Region:  If  x"  ^  j  Reject  Hq 

2  2 

Acceptance  Region:  If  x  j  xa  ^  1_m  Accept  Hq 
2  o  -> 

If  x;  <  <  xa  k-l  withhold  judgment 

where  k=  number  of  cells  and  m=  number  of  estimated  parame- 
ters.  In  this  case  k=9  and  m=2  (0 ,r),  where  0=  Lam  and 
Neta.  With  a  96*  confidence  interval, a  =  .05  and 


Xcx,k-1 


X. 05 , 8 


15.507 


o  2 

Xa,k-l-m  ~  x. 05,6 


12.592 


For  the  specific  problem 


2 

X 


_  y  (qbs  -  Expy 

=  ^  EXP 


.583 


~> 

Since  x~  < 


£. 


xa,k- 1-m 


the  hypothesis  is  accepted. 
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